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THE  EFFECT  OF  PHYSICAL  AGING  ON  THE 
CREEP  RESPONSE  OF  A  THERMOPLASTIC  COMPOSITE. 


Robert  L.  Hastie  Jr.,  Major,  USAF 

Doctor  of  Philosophy,  1991,  240  pages 
Department  of  Engineering  Science  and  Mechanics 
Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  Virginia 

(ABSTRACT) 


‘ '  The  effect  of  thermoreversible  physical  aging  on  the  linear  viscoelastic  creep  prop¬ 
erties  of  a  thermoplastic  composite  was  investigated.  Radel  X/IM7,  an  amorphous  com¬ 
posite  material  considered  for  use  in  the  next  generation  high  speed  transport  aircraft,  was 
studied.  The  operating  environment  for  the  aircraft  material  will  be  near  188°C  (370°F) 
with  a  service  life  in  excess  of  60,000  hours  at  temperature.  Accurate  predictions  of  the 
viscoelastic  properties  of  the  material  are  essential  to  insure  that  design  strength  and 
stiffness  requirements  are  met  for  the  entire  service  life. 

The  effect  of  physical  aging  on  the  creep  response  was  studied  using  momentary 
tensile  creep  tests  conducted  at  increasing  aging  times  following  a  rapid  quench  from 
above  the  glass  transition  temperature  (Tg)  to  a  sub-Tg  aging  temperature,  ^the  agjng 
time  increased,  the  creep  response  of  the  material  significantly  decreased,  ^e  tensile 
creep  compliance  data  for  each  aging  time  were  fit  to  the  empirical  equation  for  the  creep 
compliance  D(t); 


D(t)  =  Doe 


where  D^,  r^,  and  m  are  fitting  parameters  determined  using  a  nonlinear  fitting  program 

based  on  the  Levenberg-Marquardt  finite  difference  algorithm.  The  short-term  creep 
compliance  curves,  obtained  at  various  aging  times,  were  then  shifted  to  form  a  momentary 
master  compliance  curve.  The  double-logarithmic  aging  shift  rate  p  and  its  dependence 
on  sub-Tg  aging  temperature  were  determined.  The  aging  characterization  process  was 
conducted  on  unidirectional  specimens  with  0,  90,  and  45  degree  fiber  direction  orienta¬ 
tions.  This  permitted  the  calculation  of  the  complete  principal  compliance  matrix  for  the 
composite  material. 

The  effect  of  physical  aging  becomes  more  apparent  during  long-term  tests  when 
creep  and  aging  occur  simultaneously.  This  results  in  a  gradual  stiffening  and  decrease 
in  the  creep  response  with  increased  time.  Predictions  based  solely  on  the  Time- 
Temperature  Superposition  Principle  would  significantly  over-predict  the  creep  response 
if  physical  aging  effects  were  ignored.  Theoretical  predictions  for  long-term  creep 
compliance  were  made  using  an  effective  time  theory  and  compared  to  long-term  exper¬ 
imental  data  for  each  fiber  orientation.  Finally,  experimental  results  of  a  long-term  test  of 
a  30  degree  fiber  angle  orientation  specimen  were  compared  to  theoretical  predictions 
obtained  by  transforming  the  principal  compliance  matrix  to  the  30  degree  orientation. 
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The  effect  of  thermoreversible  physical  aging  on  the  linear 
viscoelastic  creep  properties  of  a  thermoplastic  composite  was 
investigated.  Radel  X/IM7,  an  amorphous  composite  material 
considered  for  use  in  the  next  generation  high  speed  transport 
aircraft,  was  studied.  The  operating  environment  for  the 
aircraft  material  will  be  near  188°C  (370°F)  with  a  service 
life  in  excess  of  60,000  hours  at  temperature.  Accurate 
predictions  of  the  viscoelastic  properties  of  the  material  are 
essential  to  insure  that  design  strength  and  stiffness 
requirements  are  met  for  the  entire  service  life. 

The  effect  of  physical  aging  on  the  creep  response  was 
studied  using  momentary  tensile  creep  tests  conducted  at 
increasing  aging  times  following  a  rapid  quench  from  above  the 
glass  transition  temperature  (Tg)  to  a  sub-Tg  aging  temperature. 
As  the  aging  time  increased,  the  creep  response  of  the  material 
significantly  decreased.  The  tensile  creep  compliance  data 
for  each  aging  time  were  fit  to  the  empirical  equation  for  the 
creep  compliance  D(t): 


where  D^,  to,  and  m  are  fitting  parameters  determined  using  a 

nonlinear  fitting  program  based  on  the  Levenberg-Marquardt 
finite  difference  algorithm.  The  short-term  creep  compliance 
curves,  obtained  at  various  aging  times,  were  then  shifted  to 
form  a  momentary  master  compliance  curve.  The  double- 
logarithmic  aging  shift  rate  n  and  its  dependence  on  sub-T, 
aging  temperature  were  determined.  The  aging  characterization 
process  was  conducted  on  unidirectional  specimens  with  0,  90, 
and  45  degree  fiber  direction  orientations.  This  permitted 
the  calculation  of  the  complete  principal  compliance  matrix 
for  the  composite  material.  The  effect  of  physical  aging 
becomes  more  apparent  during  long-term  tests  when  creep  and 
aging  occur  simultaneously.  This  results  in  a  gradual 
stiffening  and  decrease  in  the  creep  response  with  increased 
time.  Predictions  based  solely  on  the  Time-Temperature 
Superposition  Principle  would  significantly  over-predict  the 
creep  response  if  physical  aging  effects  were  ignored. 
Theoretical  predictions  for  long-term  creep  compliance  were 
made  using  an  effective  time  theory  and  compared  to  long-term 
experimental  data  for  each  fiber  orientation.  Finally, 
experimental  results  of  a  long-term  test  of  a  30  degree  fiber 
angle  orientation  specimen  were  compared  to  theoretical  pre¬ 
dictions  obtained  by  transforming  the  principal  compliance 
matrix  to  the  30  degree  orientation. 


ACKNOWLEDGEMENTS 


The  author  gratefully  acknowledges  the  numerous  individuals 
who  have  made  this  study  possible.  Each  of  you  have  made 
your  own  special  contribution  to  help  me  over  this  hurdle  in 
my  life. 

I  especially  thank  Colonel  Cary  Fisher,  Head  of  the 
Department  cf  Engineering  Science  and  Mechanics,  USAF  Academy, 
for  the  opportunity  to  continue  my  education  and  for  expressing 
his  confidence  in  me  and  my  abilities.  Your  respected  voice 
of  encouragement  reached  a  long  way  to  boost  my  morale. 

My  deepest  thanks  go  to  my  Committee  chairman  Dr.  Don 
Morris  for  his  advice,  patience,  and  never  ending  enthusiasm. 
You  provided  the  shining  light  for  all  those  dark  moments  that 
occur  during  experimental  work. 

I'd  also  like  to  thank  Dr.  E.G.  Henneke  II,  Dr.  D.T.  Mook, 
Dr.  W.W.  Stinchcomb,  and  Dr.  E.R.  Johnson  for  serving  as 
committee  members.  It  has  been  a  pleasure  knowing  and  working 
with  each  of  you. 

I  thank  Dr.  Charles  Harris  and  Dr.  Thomas  Gates  from  NASA 
Langley  Research  Center  for  their  patience  during  the  initial 
stages  of  the  project  and  for  providing  the  composite  material. 

A  special  thanks  goes  to  everyone  who  has  helped  me  with 
the  details  of  conducting  my  experiments.  Bob  Simmons  provided 
the  data  acquisition  system.  Dr.  D.  Dillard,  Yeou  Shin  Chang, 


IV 


and  Danny  Reed  gave  their  advice  and  shared  the  Adhesive  Science 
lab  equipment.  Dr.  Hyer  donated  the  NASA  oven  controller  and 
data  acquisition  board  for  use  in  this  study.  Rob  Swain,  Jack 
Lesko,  Greg  Carman,  and  Ahmad  Razvan  provided  friendship  and 
long  term  use  of  MRG  equipment.  George  Lough  provided  the 
strain  gage  conditioning  amplifiers.  Bob  Davis,  Bill  Shaver, 
and  Archie  Montgomery  squeezed  in  my  no  notice  machining  needs. 
Randy  Waldan  did  the  welding  of  the  numerous  thermocouples. 
Min-Chung  Li  helped  determine  the  fiber  volume  fraction  of  the 
composite.  Paul  Vail,  from  the  Chemistry  department,  conducted 
the  DSC  and  TGA  scans.  Pedro  Morales  shared  the  office  in  309 
Femoyer  with  me.  Robert  Canfield  and  Mike  Smith  provided 
timely  comments  on  dissertation  text. 

To  my  wife,  Vicki,  and  daughter  Jessica,  your  love, 
understanding,  and  little  hugs  gave  me  the  will  to  continue. 
I'll  always  be  there  for  you,  as  you  were  here  for  me.  We  can 
make  it  through  anything  as  long  as  we  have  each  other. 

To  my  mother,  Joan,  and  my  sister  Sheryl,  thanks  for  all 
your  love  and  support  during  my  periods  of  shaken  confidence. 
It  was  enjoyable  living  close  enough  so  we  could  have  weekend 
visits  together. 

Finally,  I'd  like  to  dedicate  this  work  to  my  deceased 
father.  Bob  Hastie.  He  sacrificed  everything  during  his  life 
so  that  my  life  would  be  better.  No  deeds  or  acts  can  ever 
fully  express  my  gratitude... 


V 


TABLE  OF  CONTENTS 


1  INTRODUCTION  .  1 

1 . 1  Literature  Review  .  5 

1.1.1  Creep  Characterization  at  VPI  &  SU  . 5 

1.1.2  Physical  Aging  Effects  .  12 

1.2  Overview  of  Proposed  Study  .  16 

2  BACKGROUND  INFORMATION  .  19 

2.1  Viscoelastic  Materials  .  19 

2.1.1  Common  Tests  . 20 

2.1.2  Linear  Viscoelasticity  .  22 

2. 1.2.1  Definition  . 22 

2. 1.2.2  Boltzmann's  Superposition  Principle  .  23 

2. 1.2. 3  Mechanical  Models  . 25 

2. 1.2. 4  Power  Law  Model  .  29 

2. 1.2.5  Struik  Model  .  30 

2. 1.2. 6  Time-Temperature  Superposition 

Principle  . 31 

2.1.3  Nonlinear  Viscoelasticity  .  35 

2. 1.3.1  Findley  Model  . 36 

2. 1.3.2  Schapery  Model  . 37 

2.2  Physical  Aging  Effects  .  41 

2.2.1  Basic  Polymer  Concepts  .  42 

2.2.2  Physical  Aging  Definition  .  45 

2.2.3  Effect  of  Physical  Aging  on  Material 

Properties  . 47 


2.2.4  Effect  of  Physical  Aging  on  Creep 

Compliance  . 4  8 

2.2.5  Characterizing  and  Predicting  the  Effect 

of  Physical  Aging  . 51 

2. 2. 5.1  Momentary  Compliance  Curves  . 52 

2. 2. 5. 2  Aging  Shift  Rate  Determination  .  53 

2. 2. 5. 3  TTSP  Master  Momentary  Compliance 

Curve  . 57 

2. 2. 5. 4  Effective  Time  Theory  . 59 

2.3  Composite  Material  Compliance  Relationships  .  63 

3  EXPERIMENTAL  METHODS  .  74 

3.1  Testing  Equipment  .  74 

3.2  Calibration  Techniques  .  75 

3.2.1  Oven  Temperature  . 75 

3.2.2  Quench  Rate  Determination  . 78 

3.2.3  Specimen  Suress  . 80 

3.3  Specimen  Preparation  .  84 

3.4  Strain  Gage  Installation  .  86 

3.5  Data  Collection  .  92 

3.6  Aging  Creep  Test  Procedures  .  96 

3.7  Data  Reduction  Procedures  .  98 

4  CHARACTERIZATION  RESULTS  AND  DISCUSSION  .  101 

4.1  Room  Temperature  Material  Response  .  101 

4.2  90  Degree  Specimen  Tests,  S22  Determination  ..  103 

4.2.1  Linearity  and  Sequencing  Effects  .  104 

4.2.2  Aging  Shift  Rate  Determination  .  105 

4.2.3  TTSP  Master  Momentary  Compliance  Curve  ..  106 


Vll 


4.3  45  Degree  Specimen  Tests,  '^66  Determination  . .  108 

4.3.1  Linearity  and  Sequencing  Effects  .  108 

4.3.2  Aging  Shift  Rate  Determination  .  109 

4.3.3  TTSP  Master  Momentary  Compliance  Curve  ..  110 

4.4  0  Degree  Specimen  Tests  .  112 

4.4.1  5,1  Determination  . 113 

4.4.2  5,2  Determination  . 114 

4.5  Temperature  Dependence  of  Aging  Shift  Rate  ..  115 

4.6  Temperature  Dependence  of  Shift  Factors  .  116 

4.7  Quench  Rate  Effect  . 117 

5  LONG  TERM  CREEP  RESULTS  AND  DISCUSSION  . 120 

5.1  90  Degree  Specimen  Test,  S22  Comparison  . 120 

5.2  45  Degree  Specimen  Tests,  Comparison  . 122 

5.3  30  Degree  Specimen  Test,  S„  Comparison  .  124 

5.4  Overall  S„  Comparison  .  125 

6  SUMMARY,  RECOMMENDATIONS,  AND  CONCLUSIONS  .  127 

6.1  Summary  .  127 

6.2  Recommendations  .  12  9 

6.3  Conclusions  .  132 

REFERENCES  .  133 

TABLES  . 141 

FIGURES  .  146 

VITA  . 225 


viii 


TABLE  OF  TABLES 


Table  Page 

4.1  Comparison  of  Commercial  Data  of  Radel  X/T650-42 

and  Radel  X/IM7  Composite  Materials .  142 

4.2  Comparison  of  Commercial  T650-42  and  IM7  Fiber 

Properties .  143 

4.3  Shift  Factors  Required  for  S22  Master  Momentary 

Compliance  Curves  Referenced  to  /,  =  81  hours  ...  144 

4.4  Shift  Factors  Required  for  ^66  Master  Momentary 

Compliance  Curves  Referenced  to  t,  =  81  hours  .  .  .  145 


IX 


TABLE  OF  FIGURES 


Figure  Page 

2.1  Stress  and  Strain  Response  for  a  Creep/Recovery 

Test .  147 

2.2  Stress  and  Strain  Response  for  a  Stress  Relaxation 

Test .  14  8 

2.3  Stress  and  Strain  Response  for  a  Constant  Strain 

Rate  Test .  149 

2.4  Graphical  Representation  of  the  Requirements  for 

a  Linear  Viscoelastic  Material .  150 

2.5  Stress  History  Approximation  Using  a  series  of 

Step  Stress  Inputs  .  151 

2.6  Representation  of  Maxwell  and  Kelvin  Mechanical 

Analog  Models .  152 

2.7  Representation  of  Generalized  Kelvin  Model .  153 

2.8  Time-Temperature  Superposition  Principle  Illus¬ 
tration  .  154 

2.9  Visualization  of  Two-Dimensional  Hard-Sphere 

Model  of  Free  Volume  [75] .  155 

2.10  Illustration  of  the  Qualitative  Free  Volume 

Concept .  156 

2.11  Typical  Compliance  Response  of  Thermoplastic  and 

Thermoset  Materials .  157 

2.12  The  Origin  of  Aging  Explained  from  Free  Volume 

Concepts  [59] .  158 

2.13  Tensile  Creep  Compliance  Curves  for  PVC  Quenched 

from  90  to  20°C  [59]  .  159 

2.14  Calculation  of  ]J.  from  the  Shifting  of  the  Creep 
Curves  of  Rigid  PVC  Quenched  from  90°C  to  0,  20, 

and  40°C  [59]  .  160 

2.15  The  Aging  Shift  Rate,  M-,  vs.  Temperature  for 

Various  Polymers  [59] .  161 


X 


2.16  Illustration  of  the  Sequence  of  Creep  and  Recovery 

Tests  for  Determining  Aging  Effects  [59]  .  162 

2.17  Illustration  of  Aging  Shift  Rate  Determination 

by  Shifting  Momentary  Creep  Compliance  Curves  .  .  163 

2.18  Illustration  of  Effective  Time  Variables .  164 

2.19  Coordinate  Systems  Used  to  Define  Material 

Properties .  165 

2.20  Strain  Gage  Rosette  Placement  on  45  degree 

Specimen .  166 

3.1  Illustration  of  Equipment  Used  for  Creep  Testing.  167 

3.2  Illustration  of  Specimen  Attachment  to  Load  Frame .  168 

3.3  Stabilization  Temperature  vs.  Quench  Time .  169 

3.4  Specimen  Temperature  vs.  Time  Data  Used  for  Quench 

Rate  Calculation .  170 

3.5  Calibration  of  Output  Voltage  vs.  Specimen  Load.  171 

3.6  Calibration  of  Specimen  Load  vs.  Applied  Load  on 

Arm .  172 

3.7  Hot  Press  Cure  Cycle .  173 

3.8  DSC  Trace  of  the  Radel  X/IM7  Composite .  174 

3.9  TGA  Trace  of  the  Radel  X/IM7  Composite .  175 

3.10  Illustration  of  Full  Bridge  Wheastone  Circuit,.  176 

3.11  Strain  Gage  Curing  Cycles .  177 

3.12  Initial  Zero  Drift  Measurements . 178 

3.13  Zero  Drift  Measurements  after  Resoldering .  179 

3.14  Zero  Drift  Measurements  after  Remounting  Dummy 

Specimen .  180 

3.15  Comparison  of  Line  Noise  Between  the  Vishay  2100 

and  2210  Signal  Conditioner  Amplifiers .  181 

3.16  Linear  Fit  to  Recovery  Voltages .  182 

4.1  Stress  Strain  Curve  for  Transverse  Direction...  183 


XI 


4.2  Stress  Strain  Curve  for  Fiber  Direction .  184 

4.3  Linearity  Check  of  Transverse  Direction  Com¬ 
pliance  at  a  Constant  Temperature  of  392°F  (200°C) 

and  /«  =  3  hours .  185 

4.4  Check  of  Test  Sequencing  Effect  on  Transverse 
Direction  Compliance  at  a  Constant  Temperature 

of  338°F  (170°C)  and  t,  =  9  hours .  186 

4.5  Transverse  Direction  Compliance  for  a  Constant 

Temperature  of  338°F  (170°C)  at  Various  Aging 

Times .  187 

4.6  Transverse  Direction  Compliance  for  a  Constant 

Temperature  of  365°F  (185°C)  at  Various  Aging 

Times .  188 

4.7  Transverse  Direction  Compliance  for  a  Constant 

Temperature  of  383°F  (195°C)  at  Various  Aging 

Times .  189 

4.8  Transverse  Direction  Compliance  for  a  Constant 

Temperature  of  392°F  {200°C)  at  Various  Aging 

Times .  190 

4.9  Transverse  Direction  Compliance  for  a  Constant 

Temperature  of  399. 2°F  (204°C)  at  Various  Aging 
Times .  191 

4.10  Transverse  Direction  Master  Momentary  Compliance 
Curves  at  Various  Temperatures  Shifted  to  the  81 

hour  Aging  Time .  192 

4.11  Aging  Shift  Rate  for  the  90  degree  Test  Sequence 

at  Various  Temperatures .  193 

4.12  Transverse  Direction  TTSP  Master  Momentary  Com¬ 
pliance  Curve  Fit  with  =  0.4167  and  Referenced 

to  T  =  338°F  (170°C),  t,  =  9  hours .  194 

4.13  Transverse  Direction  TTSP  Master  Momentary  Com¬ 
pliance  Curve  Graphically  Fit  and  Referenced  to 

T  =  338°F  (170°C),  r.  =  9  hours .  195 

4.14  Transverse  Direction  TTSP  Master  Momentary  Com¬ 
pliance  Curve  Graphically  Fit  and  Referenced  to 

T  =  338°F  (170°C),  t,  =  81  hours .  196 


xii 


4.15  Comparison  of  Transverse  Direction  TTSP  Master 
Momentary  Compliance  Curves  Referenced  to  T  = 

338°F  (170°C),  r.  =  9  hours .  197 

4.16  Linearity  Check  of  Shear  Direction  Compliance  at 

a  Constant  Temperature  of  365°F  (185°C)  and  /,  =  3 
hours .  198 

4.17  Check  of  Test  sequencing  Effect  on  Shear  Direction 

Compliance .  199 

4.18  Shear  Direction  Compliance  for  a  Constant  Tem¬ 
perature  of  338°F  (170°C)  at  Various  Aging  Times.  200 

4.19  Shear  Direction  Compliance  for  a  Constant  Tem¬ 
perature  of  365°F  (185°C)  at  Various  Aging  Times.  201 

4.20  Shear  Direction  Compliance  for  a  Constant  Tem¬ 
perature  of  383°F  {195°C)  at  Various  Aging  Times.  202 

4.21  Shear  Direction  Master  Momentary  Compliance 
Curves  at  Various  Temperatures  Shifted  to  the  81 

hour  Aging  Time .  203 

4.22  Aging  Shift  Rate  for  the  45  degree  Test  Sequence 

at  Various  Temperatures .  204 

4.23  Shear  Direction  TTSP  Master  Momentary  Compliance 
Curve  Fit  with  =  0.4564  and  Referenced  to  T 

=  338°F  (170°C),  t,  =  9  hours .  205 

4.24  Shear  Direction  TTSP  Master  Momentary  Compliance 

Curve  Graphically  Fit  and  Referenced  to  T  =  338°F 
(170°C),  f,  =  9  hours .  206 

4.25  Shear  Direction  TTSP  Master  Momentary  Compliance 

Curve  Graphically  Fit  and  Referenced  to  T  =  338°F 
(170°C),  t,  =  81  hours .  207 

4.26  Comparison  of  Shear  Direction  TTSP  Master 
Momentary  Compliance  Curves  Referenced  to  T  = 

338°F  (170°C),  t.  =  9  hours .  208 

4.27  Linearity  Check  of  Fiber  Direction  Compliance  at 

a  Constant  Temperature  of  338°F  (170°C)  and  /,  =  3 
hours .  209 

xiii 


4.28  Fiber  Direction  Compliance  for  a  Constant  Tem¬ 
perature  of  338°F  (170°C)  at  Various  Aging  Times.  210 

4.29  Fiber  Direction  Compliance  for  a  Constant  Tem¬ 
perature  of  399  .  2°F  {204°C)  at  Various  Aging  Times .  211 

4.30  Linearity  Check  of  Fiber/Transverse  Coupling 

Compliance  at  a  Constant  Temperature  of  338°F 
(170°C)  and  =  3  hours .  212 

4.31  Fiber/Transverse  Coupling  Compliance  for  a  Con¬ 

stant  Temperature  of  338°F  (170^0  at  Various  Aging 
Times .  213 

4.32  Fiber/Transverse  Coupling  Compliance  for  a  Con¬ 

stant  Temperature  of  399. 2°F  (204°C)  at  Various 
Aging  Times .  214 

4.33  Temperature  Dependence  of  the  Aging  Shift  Rate.  215 

4 . 34  Temperature  Dependence  of  the  Transverse  Direc¬ 
tion  TTSP  Horizontal  Shift  Factors .  216 

4.35  Temperature  Dependence  of  the  Shear  Direction 

TTSP  Horizontal  Shift  Factors .  217 

4.36  Comparison  of  Air  Quenched  and  Oven  Cooled 

Transverse  Direction  Master  Momentary  Compliance 
Curves  at  383°F  (195°C)  with  r,  =  81  hours .  218 

4.37  Comparison  of  Air  Quenched  and  Oven  Cooled  Shear 

Direction  Master  Momentary  Co.apliance  Curves  at 
383°F  (195°C)  with  =  81  hours .  219 

5.1  Transverse  Direction  Long-term  Compliance  Pre¬ 

dictions  and  Test  Data  at  338°F  (170°C)  and  /,  =  9 
hours .  220 

5.2  Shear  Direction  Long-term  Compliance  Predictions 

and  Test  Data  at  338°F  (170°C)  and  t,  =  9  hours..  221 

5.3  Shear  Direction  Long-term  Compliance  Predictions 

and  Test  Data  at  338°F  (170°C)  and  t,  =  3  hours..  222 

5.4  Load  Direction  Long-term  Compliance  Predictions 
for  the  30  degree  Specimen  at  338°F  (170°C)  and  t, 

=  9  hours .  223 


xiv 


5.5  Load  Direction  Long-term  Compliance  Predictions 
and  Test  Data  at  338°F  (170°C)  and  /,  =  9  hours,. 


224 


XV 


1  INTRODUCTION 


The  use  of  composite  materials  has  dramatically  increased 
during  recent  years.  Today,  numerous  industries  use  composite 
materials  to  fabricate  products.  Composite  materials  are  used 
in  products  ranging  in  diversity  from  aircraft  structures, 
automobile  parts,  boat  hulls,  submarines,  to  artificial  hip 
joints,  sporting  equipment,  and  body  armor.  The  potential 
advantages  of  composite  materials  over  conventional  materials 
are  numerous.  Two  important  advantages,  especially  in  weight 
sensitive  applications,  are  the  improved  strength-to-weight 
and  strength-to-stif fness  ratios  composites  materials  have 
over  metal  alloys.  The  aerospace  industry  initially  capitalized 
on  these  material  advantages  by  building  replacement  parts  for 
aircraft  using  composite  materials.  The  replacement  parts 
were  lighter  in  weight  and  contained  fewer  pieces  than  their 
metal  counterparts.  More  importantly,  the  aerospace  industry 
gained  valuable  experience  working  with  composite  materials 
during  the  design,  fabrication,  and  life  validation  of  the 
replacement  parts.  This  positive  experience  with  composites 
helped  convince  aerospace  industry  leaders  to  start  utilizing 
composite  materials  in  new  aircraft  designs. 
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Another  advantage  of  composites  is  the  ability  to 
"tailor-design"  the  directional  properties  of  the  material. 
Single-ply  direction  laminated  fiber  reinforced  plastic  com¬ 
posites  are  highly  orthotropic  in  stiffness  and  strength.  The 
properties  of  a  composite  material  can  vary  from 
quasi-isotropic,  orthotropic,  to  anisotropic,  depending  upon 
the  ply  orientations  used  in  the  laminate.  Conventional  metal 
alloys  are  usually  considered  isotropic  in  both  stiffness  and 
strength.  The  X-29  forward  swept  wing  aircraft  is  an  example 
of  how  the  orthotropic  nature  of  composites  was  utilized  in 
the  aircraft  design  to  solve  an  aerodynamic  instability  problem. 
The  aircraft  wing  structure  was  fabricated  using  a  selective 
layering  of  graphite-epoxy  plies  of  various  thickness  and  angle 
orientations.  The  resulting  forward  swept  wing  twists  the 
leading  edge  down  when  the  wing  bends  upward.  This  decreases 
the  angle  of  attack  of  the  wing  and  keeps  the  structure 
aerodynamically  stable. 

The  military,  in  its  quest  for  state  of  the  art  aerospace 
equipment,  is  a  primary  driver  and  funding  source  for  composite 
material  research  and  development.  Recent  aircraft  development 
programs,  like  the  F-117  stealth  fighter  and  B-2  bomber, 
extensively  utilized  composite  materials.  The  advanced  tac¬ 
tical  fighter  (ATF) ,  currently  in  the  demonstrator  phase  of 
development,  is  projected  to  be  constructed  using  50%  to  70%, 
by  weight,  composite  materials. 
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The  weight  savings  and  material  direction  property  tai¬ 
loring  advantages  must  be  weighted  against  the  possible  dis¬ 
advantages  of  composites.  A  key  disadvantage  of  polymer-based 
composites  is  that  they  exhibit  time  dependent  properties. 
These  time  dependent  or  viscoelastic  properties  must  be  fully 
characterized  in  order  to  predict  the  long-term  response  of 
the  material.  The  viscoelastic  characterization  is  complicated 
by  several  factors.  One  factor  is  the  time  required  to  perform 
arduous  long-term  characterization  tests.  Fortunately,  the 
test  time  has  been  successfully  reduced  by  applying  accelerated 
testing  techniques,  like  the  Time-Temperature  Superposition 
Principle  (TTSP)  [25,68].  Another  complicating  factor  is  the 
number  of  material  constants  required  to  describe  the  behavior 
of  orthotropic  composite  materials.  Metal  alloys,  normally 
considered  isotropic,  require  only  two  independent  material 
variables,  such  as  Young's  modulus  and  the  shear  modulus  to 
characterize  the  material.  Orthotropic  composite  materials 
require  nine  independent  material  variables  [29]  for  a  full 
three  dimensional  characterization.  Typically,  composites  are 
used  in  the  form  of  thin  sheets  with  the  loading  applied 
in-plane.  Under  these  conditions,  i.e.,  a  two  dimensional 
state  of  plane  stress,  the  number  of  independent  material 
constants  is  reduced  from  nine  to  four.  Additional  complication 
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arises  from  the  possible  dependence  of  the  four  material 
constants  on  temperature,  moisture,  ultraviolet  radiation, 
aging  time,  aging  temperature,  and  stress  level. 

This  study  concentrated  on  characterizing  the  effect  of 
thermoreversible  physical  aging  on  the  viscoelastic  creep 
response  of  a  composite  in  its  glassy  state.  Specifically, 
the  effects  of  aging  temperature  and  aging  time  on  the  creep 
response  of  Radel  X/IM7,  a  thermoplastic  composite,  were 
studied.  The  goal  was  to  predict  the  long-term  creep  response 
of  unidirectional  off-axis  composites  as  the  material  expe¬ 
riences  simultaneous  physical  aging  and  creep  deformation 
associated  with  an  applied  load. 

While  the  details  are  discussed  in  section  2.2,  the  concept 
of  thermoreversible  physical  aging  is  briefly  summarized  here. 
A  polymer,  quenched  from  the  rubbery  state  above  the  glass 
transition  temperature,  Tg,  to  some  lower  temperature,  does 
not  achieve  thermodynamic  equilibrium  [31] .  This  is  due  to 
the  rapid  increase  in  viscosity  as  Tg  is  approached  from  above. 
As  the  viscosity  increases,  the  motion  of  polymer  chains  are 
impaired  and  the  polymer  is  unable  to  reach  the  equilibrium 
degree  of  packing  before  solidification.  A  characteristic  of 
a  polymer,  in  the  quenched  state,  is  that  an  excess  of  trapped 
free  volume  exists  in  the  structure.  As  a  result,  the  polymer 
exhibits  volume  relaxation  with  time  as  the  structure  tends 
toward  the  equilibrium  state.  Struik  [62]  referred  to  this 
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slow  gradual  approach  to  volume  equilibrium  as  physical  aging. 
During  the  physical  aging  process,  many  of  the  physical  and 
mechanical  properties  of  the  material  change  [67],  These 
changes  include  an  increase  in  density,  tensile  yield  stress, 
and  elastic  modulus.  Also,  a  decrease  occurs  in  impact  strength, 
fracture  energy,  ultimate  elongation,  and  creep  rate.  This 
study  concentrates  specifically  on  investigating  how  the  creep 
response  of  a  composite  material  is  affected  by  physical  aging. 

1 . 1  Literature  Review 

The  literature  review  is  divided  into  two  sections.  First, 
a  review  of  the  creep  characterization  work,  completed  at 
Virginia  Polytechnic  Institute  and  State  University  (VPI  & 
SU) ,  is  presented.  This  section  outlines  the  results  of 
previous  efforts  to  characterize  the  creep  response  of  composite 
materials  and  illustrates  the  need  for  the  current  research 
effort.  Next,  a  summary  of  the  effects  physical  aging  has  on 
polymer-based  materials  is  presented.  Polymer  based  composites 
are  expected  to  show  similar  aging  effects,  since  the  polymer 
matrix  constituent  mainly  controls  the  viscoelastic  response 
of  the  composite. 

1.1.1  Creep  Characterization  at  VPI  6  SU 

Over  the  past  ten  years,  extensive  work  has  been  completed 
at  VPI  &  SU  to  improve  the  understanding  of  the  viscoelastic 
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response  of  polymer-based  composite  materials.  The  focus  of 
this  work  has  been  to  characterize  the  creep  compliance  and 
the  time  to  rupture  of  various  composite  material  systems. 

Yeow  et  al  [74]  initiated  the  viscoelastic  characterization 
work  at  VPI  &  SU  with  the  investigation  of  the  creep  response 
of  T300/934  graphite/epoxy  composite  material.  He  used  the 
Time-Temperature  Superposition  Principle  (TTSP)  to  shift  the 
results  of  short-term  (16  minute)  creep  tests,  at  various 
temperatures,  to  form  transverse,  ^221  and  shear,  S^,  master 
creep  compliance  curves.  Among  his  noteworthy  findings  was 
the  observation  that  the  fiber  direction  compliance,  5,,,  and 
fiber/transverse  coupling  compliance  term,  S,2,  were  both 
t ime- independent .  Also,  from  the  results  of  Yeow's  work, 
Morris  et  al  [47]  reported  the  coupling  compliance  terms  5,2 
and  52,  were  symmetric  with  only  a  small  difference  in  magnitude. 

Griffith  et  al  [26]  continued  the  viscoelastic  charac¬ 
terization  of  the  T300/934  composite  material.  He  generated 
master  compliance  curves  for  the  S22  and  directions  using 
both  temperature  and  stress  as  shifting  parameters.  The  fiber 
dominated  5„  compliance  term  was  found  to  have  a  slight  decrease 
with  temperature,  but  no  stress  dependence.  The  slight  tem¬ 
perature  dependence  was  attributed  to  the  viscoelastic  nature 
of  the  strain  gage  adhesive.  The  master  compliance  curves 
were  used  to  generate  compliance  predictions  for  30  and  60 
degree  off-axis  lamina  specimens.  The  predictions  had  fairly 
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good  agreement  with  the  test  data  up  to  a  test  time  of  about 
4  days.  After  this  period,  the  prediction  and  test  data  tended 
to  digress.  Griffith  also  expressed  concern  about  possible 
effect  of  physical  aging  on  the  T300/934  material.  In  order 
to  avoid  the  physical  aging,  he  heated  the  specimens  above  Tg 
and  then  cooled  them  at  a  very  slow  rate  (several  °F  per  hour) 
to  a  point  reasonably  below  T,.  Griffith  assumed  this  annealing 
process  would  freeze  the  equilibrium  amount  of  free  volume 
into  the  specimen  and  minimize  the  effects  of  physical  aging. 

Dillard  et  al  [13,14]  continued  working  with  the  T300/934 
graphite/epoxy  composite  material.  He  used  the  nonlinear 
Findley  power  law  to  represent  the  S22  and  creep  compliance 
terms  at  a  constant  temperature  of  160°C  (320°F)  .  The  power 
law  compliance  representation  was  then  incorporated  into  an 
incremental  numerical  computer  program,  based  on  lamination 
theory,  called  VISLAP.  The  program  was  used  to  predict  the 
creep  compliance  and  rupture  time  for  general  laminates  at 
160°C  (320°F)  .  Predictions  for  matrix  dominated  laminates, 
at  several  different  stress  levels,  agreed  with  experimental 
data.  Fiber  dominated  laminate  predictions  were  erroneously 
bounded  by  the  lamination  theory  assumption  that  no  interlaminar 
deformations  exist  [13] .  The  absence  of  interlaminar  defor¬ 
mations  fixed  the  predicted  laminate  deformation  at  the  value 
set  by  the  truss  network  formed  by  the  time-independent  fiber 
ply  directions  in  the  laminate. 
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Caplan  et  al  [7]  applied  nonlinear  viscoelastic  models 
developed  by  Findley  and  Schapery  to  polycarbonate  creep  and 
recovery  data  obtained  at  various  temperatures  and  stress 
levels.  The  creep  tests  lasted  30  minutes  and  were  followed 
by  a  60  minute  recovery  period.  The  Findley  power  law  model 
parameters  were  fit  to  the  data  using  a  linear  regression 
technique.  The  predictions  resulting  from  the  model  deviated 
from  the  experimental  data  at  longer  creep  times.  The  Schapery 
model  was  then  used  to  fit  the  data  using  a  graphical  technique. 
A  reasonable  fit  was  obtained  after  a  correction  for  unre¬ 
coverable  strain  effects.  No  attempt  was  made  during  this 
study  to  compare  the  predictive  capability  of  either  model 
beyond  the  30  minute  short-term  tests. 

Bertolotti  [4]  developed  a  numerical  procedure  for  eval¬ 
uating  the  parameters  in  Schapery' s  nonlinear  viscoelastic 
model.  The  procedure  was  based  on  the  Levenberg-Marquardt 
algorithm  [28] .  This  algorithm  finds  the  minimum  of  the  sum 
of  the  squared  error  between  the  data  points  and  a  nonlinear 
function.  This  numerical  fitting  procedure  provided  an  easy 
and  repeatable  means  of  fitting  the  Schapery  model  parameters 
to  actual  test  data. 

Hiel  et  al  [27]  continued  work  on  the  viscoelastic  char¬ 
acterization  of  T300/934  graphite/epoxy  composite  material 
following  Dillard.  He  used  the  numerical  procedure  developed 
by  Bertolotti  to  fit  the  Schapery  model  uo  experimental  data 
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generated  at  various  temperatures  and  stress  levels.  The 
stress  dependence  of  the  Schapery  model  parameters  was 
determined  for  the  shear  compliance  term.  It  was  also 
reported  that  the  transverse  compliance  S22  term  was  independent 
of  stress  during  the  short  time  creep  (1  hour)  test  conducted 
during  this  study.  However,  when  compared  to  long-term  test 
data  previously  generated  by  Griffith  et  al  [26],  it  was 
determined  that  the  stress  independence  of  the  S22  term  required 
further  examination.  Long-term  tests  to  verify  compliance 
predictions  derived  using  the  Schapery  model  were  not  conducted 
during  this  study. 

Yen  [73]  continued  the  creep  characterization  work  at  VPI 
&  SU  but  switched  from  orientated  laminated  composite  materials 
to  a  random  fiber  composite  material  known  as  SMC-R50.  He 
developed  a  numerical  procedure,  similar  to  Bertolotti' s,  to 
fit  the  Findley  model  parameters  to  creep  data.  It  was  found 
that  the  time  exponent,  "n"  in  the  Findley  equation,  depended 
on  the  duration  of  the  creep  test,  but  approached  an  asymptotic 
value  for  long  tests.  The  predictions  derived  form  the  Findley 
equation  accurately  predicted  long-term  creep  response.  In 
addition,  the  creep  response  due  to  multiple  step  loadings  was 
predicted  using  a  modified  superposition  principle  with  suc¬ 
cess.  Yen  noted  that  the  random  fiber  composite  required 
mechanical  conditioning  prior  to  each  creep  and  recovery  test 
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in  order  to  obtain  repeatable  results.  Additionally,  Yen  found 
the  results  from  the  conditioned  specimen  tests  were  within 
the  data  scatter  band  of  unconditioned  specimens. 

Tuttle  et  al  [68]  investigated  the  viscoelastic  creep 
response  of  T300/5208  graphite  epoxy  composite  laminates.  He 
used  a  480/120  minute  creep/recovery  test  cycle  on  0,  10,  and 
90  degree  specimens,  at  various  stress  levels,  with  the  tem¬ 
perature  held  constant  at  149°C  (353°F)  .  The  Schapery  nonlinear 
viscoelastic  model  parameters  were  successfully  fit  to  the 
test  data  using  the  numerical  technique  developed  by  Bertolotti  . 
Tuttle  noted  that  long-term  predictions  from  the  Schapery  model 
were  most  sensitive  to  errors  in  the  exponent  "n".  Long-term 
creep  tests  of  unidirectional  laminates,  to  verify  the  pre¬ 
dictive  capability  of  the  Schapery  model,  were  not  conducted 
during  this  study.  Instead,  long-term  tests,  lasting  10^ 
minutes,  were  conducted  on  two  multiangle  symmetric  laminates. 
Analytical  creep  compliance  predictions  were  obtained  using 
the  VISLAP  program  written  by  Dillard.  The  program  was 
modified  to  utilize  the  Schapery  nonlinear  viscoelastic  model 
instead  of  the  Findley  power  law.  The  long-term  testing 
revealed  that  the  predicted  viscoelastic  response  was  sig¬ 
nificantly  less  than  the  measured  response.  This  discrepancy 
was  not  definitely  explained,  but  was  thought  to  be  due  to 
stress  interaction  effects  within  individual  plies  or  due  to 
interlaminar  shear  deformations  effects. 
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Owens  [51]  conducted  an  experimental  investigation  of  the 
material  properties  of  graphite/polyphenylene  sulfide,  a 
thermoplastic  composite  system.  His  creep  tests,  conducted 
on  (±462)  s  specimens  at  177°C  (350°F) ,  produced  results  of 
particular  interest  to  the  present  study.  The  testing  consisted 
of  two  72  hour  creep  tests  run  sequentially.  The  first  creep 
test  was  followed  by  a  recovery  period  until  the  strain  returned 
to  zero.  The  second  test  was  then  started  with  the  temperature 
continuing  to  be  held  at  177°C  (350°F)  .  Owens  observed  that 
the  stiffness  significantly  increased  between  the  two  creep 
tests.  Unable  to  explain  this  phenomenon,  the  test  sequence 
was  repeated  with  a  new  specimen.  Again,  the  stiffness  was 
found  to  increase.  It  is  well  known  (see  Struik  [59])  that 
polymer  based  materials  increase  in  stiffness  as  they  physically 
age.  Thus,  physical  aging,  occurring  simultaneously  with  the 
creep  deformations  during  the  first  test,  offers  an  explanation 
for  the  increase  in  stiffness  during  the  second  creep  test. 

Gramoll  et  al  [25]  studied  the  thermoviscoelastic  response 
of  Kevlar  49/Fiberlite  7714A  epoxy  composite  lamina  and  lam¬ 
inates.  Both  the  matrix  and  fiber  in  this  material  system 
exhibit  viscoelastic  behavior.  Thus,  all  four  material 
compliance  terms  were  time  dependent.  The  viscoelastic  response 
was  determined  using  short-term  (20  minute)  tests  at  various 
temperatures  and  stress  levels.  Gramoll  used  TTSP  and  the 
Findley  power  law  to  model  the  linear  viscoelastic  behavior. 
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Nonlinear  effects  were  modeled  with  a  new  quadratic  power  law 
developed  during  the  study  to  facilitate  numerical  implemen¬ 
tation.  A  new  program  called  VCAP  (Viscoelastic  Composite 
Analysis  Program)  was  written  to  predict  the  viscoelastic 
response  of  composite  laminates.  The  program  was  based  on  a 
nonlinear  differential  equation  formulation  using  a  series  of 
Kelvin  elements  (i.e.  Prony  series)  to  represent  the  viscoe¬ 
lastic  material  behavior.  Numerical  predictions  for  unidi¬ 
rectional,  two  direction  and  three  direction  laminates  agreed 
very  well  with  the  4  week  long  test  data  generated  during  the 
study. 

1.1.2  Physical  Aging  Effects 

Interest  in  the  effects  of  physical  aging  on  polymeric 
material  properties  dates  back  to  early  work  by  Kovacs  [31]. 
His  research  showed  that  the  volume  of  a  polyvinylacetate 
sample,  after  a  quench  from  far  above  Tg  to  a  constant  tem¬ 
perature  below  Tg,  slowly  decreased  with  time.  Struik  [59] 
studied  more  than  35  amorphous  polymers  and  found  that  the 
small-strain  mechanical  properties  changed  dramatically  during 
this  volume  relaxation  process.  Struik  reported  the  effect 
of  physical  aging  on  both  the  elastic  and  viscoelastic  material 
properties.  Among  his  findings  were  the  observations  that 
during  physical  aging  the  material  increases  in  stiffness  and 
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becomes  more  brittle.  He  was  able  to  characterize  the  rate 
at  which  physical  aging  occurs  and  developed  a  model  to  predict 
the  long-term  viscoelastic  response  of  amorphous  polymers. 

Fundamental  to  Struik's  model  was  the  observation  that 
momentary  creep  compliance  curves  exhibit  an  universal  shape 
in  the  temperature  range  associated  with  physical  aging. 
Momentary  compliance  curves  are  formed  using  results  of  creep 
tests  where  the  test  duration  remains  small  compared  to  the 
time  measured  starting  after  the  quench  from  above  T,  and 
continued  till  the  start  of  the  creep  test,  referred  to  as  the 
aging  time.  Momentary  tests  are  essentially  a  "snapshot"  of 
the  compliance  properties  taken  at  a  constant  aging  time.  The 
existence  of  this  universal  shape  of  the  momentary  creep 
compliance  curves  has  been  reported  by  other  investigators 
[63] .  The  universal  shape  implies  that  the  relaxation  spectrum 
does  not  change  with  temperature  or  aging  time  [6] .  Thus,  the 
relaxation  times  contained  in  the  relaxation  spectrum  are  all 
identically  shifted  by  temperature  or  aging  time  [70]  .  The 
effect  of  temperature  or  aging  time  on  short-term  momentary 
creep  compliance  curves  was  seen  as  a  horizontal  shift  of  the 
curve  on  the  log-time  scale.  The  marked  curvature  of  the 
compliance  curves  allows  the  shift  factor  for  each  short-term 
test  to  be  unambiguously  determined.  A  master  momentary  creep 
compliance  curve  was  formed  by  shifting  several  short-term 
curves  obtained  at  different  temperatures  or  aging  times.  The 
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actual  long-term  creep  response  was  obtained  using  an  effective 
time  theory  developed  by  Struik  with  the  derivation  details 
found  in  section  2. 2. 5. 4.  The  theory  is  based  on  the  assumption 
that  the  relaxation  times  all  have  the  same  aging  dependence. 
The  aging  dependence  is  characterized  by  the  aging  rate  of  the 
material.  The  theory  is  applied  by  replacing  the  current  time 
t,  used  in  the  momentary  master  creep  curve,  with  an  effective 
time  X.  The  effective  time  accounts  for  the  shifting  of  the 
relaxation  times  as  the  material  ages. 

Booij  et  al  [6]  noted  that  the  physical  aging  effects  are 
influenced  by  the  rate  the  material  is  cooled  from  above  Tg. 
This  influence  is  due  to  the  free  volume  dependence  on  the 
cooling  rate.  At  faster  cooling  rates,  the  difference  between 
the  material  volume  state  and  the  equilibrium  volume  state  is 
larger  than  it  is  at  slower  cooling  rates.  The  larger  volume 
difference  equates  to  a  greater  change  in  the  physical  prop¬ 
erties  of  the  material  during  the  physical  aging  process. 

During  physical  aging,  the  enthalpy  of  the  material  exhibits 
the  same  relaxation  phenomena  as  volume  [38] .  The  enthalpy 
and  heat  capacity,  i.e.,  first  derivative  of  enthalpy  with 
respect  to  temperature,  at  a  constant  pressure,  are  thermo¬ 
dynamic  quantities  that  can  be  monitored  using  differential 
scanning  calorimetry  (DSC)  equipment.  The  effect  of  physical 
aging  on  these  thermodynamic  properties  has  been  reported  by 
several  investigators  [10,49,67].  They  noted  a  correlation 
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between  aging  time  and  the  location  of  the  endothermic  peak 
as  the  temperature  was  raised  during  a  DSC  test.  This  cor¬ 
relation  was  used  as  a  method  to  monitor  the  extent  of  physical 
aging  in  the  material. 

The  effects  of  physical  aging  continue  until  the  material 
reaches  volume  equilibrium.  The  time  required  to  reach  volume 
equilibrium  depends  on  the  aging  temperature.  At  aging  tem¬ 
peratures  less  than  T,  .ninus  27°F  (15°C)  ,  the  material  will  not 
reach  thermodynamic  equilibrium  in  a  practical  time  frame  [66]  . 
Therefore,  polymer  products  used  in  the  aging  temperature  range 
can  be  expected  to  change  physical  properties  during  a  normal 
service  life  of  10  -  50  years  [59] .  Predicting  the  physical 
aging  effect  on  the  material  properties  has  prompted  consid¬ 
erable  research.  Several  recent  papers  have  reported  physical 
aging  effects  on  neat  resin  and  composite  laminate  materials. 

Vleeshouwers  et  al  [70]  recently  reported  the  effects  of 
physical  aging  on  the  creep  and  stress  relaxation  properties 
of  cured  EPON  828/Epoxy  in  the  linear  viscoelastic  range.  They 
successfully  measured  the  aging  rate  at  several  temperatures 
below  Tg.  The  short-term  momentary  creep  and  stress  relaxation 
curves  were  successfully  shifted  to  form  master  curves. 
Long-term  test  results  were  not  reported  to  verify  the  accuracy 
of  the  master  curves. 
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Work  by  Wang  and  Ogale  [71]  showed  the  effect  of  physical 
aging  on  the  relaxation  and  dynamic  mechanical  moduli  properties 
of  neat  epoxy  and  an  epoxy/carbon  fiber  reinforced  composite 
at  a  constant  test  temperature.  It  was  reported  that  both  the 
epoxy  and  composite  had  the  same  aging  shift  rate.  Their 
conclusion  was  that  the  fiber  reinforcing  did  not  alter  the 
aging  characteristics  of  the  epoxy. 

The  most  recent  physical  aging  study  on  composites  was 
performed  by  Sullivan  [66].  He  worked  with  neat  Dow  Derakane 
470-36  resin  and  an  unidirectional  glass  reinforced  composite. 
The  thermoset  resin  and  composite  were  both  found  to  have  the 
same  aging  shift  rate.  This  observation  is  consistent  with 
the  work  of  Wang  and  Ogale  [71] .  Momentary  creep  master  curves 
were  constructed  for  the  neat  resin  and  the  four  composite 
compliance  directions.  The  and  5,2  directions  were  found 
to  be  independent  of  physical  aging.  Long-term  predictions, 
using  Struik's  effective  time  theory,  agreed  with  experimental 
data  up  to  about  200  hours  of  creep.  After  this  time  period, 
the  predictions  followed  the  same  data  trend  but  differed  in 
the  magnitude  of  the  creep  compliance. 

1 . 2  Overview  of  Proposed  Study 

The  goal  of  this  study  was  to  determine  the  effect  of 
physical  aging  on  a  thermoplastic  composite  material.  Radel 
X/IM7  composite  material  was  chosen  for  evaluation.  This 
material  is  currently  being  considered  by  NASA  Langley  for  use 
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in  the  next  generation  high  speed  transport  aircraft.  The 
operating  environment  for  this  aircraft  will  be  near  188°C 
(370°F)  with  the  service  life  of  60,000  to  120,000  hours  at 
temperature.  Accurate  predictions  of  the  viscoelastic  prop¬ 
erties  of  the  material  are  essential  to  insure  the  material 
meets  the  design  strength  and  stiffness  requirements  for  the 
entire  service  life. 

Three  project  objectives  were  defined  to  meet  the  goal  of 
this  study.  The  first  objective  was  to  develop  the  experimental 
procedures  required  to  conduct  physical  aging  creep  experi¬ 
ments.  Several  experimental  problems  were  encountered  during 
preliminary  creep  testing.  In  particular  problems  with  strain 
measurement,  data  acquisition,  oven  thermal  stability,  and 
specimen  quenching  procedures  were  experienced.  The  exper¬ 
imental  procedures  developed  to  solve  these  problems  were 
documented  to  ensure  future  aging  experiments  are  conducted 
in  a  repeatable  and  consistent  manner. 

The  second  objective  was  to  apply  the  developed  experimental 
procedures  to  characterize  the  influence  of  sub  Tg  physical 
aging  on  the  creep  response  of  Radel  X/IM7  composite  material. 
Short-term  creep  tests  were  used  to  determine  the  aging  rate 
and  its  dependence  on  aging  temperature.  In  addition,  various 
ply  angle  orientations  were  tested  to  determine  the  effect  of 
the  physical  aging  on  each  of  the  four  compliance  matrix  terms 
(S|i,  5)2,  522,  and  S^)  in  the  linear  viscoelastic  range. 
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The  third  and  final  objective  was  to  conduct  long-term 
creep  tests  to  verify  the  influence  of  physical  aging  on  the 
creep  response  of  the  Radel  X/IM7  composite  material.  The 
long-term  tests  were  initially  conducted  on  each  of  the  ply 
angle  orientations  tested  during  determination  of  the  com¬ 
pliance  matrix  terms.  The  test  results  were  then  compared 
with  theoretical  predictions  calculated  using  Struik's 
effective  time  theory.  Finally,  a  30  degree  off-axis  specimen 
was  tested  and  compared  to  the  creep  compliance  prediction 
obtained  by  transforming  the  compliance  matrix  to  the  specimen' s 
angle  orientation. 
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2  BACKGROUND  INFORMATION 


Viscoelasticity,  as  its  name  implies,  is  a  combination  of 
the  studies  of  viscous  fluids  and  elasticity.  Viscoelastic 
materials  are  characterized  by  having  time-dependent  mechanical 
behavior  with  a  memory  effect.  The  discussion  herein  will 
describe  the  terminology  and  analysis  models  used  to  describe 
viscoelastic  material  behavior.  Also,  the  effects  of  physical 
aging  on  the  viscoelastic  response  of  the  material  are  dis¬ 
cussed.  Additional  information  on  the  subject  is  located  in 
several  reference  books  [1,11,16,22,23,59]. 

2 . 1  Viscoelastic  Materials 

Viscoelastic  materials  exhibit  time-dependent  mechanical 
behavior  with  memory  effects.  This  type  of  material  response 
is  demonstrated  when  a  viscoelastic  tensile  specimen  is  loaded 
with  a  constant  force.  A  viscoelastic  material  will  initially 
elongate  when  the  load  is  applied,  as  in  any  normal  elastic 
material,  but  will  also  continue  to  elongate  with  time  while 
the  load  remains  applied.  This  continuing  time-dependent 
elongation  is  referred  to  as  creep.  During  this  creep  process, 
the  apparent  Young's  modulus  E  will  steadily  decrease  as  the 
specimen  elongates  with  time.  Upon  removal  of  the  load,  the 
viscoelastic  material  will  remember  its  original  configuration 
and  will  tend  to  recover  to  that  configuration  with  time. 
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2.1.1  Common  Tests 


Three  types  of  tests  are  commonly  used  to  characterize  the 
mechanical  response  of  a  viscoelastic  material.  The  tests 
used  are  the  creep/creep  recovery  test,  the  stress  relaxation 
test,  and  the  constant  strain  rate  test.  The  least  complex 
test  to  conduct  of  the  three  is  the  two  part  creep/creep 
recovery  test.  The  first  part,  known  as  the  creep  portion, 
is  defined  by  the  response  up  to  time  ti.  During  this  time 
segment,  a  constant  uniaxial  step  input  load  is  applied  to  the 
specimen,  as  shown  in  figure  2.1.  The  viscoelastic  material 
responds  with  an  initial  elastic  strain  followed  by  a 
continually  increasing  creep  strain  E^,  also  shown  in  figure 
2.1.  The  second  part,  known  as  the  recovery  portion,  starts 
after  time  ti,  when  the  input  load  is  removed.  During  this 
portion,  a  reverse  elastic  strain  is  followed  by  recovery 
of  a  portion  of  the  creep  strain  at  a  continuously  decreasing 
rate,  as  shown  in  the  material  response  illustration  in  figure 
2.1.  The  viscoelastic  material  may  recover  to  the  entire  creep 
strain  if  sufficient  time  is  allowed  for  recovery.  This  is 
the  simplest  test  of  the  three  tests  to  conduct,  because  the 
applied  force  is  either  a  constant  or  zero.  Additionally,  the 
creep  compliance  D(t)  is  easily  calculated  by  dividing  the 
strain,  measured  at  the  center  of  the  specimen,  by  the  applied 
constant  stress  value.  Specimens  for  this  test  are  designed 
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with  adequate  length  to  insure  the  effect  of  the  end  constraint, 
provided  by  the  gripping  fixtures,  has  a  minimal  effect  on  the 
stress  at  the  center  of  the  specimen. 

The  next  increasingly  complex  test  which  can  be  conducted 
is  the  stress  relaxation  test.  This  test  measures  the  stress 
relaxation  in  the  specimen  as  it  is  subjected  to  a  constant 
uniaxial  strain  input  load,  shown  as  in  figure  2.2.  The 
input  load  produces  an  initial  stress  0^,  which  slowly  relaxes 
with  increased  time,  as  shown  in  the  material  response 
illustration  in  figure  2.2.  The  increased  complexity  of  this 
test  is  partly  due  to  the  equipment  required  to  hold  a  constant 
strain  on  the  test  specimen.  In  order  to  accomplish  this,  the 
testing  machine  must  be  able  to  decrease  the  load  applied  to 
the  specimen  as  time  increases.  Another  complication  is  the 
requirement  to  minimize  the  effect  of  end  constraint  provided 
by  the  grips  on  a  anisotropic  specimen.  The  assumption  of 
constant  strain  for  the  length  of  the  specimen  cannot  be 
satisfied  at  the  specimen  attachment  points  unless  proper  grip 
rotation  is  provided. 

The  final  test  discussed  herein  is  the  constant  strain 
rate  test.  This  test  was  suggested  by  Smith  [58]  to  determine 
the  relaxation  modulus  E(t)  of  elastomers.  During  this  test, 
the  viscoelastic  material  is  loaded  with  an  uniaxial  strain 
that  is  increased  as  some  constant  rate  R,  as  shown  in  figure 
2.3.  The  axial  load  required  to  produce  the  constant  strain 
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rate  is  monitored  during  the  test  and  used  to  calculate  the 
stress  response.  Viscoelastic  materials  respond  to  the  loading 
with  a  nonlinear  increase  in  stress  with  increasing  time,  as 
shown  in  the  material  response  illustration  of  figure  2.3. 
The  relaxation  modulus  is  found  by  calculating  the  instantaneous 
slope  of  the  stress  strain  curve.  This  test  is  also  more 
difficult  to  conduct  in  view  of  the  increased  complexity  of 
the  test  machine  required  to  produce  a  constant  strain  rate 
in  the  specimen. 

2.1.2  Linear  Viscoelasticity 

2 . 1 . 2 . 1  Definition 

Viscoelastic  materials,  like  elastic  materials,  may  be 
classified  as  either  linear  or  nonlinear.  Findley  et  al  [22] 
and  Schapery  [57]  both  define  a  viscoelastic  material  as  linear 
when  stress  is  proportional  to  strain  at  a  given  time  and  when 
linear  superposition  holds.  Mathematically  these  requirements 
are  represented  by  the  equations: 

e[co(01  =  ce[a(t)]  (2.1) 

e[a,(t)  +  =  E[a,(t)] + eLOjCt-t,)]  (2.2) 
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where  <5  is  the  input  stress,  e  is  the  output  strain  and  c  is 
a  constant.  The  first  requirement,  equation  2.1,  states  that 
the  strain  response  of  a  scalar  times  the  input  stress  equals 
the  scalar  times  the  strain  response  of  the  original  input 
stress.  This  is  illustrated  in  figure  2.4  (a).  The  second 
requirement,  equation  2.2,  states  that  the  combined  strain 
response  of  two  different  loadings  applied  at  different  times 
is  the  sum  of  the  strain  outputs  of  each  load  acting  separately 
and  is  illustrated  in  figure  2.4  (b)  .  This  second  requirement 
is  also  commonly  referred  to  as  the  Boltzmann's  Superposition 
Principle . 

2. 1.2. 2  Boltzmann's  Superposition  Principle 

The  concept  of  summing  the  effects  of  various  input  loads 
can  be  extended  to  arbitrary  stress  or  strain  histories.  This 
is  accomplished  by  approximating  the  arbitrary  stress  history 
with  a  series  of  steps  of  stress  inputs  as  shown  in  figure 
2.5.  Mathematically  this  is  represented  by: 

a(t)  =  AaoH(t)  +  Aa,H(t  -  t,)  +  AajHCt  -  Xj)  +  •  •  •  +  AajHCt  -  Xj)  (2.3) 


where  the  Heaviside  function  H(t)  is  defined  as: 


H(t)  = 


1 0,  when  t  <  0 
[  l,when  t^O 


(2.4) 
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with  the  strain  response  to  this  stress  input  given  by: 


e(t)  =  D(t)Aao  +  D(t  -  'C,)Aa,  +  D(t  -  +  •  •  •  +  D(t  -  ti)Aai  (2.5) 


where  D(t-Xi)  is  the  compliance  for  each  time  step.  If  the  time 

intervals  between  steps  shrinks  to  an  infinitesimal  amount, 
an  integral  expression  of  the  total  strain  is  given  by: 


e(t)=  r‘D(t-x)H(t-t)d[o(T)]. 

Jo 


(2.6) 


Assuming  the  stress  history  is  differentiable  and  knowing  that 
the  Heaviside  function  is  unity  over  the  range  of  integration, 
equation  2.6  reduces  to  the  following  form: 


e(t)  =  J^  D(t-T)^^dx. 


(2.7) 


Similarly,  an  expression  for  the  stress  response  to  an  arbitrary 
strain  history  is  given  by: 


a(t)  =  J^  E(t-x)^^dx 

where  E(t)  is  the  relaxation  modulus. 


(2.8) 
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2 . 1 . 2 . 3  Mechanical  Models 


Many  types  of  mechanical  analog  models  can  be  used  to  model 
viscoelastic  behavior.  Mechanical  models  are  constructed  using 
various  combinations  of  linear  spring  and  linear  dashpot 
elements.  The  stress  and  strain  in  a  linear  spring  are  related 
by  the  equation: 


<y.pnn«  =  <5.  =  Ee 


(2.9) 


where  E  is  interpreted  as  a  linear  spring  constant  or  Young's 
modulus.  The  stress  in  a  linear  dashpot  is  related  to  the 
strain  rate  in  the  dashpot  by: 

d£ 

=  =  =  (2.10) 

where  is  called  the  coefficient  of  viscosity.  The  two 
simplest  mechanical  models  are  derived  by  combining  a  single 
spring  and  single  dashpot  element.  The  Maxwell  model  connects 
the  two  elements  in  series,  whereas  the  Kelvin  model  uses  a 
parallel  arrangement.  Both  models  are  graphically  represented 
in  figure  2.6. 

The  Maxwell  model  is  examined  in  more  detail  first.  In 
this  model,  the  total  strain  and  total  strain  rate  are  calculated 
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by  summing  the  component  carried  by  each  of  the  two  elements. 
These  relationships  and  the  stress  definitions  in  equations 
2.9  and  2.10,  are  used  to  derive  the  following  stress-strain 
equation: 


a  a 
£=-+-. 
E  \i 


(2.11) 


This  differential  equation  can  be  solved  for  various  stress 
or  strain  inputs.  Classical  solutions  include  the  strain 
response  to  a  constant  stress  input  a  =  ao,  applied  at  t  =  0, 
given  by  the  following  relationship; 


(2.12) 


or  the  stress  response  to  a  constant  strain  input  e=0o/  applied 
at  t  =  0,  given  by  the  relationship; 

o(t)  =  Eeoe"‘^‘''‘>  =  EEoe’^‘'^>  (2.13) 

where  X  =  \L/E  is  called  the  relaxation  time. 

The  Kelvin  model  uses  a  parallel  arrangement  of  the  spring 
and  dashpot  elements.  In  this  model,  the  total  stress  carried 
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by  the  model  is  calculated  by  summing  the  amount  of  stress 
carried  by  each  of  the  two  elements.  This  relationship  is 
used  with  equations  2.9  and  2.10  to  derive  the  following 
stress-strain  equation: 


.  E  a 

e+-e=-.  (2.14) 

This  differential  can  be  solved  for  creep  loading  with  C  =  ao, 
applied  at  t=0.  The  resulting  strain  solution  is: 


e(t)=|:(l-e-<“"’)=|(l-e-"'‘>)  (2.15) 

where  X=n/E  is  called  the  retardation  time. 

The  use  of  the  Maxwell  and  Kelvin  models  to  accurately 
represent  viscoelastic  materials  is  limited.  This  is  due  to 
the  inability  of  the  models  to  describe  the  entire  actual 
material  response.  For  example,  the  Maxwell  model  does  not 
show  time-dependent  recovery  after  a  constant  stress  input 
load  is  removed.  Additionally,  upon  loading  or  unloading,  the 
Kelvin  model  fails  to  exhibit  time-independent  strain  that 
usually  occurs  with  actual  materials.  More  complex  mechanical 
models  have  been  developed  by  combining  Maxwell  and  Kelvin 
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models  in  various  combinations.  These  variations  will  not  be 
described  here,  but  may  be  found  in  a  book  by  Findley  et  al 
[22]  . 

One  complex  model  that  is  widely  used  for  viscoelastic 
modeling  where  the  stress  history  has  been  prescribed  is  the 
generalized  Kelvin  model,  illustrated  in  figure  2.7.  This 
model,  sometimes  referred  to  as  the  Prony  series,  consists  of 
n  Kelvin  elements  connected  in  series.  The  Prony  series  offers 
the  advantage  that  each  Kelvin  element  can  be  solved  inde¬ 
pendently.  The  total  solution  is  obtained  by  summing  each  of 
the  individual  Kelvin  element  solutions.  The  Prony  series 
solution  for  creep  loading  with  a  =  ao,  applied  at  t  =  0,  is  written 
as : 


e(t)  =  o,  i  l(i  =  a.  i  l(i  -  c-"'*'’). 

i  =  l  Ej '  '  i  =  l  Ej '  ' 


(2.16) 


An  initial  elastic  time-independent  strain  can  be  added  to  the 
model  by  placing  a  single  spring  element  in  series  with  the 
Kelvin  elements.  The  solution  for  creep  loading  becomes: 


(2.17) 
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where  E,  is  the  spring  constant  of  the  additional  single  spring 


element . 

The  obvious  disadvantage  to  the  Prony  model  is  the  large 
number  of  material  constants  required  to  describe  the  material . 
The  model  requires  both  E  and  |i  for  each  Kelvin  element  used 
in  the  summation.  Gramoll  et  al  [25]  used  the  Prony  series 
in  his  work  and  suggests  spacing  the  retardation  times  A^  =  |Xi/Ei 
at  about  one  per  each  decade  of  time  that  is  modeled.  Gramoll 
solved  for  the  numerous  material  constants  in  the  Prony  series 
by  initially  fitting  the  viscoelastic  data  with  a  power  law 
model.  Then,  he  utilized  a  numerical  procedure  to  fit  the 
Prony  series  constants  to  the  power  law  model  for  each  time 
decade. 

2 . 1 . 2 . 4  Power  Law  Model 

Findley  et  al  [17,19,20,22,33,34]  observed  that  the  fol¬ 
lowing  povjer  function  of  time  empirically  fit  creep  data  of 
many  linear  rigid  plastics  over  a  wide  time  span: 

e(t)  =  eo+mt"  (2.18) 

where  Eq, m, n  are  material  constants  of  the  model.  Findley  et 

al  [22]  have  reported  that  the  material  constants  £o  and  m  are 
functions  of  stress  and  temperature,  whereas  n  is  independent 
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of  stress  and  nearly  independent  of  temperature. 

Several  investigators  at  VPI  &  SU  have  used  the  Findley 
power  law  model  to  represent  viscoelastic  behavior.  Notably, 
Dillard  el  al  [13]  found  that  the  coefficients  of  the  power 
law  equation  may  become  singular  when  fit  to  creep  data  using 
a  fitting  procedure  developed  by  Boiler  [5]  .  Additionally, 
Dillard  noted  how  small  errors  in  experimental  data  may  result 
in  large  errors  in  the  determined  power  law  material  constants. 
Yen  [73]  found  the  material  constants  using  a  Levenberg- 
Marquardt  algorithm,  available  on  IMSL  [28]  FORTRAN  subrou¬ 
tines,  to  minimize  the  sum  of  the  squares  of  error  between  the 
data  points  and  the  power  law  function.  This  method  provided 
an  accurate  and  consistent  method  for  determining  the  power 
law  constants. 

2 . 1 . 2 . 5  Struik  Model 

Struik  [59],  in  his  studies  of  physical  aging  of  polymers, 
experimentally  determined  that  momentary  creep  curves,  in  the 
small  strain  behavior  range  of  polymeric  materials,  have  the 
same  shape  when  plotted  on  double-logarithmic  scale.  The 
requirement  for  momentary  creep  testing,  i.e.,  when  the  creep 
test  duration  is  short  compared  to  the  previous  aging  time, 
is  fully  explained  in  section  2.2.  Struik  found  the  shape  of 
the  creep  compliance  curve  could  be  described  by  the  equation: 
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(2.19) 


D(t)  =  Doe 


where  Dq  represents  the  compliance  immediately  after  stress 

application,  to  is  variable  proportional  to  the  mechanical 
relaxation  time  and  m  is  a  material  constant.  In  uniaxial 
tensile  creep  tests,  the  applied  stress  a  is  constant  and  is 
related  to  the  compliance  and  strain  by  the  equation: 


(2.20) 


Using  equation  2.20,  equation  2.19  is  rewritten  in  terms  of 
strain  response  as: 


E(t)  =  eoe 


(2.21) 


where  eo  =  oDo.  The  analytical  representation  of  the  compliance 

response  (equation  2.19)  is  useful  for  determining  logarithmic 
shift  factors  used  in  time-aging  shifting  and  time-temperature 
superposition . 

2 . 1 . 2 . 6  Time-Temperature  Superposition  Principle 

The  goal  of  accelerated  characterization  of  a  material  is 
to  develop  a  method  that  uses  short-term  testing  results  for 
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the  prediction  of  long-term  behavior.  Time-Temperature 
Superposition  Principle  (TTSP) ,  also  known  as  The  Method  of 
Reduced  Variables,  uses  temperature  as  a  parameter  to  accelerate 
time  dependent  processes,  like  creep.  Leaderman  [35]  introduced 
TTSP  in  1943  as  an  empirical  observation  about  the  shape  of 
creep  compliance  curves  at  elevated  temperatures.  He  noted 
that  compliance  curves  at  various  temperatures  are  basically 
the  same  shape,  but  displaced  along  the  time  axis.  The  principle 
is  illustrated  by  shifting  several  compliance  curves,  each 
measured  at  a  different  temperature,  along  the  time  axis  to 
form  the  master  compliance  curve  shown  in  figure  2.8.  The 
amount  of  horizontal  shifting,  labeled  Et,  and  Ejj  respectively 
for  each  curve  in  figure  2.8,  is  a  function  of  the  temperature 
differential  between  the  shifted  curve  and  the  curve  chosen 
as  the  reference. 

Two  functional  forms  for  the  shift  factor  function  are 
widely  used  and  work  for  many  materials.  First,  when  the 
temperature  is  less  than  the  glass  transition  temperature 
(T<Tg),  the  Arrhenius  equation  [16]  is  usually  applicable  and 
is  written  as: 


1  _  ah  f  1  1  ' 

2.303  Tr^ 


(2.22) 
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where  AH  is  the  constant  activation  energy  (per  mole) ,  R  = 
1.987  (cal/mole  K)  is  the  universal  gas  constant,  and  both  the 
temperature  T  and  reference  temperature  Tr  are  in  degrees 
Kelvin. 

The  second  equation,  called  the  Williams-Landel-Ferry  (WLF) 
equation  [72],  normally  applies  when  the  temperature  is  greater 
than  the  glass  transition  temperature  (T>Tg)  .  The  WLF  equation 
is  given  as: 


,  -C,(T-Tr) 

(Q  +  T-Tr) 


(2.23) 


where  C,  and  Q  are  constants,  and  the  reference  temperature 

Tr  is  generally  equal  to  the  glass  transition  temperature  or 
higher . 

The  shift  factor  ax  relates  the  master  compliance  response 

at  the  reference  temperature  to  the  response  at  any  other 
temperature  using  the  "reduced  time"  approach.  The  compliance 
at  any  temperature  T  follows  the  relationship: 

D(T,t)  =  D(TR,C)  (2.24) 


where  t  is  the  actual  time  measured  from  the  start  of  load 
application,  Tr  is  the  reference  temperature  and  ^  is  the 
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"reduced  time".  During  transient  temperature  conditions  i.e., 
when  aT{T(x)}  is  a  function  of  temperature  and  time,  the  reduced 
time  is  related  to  the  actual  time  by: 

aT{T(x)} 

where  x  is  an  arbitrary  time  variable  prior  to  t.  This 
relationship  is  simplified  for  constant  temperature  testing 
to : 


(2.26) 


where  aj  does  not  change  with  time.  The  relationship  between 

the  reference  creep  curve  and  a  different  temperature  creep 
curve  becomes  more  obvious  when  the  log  of  equation  2.26  is 
taken. 


log  C  =  log  t- log  Or  (2.27) 

It  is  easily  noted  that  the  new  temperature  creep  curve  is  a 
shift  of  the  reference  curve  in  log  time  by  an  amount  of  loga-f. 


34 


The  class  of  material  where  the  viscoelastic  response  at 
various  temperatures  can  be  manipulated  with  translation  or 
horizontal  shifting  to  form  smooth  master  curve  is  referred 
to  as  a  "Thermorheologically  Simple  Material"  (TSM) .  Often, 
simple  horizontal  shifts  fail  to  produce  a  smooth  master  curve. 
In  these  cases,  additional  vertical  shifts  of  the  data  may  be 
required  to  obtain  a  smooth  master  curve.  When  this  additional 
shifting  is  required,  due  to  moisture,  stress  or  other 
influences,  the  material  is  referred  to  as  a  "Thermorheolog¬ 
ically  Complex  Material"  (TCM) .  Griffith  et  al  [26]  reviews 
the  shifting  procedures  for  TCM.  He  reported  that  the  graphical 
normalization  procedure  provided  the  best  average  material 
properties . 

2.1.3  Nonlinear  Viscoelasticity 

Materials  that  fail  to  meet  the  linearity  requirements 
stated  by  equations  2.1  and  2.2  are  called  nonlinear.  The 
demarcation  between  linear  and  nonlinear  is  somewhat  arbitrary 
depending  on  the  accuracy  required  in  the  analysis.  Testing 
variables,  i.e.,  stress  level,  temperature,  moisture  and  aging 
time,  can  each  influence  the  nonlinear  viscoelastic  response 
of  a  material.  The  Findley  and  Schapery  models  have  gained 
popularity  for  modeling  nonlinear  effects  and  are  discussed 
in  further  detail. 
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2. 1.3.1  Findley  Model 


Findley  et  al  [18,21]  proposed  modifying  the  simple  power 
law  function  for  strain  response  given  as: 

e(t)  =  £^+m^  (2.18) 


where  Eo, m  and  n  are  material  constants,  to  include  nonlinear 

stress  effects  by  making  Eq  and  m  functions  of  stress.  The 
function  chosen  was  in  the  form  of  a  hyperbolic  sine  such  that: 


Eo  =  £'oSinh 


m  =  in'sinh 


a 


n  =  constant  ( independent  of  stress  level ) 


(2.28a) 

(2.28b) 

(2.28c) 


where  e'o,  m',  and  an,  are  material  constants  and  may  be  functions 

of  the  testing  variables.  Substitution  of  equations  2.28  a, 
b,  and  c  into  equation  2.18  yield  the  Findley  nonlinear  model 
written  as: 


£(t): 


+  m'sinh 

t". 


(2.29) 
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Findley  et  al  [18,21]  and  others  [14,73]  have  used  this 
empirical  fitting  approach  to  successfully  model  various 
material  systems.  An  advantage  of  the  Findley  model  is  that 
the  five  material  constants  (e'o,  m',  Oe,  an,,  n)  can  be  solved  for 
using  creep  test  strain  measurements  without  the  recovery 
measurements . 

Gramoll  et  al  [25]  noted  a  problem  with  the  Findley  nonlinear 
model  at  high  stress  levels.  The  problem  involved  numerical 
difficulties  when  the  sinh  function  rapidly  increases  for 
values  larger  than  1.  Gramoll  solved  the  numerical  problem 
by  using  an  alternative  model  which  he  referred  to  as  the 
"Quadratic  Power  Law".  This  model  replaced  the  sinh  functions 
with  quadratic  functions  such  that; 


eo  =  e'oSinh  — 


m 


/  (  ^ 
=  m  sinh  — 


=  E'o(l+gO^) 


=  m'(l+/o^). 


(2.30a) 

(2.30fc) 


The  constants  g  and  /  are  determined  by  a  least-squared  fit 

of  Eq  and  m  data  from  a  series  of  creep  tests,  at  various  stress 
levels,  to  the  quadratic  function. 

2 . 1 . 3 . 2  Schapexy  Model 

Another  popular  nonlinear  viscoelastic  model  is  the 
Schapery  integral  model  derived  from  irreversible  thermodynamic 
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theory  [54,56].  The  Schapery  model  has  been  successfully 
applied  to  many  materials  [3,37,45].  Work  at  VPI  &  SU  has 
also  used  the  Schapery  model.  Caplan  et  al  [7]  modeled  the 
nonlinear  creep  response  of  polycarbonate  material.  Kiel  et 
al  [27]  characterized  T300/934  graphite-epoxy  composites. 
Tuttle  et  al  [68]  studied  the  nonlinear  effects  in  T300/5208 
graphite-epoxy  composites. 

The  Schapery  model's  popularity  is  due  to  its  single 
time-integral  form  to  represent  nonlinear  constitutive  rela¬ 
tionships.  The  constitutive  equation  for  a  isotropic  material 
under  uniaxial  loading  with  isothermal  conditions  can  be  written 
as : 

with 

_  r  dt' 

where  Dq  and  AD(\|/)  are  the  initial  and  transient  components  of 

compliance,  gotgi^gi  and  a,,  are  stress  dependent  material  func¬ 
tions  - 


(2.31) 

(2.32a) 

(2.32b) 
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The  strain  response  for  a  uniaxial  creep/creep  recovery 
test  can  be  derived  by  integrating  equation  2.31.  The 
creep/creep  recovery  stress  input,  shown  in  figure  2.1,  can 
be  written  as: 


ao[H(t)-H(T-t,)] 


(2.33) 


where  H(x)  is  the  Heaviside  step  function.  Schapery  [55]  notes 
the  time  dependent  creep  compliance  function  can  be  represented 
with  a  power  law  function  as: 

AD(\k)  =  D,v".  (2.34) 


Substituting  equations  2.33  and  2.34  into  equation  2.31  and 
noting  that  is  zero,  except  at  T  =  0  and  T  =  ti,  and  ^^  =  5(t), 
yields  the  creep  strain  response  for  0<t<t, 


ec(0  = 


glg2^l^ 


aS 


(2.35) 


and  the  creep  recovery  strain  response  for  t>t, 


e,(t)  =— [(l+a„X)"-(a„X)"] 
gi 


(2.36) 
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where 


(t-ti) 


Aec=ec(t,)-goDoO 


8182  Pi 

aS 


(2.37) 

(2.38) 


The  strain  Ae^  represents  the  transient  creep  component 

immediately  before  the  stress  is  removed. 

Several  methods  have  been  used  to  solve  for  the  parameters 
in  the  Schapery  model.  Original  procedures,  developed  by 
Schapery  [55]  and  Lou  et  al  [37],  were  based  on  graphical 
methods  requiring  horizontal  and  vertical  shifting  of  linear 
and  nonlinear  data  to  find  a  best  fit  of  the  parameters. 
Bertolotti  [4]  developed  a  computer-based  solution  using  the 
Levenberg-Marquardt  algorithm  to  find  the  model  parameters. 
Similarly,  Augl  et  al  [2]  obtained  the  parameters  using  a 
least-squares  fit  of  experimental  data  to  the  nonlinear  power 
law  formulation  of  the  Schapery  model. 

Schapery  [55]  generalized  his  nonlinear  model  (equation 
2.35)  to  the  Findley  nonlinear  model  (equation  2.29)  by  noting 
the  following  relationships  between  constants: 
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(2.39a) 


(2.39b) 

(2.39c) 
(2.39d) 

Additionally,  when  the  power  law  compliance  function  (equation 
2-34)  is  used  in  the  Schapery  model  and  the  stresses  and  strains 
are  in  the  linear  material  range,  the  Schapery  constants  become 
^0  =  ^1  =  ^2  =  a«=  1/  with  equation  2.35  reducing  to: 

e,(t)  =  DoOo+D,aot"  (2.40) 

which  is  a  form  of  the  Findley  linear  power  law. 

2.2  Physical  Aging  Effects 

This  section  describes  various  subjects  regarding  physical 
aging  and  its  effects  on  the  viscoelastic  creep  response  of 
amorphous  polymers.  Included  are  discussions  explaining  what 
physical  aging  is  and  how  it  is  related  to  the  free  volume 
state  in  the  material.  Also,  the  methods  used  to  characterize 
the  effect  of  physical  aging  on  creep  compliance  response  of 
materials  are  examined. 


e'o  =  Do(Te 
m'  =  D,o_. 
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2.2.1  Basic  Polymer  Concepts 

The  free  volume  concept  is  frequently  used  to  explain 
physical  aging  of  polymer  materials.  This  concept  is  easily 
understood  with  the  help  of  a  two-dimensional  hard-sphere 
molecule  model  used  by  Zallen  [75]  .  The  free  volume  associated 
with  three  particular  molecules  A,  B,  and  C  are  shown  as  shaded 
areas  in  figure  2.9.  The  shaded  areas  define  the  region 
accessible  to  the  mass  center  of  the  molecule  when  constrained 
by  the  repulsive  forces  of  the  neighboring  molecules .  Molecules 
B  and  C  cannot  move  to  locations  where  they  can  be  influenced 
by  different  nearest  neighbor  molecules.  They  are  therefore 
confined  to  small  oscillations  within  their  free  volume  regions. 
The  small  oscillations  are  associated  with  secondary  relax¬ 
ations  in  the  material.  Molecule  A,  with  its  larger  amount 
of  free  volume,  is  capable  of  changing  its  nearest  neighbors 
by  translating  to  a  new  location,  labeled  A'  in  figure  2.9. 
The  ability  of  molecules  to  translate,  commonly  referred  to 
as  mobility,  is  controlled  by  the  amount  of  free  volume  in  the 
material.  A  qualitative  relationship  between  mobility  and 
free  volume  is  shown  in  figure  2.10.  Note,  that  at  high 
values  of  free  volume,  the  degree  of  packing  of  the  molecules 
is  small,  and  the  mobility  of  the  molecules  is  high.  The 
mobility  of  the  molecules  effects  many  of  the  properties  of 
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the  material.  As  the  mobility  increases,  the  material  becomes 
more  compliant  and  ductile.  Also,  the  damping  ability  and 
creep  rate  of  the  material  increase. 

Another  important  concept  is  the  distinction  between 
thermoplastic  and  thermoset  classes  of  polymer  composites. 
Thermoplastic  polymers  are  characterized  by  their  thermal 
reversibility.  The  material  will  soften  upon  heating  and  flow 
when  a  stress  is  applied.  When  cooled,  the  material  reversibly 
regains  its  solid  or  rubbery  form.  Quite  differently,  thermoset 
materials  undergo  an  irreversible  curing  reaction  upon  heating 
and  will  not  regain  their  original  form  upon  cooling.  During 
the  curing  process,  the  polymer  chains  cross-link  with  strong 
covalent  bonds  that  strengthen  the  material  and  increase  the 
stiffness . 

A  schematic  of  master  creep  compliance  curves  for  typical 
thermoplastic  and  thermoset  polymers  are  shown  in  figure  2.11. 
The  basic  difference  between  the  two  classes  of  polymers  is 
their  behavior  above  T,.  Thermoplastic  materials  melt,  whereas 
thermoset  materials  char  above  T,.  Master  curves  are  con¬ 
structed  by  shifting  compliance  data,  from  tests  at  various 
temperatures,  along  the  log  time  axis  using  TTSP.  Test  segments 
are  arranged  with  lower  temperature  tests  on  the  left  side  and 
higher  temperature  tests  shifted  towards  the  right  side  of  the 
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figure  2.11.  In  this  study  the  region  of  particular  interest 
is  between  the  P  transition  and  the  a  transition  temperatures 
and  is  referred  to  as  the  glassy  region. 

The  P  transition  defines  the  highest  secondary  transition 
temperature.  Below  this  temperature  side  groups  or  parts  of 
chain  segments  can  rotate  and  vibrate,  but  the  available  thermal 
energy  is  too  small  to  allow  entire  chain  segment  to  move. 
In  the  region  above  the  P  transition  temperature  and  below  the 
a  transition  temperature,  the  motion  of  entire  chain  segment 
is  possible  but  is  limited  by  the  amount  of  free  volume  available 
for  the  segments  to  move  into.  Above  the  a  transition  tem¬ 
perature,  commonly  referred  to  as  the  glass  transition  tem¬ 
perature  (Tg) ,  the  free  volume  is  directly  related  to 
temperature.  The  thermal  energy  available  above  Tg  is  large 
enough  to  allow  the  polymer  chain  segments  to  move  easily. 
Thus,  it  is  possible  to  reach  the  equilibrium  free  volume  state 
very  rapidly  at  temperatures  above  Tg. 

The  matrix  material  of  the  composite  analyzed  during  this 
study  is  classified  as  an  amorphous  thermoplastic  polymer. 
The  amorphous  type  of  polymer  is  characterized  by  a  random 
non-ordered  arrangement  of  polymer  chains.  When  the  polymer 
chains  have  a  regular  ordered  chain  structure  packed  into  a 
crystal  lattice,  somewhere  in  the  material,  it  is  classified 
as  a  semicrystalline  polymer.  The  degree  of  crystallinity 
is  often  given  as  a  percentage  of  amorphous  structure  that  has 
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formed  into  the  crystal  lattice  structure.  Struik 
[60.61,64,  65],  Chai  et  al  [8]  and  Ogale  et  al  [48]  have  studied 
how  crystallinity  effects  physical  aging  in  semicrystalline 
materials . 

2.2.2  Physical  Aging  Definition 

The  term  aging  is  used  to  describe  the  process  that  changes 
the  properties  of  a  material  with  time.  Aging  processes  are 
generally  grouped  into  either  a  physical  aging  or  chemical 
aging  category.  Physical  aging  processes  are  characterized 
by  the  thermoreversibility  of  the  change  in  the  material 
properties.  The  material  can  be  rejuvenated  by  heating  it 
above  Tg  and  holding  there  until  a  state  of  equilibrium  free 
volume  is  reached.  Conversely,  chemical  aging  is  characterized 
by  the  irreversibly  of  the  change  in  material  properties. 
Chemical  aging  processes,  such  as  oxidation,  curing,  and 
cross-linking,  permanently  change  the  material  properties. 

Early  work  by  Kovacs  et  al  [32]  showed  that  amorphous 
solids  fail  to  reach  their  equilibrium  free  volume  state  when 
they  are  rapidly  cooled  from  above  their  T,.  The  equilibrium 
free  volume  state  is  achieved  by  cooling  the  material  at  an 
infinitely  slow  cooling  rate.  Physical  aging  is  the  process, 
occurring  below  T^,  where  the  raacromolecules  gradually  change 
their  packing  in  order  to  approach  the  equilibrium  free  volume 
state.  The  origin  of  physical  aging  is  illustrated  in  figure 
2.12.  As  the  material  is  cooled  from  above  Tg,  the  free  volume 
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linearly  decreases  with  temperature.  Also,  the  polymer  chain 
mobility  decreases  following  its  relationship  with  free  volume, 
shown  in  figure  2.11.  At  T,,  the  viscosity  rapidly  increases 
and  the  polymer  chain  motion  is  impaired.  The  polymer  is 
unable  to  reach  the  equilibrium  degree  of  packing  before 
solidification.  Continued  cooling  creates  a  free  volume 
differential  between  the  quenched  and  equilibrium  states.  As 
the  polymer  undergoes  isothermal  aging,  the  free  volume 
gradually  decreases  towards  the  equilibrium  state. 

The  process  controlling  the  isothermal  free  volume 
relaxation  is  a  nonlinear  closed-loop  scheme  symbolically 
written  as: 

(2.41) 


Equation  2.41  implies  that  the  initial  free  volume  (Vf) 

determines  the  mobility  (M)  of  the  polymer  segments.  The 
mobility  available  to  the  polymer  segments  then  determines  the 
time  rate  of  change  of  the  free  volume.  The  whole  process 
starts  after  the  quench  from  above  T,  when  the  initial  mobility 
prompts  a  change  in  free  volume  towards  the  equilibrium  value. 
The  decrease  in  free  volume  forces  the  mobility  to  also  decrease. 
This  is  followed  by  a  decrease  in  the  rate  that  the  free  volume 
relaxes  from  the  material.  An  interesting  consequence  of 
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equation  2.41  is  that  below  T,  the  mobility  can  only  asymp¬ 
totically  approach  zero.  This  is  explained  by  noting  that  a 
decrease  in  mobility  requires  a  decrease  in  free  volume.  In 
order  for  the  free  volume  to  change,  there  must  be  some  mobility. 
Therefore,  the  mobility  cannot  disappear  over  a  finite  time 
period  but  may  become  very  small  as  the  material  approaches 
the  equilibrium  free  volume  value.  As  the  free  volume  and 
mobility  decrease,  all  the  properties  of  the  material  that 
depend  on  them  will  also  change. 

2.2.3  Effect  of  Physical  Aging  on  Material  Properties 

Although  this  study  concentrates  on  the  effect  of  physical 
aging  on  the  creep  properties  of  a  material,  many  other  material 
properties  also  depend  on  mobility  or  free  volume  state  in  the 
material.  An  excellent  review  of  the  material  properties  of 
polymer  glasses  that  are  effected  by  the  nonequilibrium  volume 
behavior  is  presented  by  Tant  and  Wilkes  [67].  Therefore, 
only  a  brief  summary  is  given  here.  Associated  with  a  decrease 
in  mobility  and  free  volume  are  increases  in  density,  viscosity, 
yield  stress,  and  elastic  modulus.  An  initial  assessment  might 
suggest  that  the  increases  in  yield  stress  and  stiffness,  due 
to  physical  aging,  would  be  beneficial  during  the  life  of  an 
aerospace  structure.  Unfortunately,  physical  aging  also 
decreases  the  impact  strength  and  fracture  energy  of  the 
material.  Degradation  of  these  properties  must  be  considered 
in  aerospace  applications  susceptible  to  impact  type  damage. 
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2.2.4  Effect  of  Physical  Jlging  on  Creep  Compliance 

The  effect  of  physical  aging  on  creep  compliance  of  a 
material  appears  as  a  horizontal  shift  of  the  compliance  data 
on  a  logarithmic  time  scale.  An  illustration  of  the  aging 
effect  on  amorphous  polymer  poly  (vinyl  chloride)  (PVC)  obtained 
from  Struik  [59],  is  shown  in  figure  2.13.  The  PVC  material 
was  quenched  from  90°C  to  20°C  and  tested  at  aging  times  ranging 
form  0.03  to  1000  days.  The  creep  compliance  data  associated 
with  each  aging  time  was  horizontally  shifted  to  the  1000  day 
aging  time  to  form  a  master  compliance  curve.  The  excellent 
superimposability  of  the  individual  creep  curves  shows  that 
the  shape  of  the  curve  remains  constant  during  aging.  It  is 
also  important  to  note  that  the  individual  creep  curves,  at 
various  aging  times,  extend  over  several  decades  in  time.  This 
implies  that  the  creep  process  requires  a  wide  distribution 
of  relaxation  times  instead  of  a  single  relaxation  time. 
However,  since  the  shape  of  the  creep  curve  remains  constant, 
all  the  relaxation  times  associated  with  a  particular  aging 
time  must  be  shifted  by  the  same  factor,  usually  referred  to 
as  a.  Note,  henceforth  the  term  relaxation  time,  as  adopted 
by  Struik  [59],  is  used  generically  in  the  description  of  the 
characteristic  time  associated  with  either  stress  relaxation 
or  creep  mechanical  models. 
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Also  illustrated  in  figure  2.13  is  the  thermoreversibility 
of  physical  aging.  After  the  1000  days  of  aging,  the  sample 
was  heated  back  to  90°C  for  20  minutes  and  quenched  to  the  20°C 
aging  temperature.  The  creep  compliance  was  measured  again 
after  aging  one  day  (data  represented  by  x  symbol  in  figure 
2.13)  was  within  2%  of  the  original  measurements  made  nearly 
four  years  earlier.  Thus,  previous  aging  can  be  completely 
erased  by  heating  the  material  above  Tg  for  a  short  period  of 
time . 

The  horizontal  shift,  on  the  logarithmic  time  scale,  between 
successive  creep  curves  are  approximately  constant,  as  shown 
in  figure  2.13.  This  implies  there  is  a  relationship  between 
the  amount  of  horizontal  shift  (loga)  and  the  log  of  the  aging 
time  OogfJ.  Using  Struik's  [59]  sign  convention,  loga  is  a 
measurement  of  the  horizontal  shifting  of  the  referenced 
compliance  curve  and  is  taken  positive  for  a  shift  to  the  left. 
The  aging  can  be  characterized  by  the  double-logarithmic  shift 
rate  defined  as; 


^(loga) 


(2.42) 


where  -loga  defines  the  amount  of  horizontal  shift  required  to 
move  the  individual  compliance  curve  from  its  original  location 
to  the  shifted  position  required  to  form  the  master  momentary 
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compliance  curve  at  the  referenced  aging  time.  The  horizontal 
shift  values  used  to  shift  the  r,=l  hour  reference  aging  curve 
to  the  aging  times  shown  in  figure  2.13  are  plotted  verses 
logt,  in  figure  2.14.  The  slope  of  this  curve  is  the  shift 
rate  ^i.  It  was  calculated  at  about  unity  for  the  PVC  material. 

The  effect  of  temperature  on  the  shift  rate,  \i,  for  several 
polymer  materials,  is  shown  in  figure  2.15.  Each  material  in 
the  figure  shows  a  rapid  increase  in  |i  below  Tg  to  a  value  near 
unity.  As  the  temperature  continues  to  decrease,  p.  reaches  a 
plateau  near  unity  for  a  range  of  temperatures  and  then  decreases 
at  a  lower  temperature  which  varies  with  each  material.  The 
changes  in  |i  with  temperature  can  be  explained  using  free 
volume  theory.  Struik  [59]  describes  the  decrease  in  mobility 
associated  with  physical  aging  as  a  "self-retarding"  process. 
As  the  mobility  decreases  during  aging,  following  the  rela¬ 
tionship  given  by  equation  2.41,  the  rate  of  decrease  of 
mobility  will  also  decrease.  Struik  [59]  suggested  that  the 
mobility  associated  with  self-retarding  processes  is  inversely 
proportional  to  the  aging  time: 


M--  (2.43) 

where  is  a  constant.  Since,  the  retardation  times  associated 
with  creep  tests  are  also  inversely  proportional  to  mobility. 
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equation  2.43  indicates  that  they  will  also  increase  propor¬ 
tionally  to  the  aging  time.  Therefore,  an  increase  in  by  a 
factor  of  ten  will  also  increase  the  retardation  times  by  a 
factor  of  ten.  This  produces  a  shift  of  the  creep  curve  along 
the  logarithmic  time  axis  by  one  decade.  This  explains  why 
the  shift  rate,  |i,  is  nearly  constant  in  the  aging  temperature 
range  and  why  it  is  approximately  unity.  The  shift  rate 
decreases  to  zero  as  the  equilibrium  volume  state  is  approached 
at  temperatures  near  Tg.  At  Tg,  volume  equilibrium  is  reached 
and  physical  aging  ceases  with  ^=0.  The  aging  rate  decreases 
at  lower  temperatures  as  the  material  approaches  the  highest 
secondary  transition  temperature,  Tp.  Below  this  temperature, 
only  secondary  relaxations  exist  which  do  not  change  the  free 
volume  state  in  the  material. 

2.2.5  Characterizing  and  Predicting  the  Effect  of  Physical 
Aging 

In  this  section  the  methodology  used  to  characterize  the 
effect  of  physical  aging  on  linear  small-strain  creep  compliance 
is  discussed.  Also,  a  method  of  predicting  long-term  creep 
behavior  form  short-term  aging  characterization  tests  is 
presented  for  a  single  compliance  direction.  Further  expla¬ 
nation  of  the  process  required  to  characterize  the  compliance 
response  in  each  of  the  principal  directions  of  a  two  dimensional 
orthotropic  composite  material  is  contained  in  section  2.3. 
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2 . 2 . 5 . 1  Momentary  Compliance  Curves 

This  section  describes  the  differences  between  short-term 
(momentary)  and  long-term  creep  tests.  Momentary  tests  are 
essentially  a  "snapshot"  of  the  material  properties  taken  at 
a  constant  aging  time.  The  duration  of  the  momentary  test 
remains  small  compared  to  the  aging  time  at  the  beginning  of 
the  test.  In  this  study,  the  relationship  between  total  test 
duration  (t,)  and  aging  time  (t^)  of  was  used.  In 
contrast,  long-term  tests  have  aging  occurring  simultaneously 
with  the  creep.  The  test  duration,  therefore,  becomes  much 
larger  than  the  aging  time  at  the  beginning  of  the  test. 

During  a  short-term  test,  the  retardation  times  hardly 
change  due  to  aging  and  the  resulting  compliance  curve  is 
referred  to  as  a  momentary  curve.  A  schematic  of  a  sequence 
of  momentary  creep  and  recovery  tests  is  shown  in  figure  2.16. 
After  the  temperature  is  quenched  from  above  T,  to  a  constant 
aging  temperature  below  Tg,  shown  as  the  dashed  line  in  figure 
2.16,  the  aging  time  starts  for  the  entire  aging  test  sequence. 
Next,  successive  creep  tests  are  conducted  at  various  aging 
times.  The  recovery  strain  built  up  from  previous  creep  tests 
is  extrapolated  into  the  next  creep  period  and  subtracted  from 
the  total  strain  response.  Thus,  the  creep  strain  associated 
with  each  aging  test  can  be  singled  out.  The  extrapolation 
is  reliable  since  the  creep  test  duration  is  less  than  20%  of 
the  previous  recovery  period  time  scale. 
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The  creep  test  sequencing  procedure  can  only  be  used  with 
small  strain  levels  in  the  linear  viscoelastic  behavior  range 
where  Boltzmann's  Superposition  Principle  is  valid.  The 
validity  of  the  sequencing  procedure  can  be  checked  by  comparing 
the  compliance  results  of  an  intermediate  aging  time  creep 
test  from  a  sequence  test  with  a  second  test.  The  second  test 
measures  the  creep  response  at  the  same  intermediate  aging 
time,  but  eliminates  the  creep  tests  conducted  prior  to  the 
intermediate  aging  time.  A  reasonable  comparison  of  the 
compliance  between  the  two  tests  verifies  the  applicability 
of  the  sequencing  test  procedure.  Creep  tests  outside  the 
linear  range  require  a  more  laborious  procedure  which  requires 
a  new  quench  before  each  creep  test. 

2. 2. 5. 2  Aging  Shift  Rate  Determination 

The  aging  shift  rate,  |l,  was  determined  from  the  results 
of  a  sequence  of  aging  tests  conducted  at  a  constant  temperature. 
Struik  [59],  states  that  the  amount  of  horizontal  and  vertical 
shifting  required  to  make  a  master  compliance  curve  can  be 
determined  unambiguously,  from  the  marked  curvature  of  the 
creep  curves.  In  practice,  matching  the  curvatures  of  the 
momentary  creep  curves  is  a  precarious  process.  Historically, 
shift  factors  have  been  determined  by  shifting  creep  curves 
by  hand.  This  process  is  a  time  consuming  and  tedious  job 
utilizing  a  trial  and  error  procedure.  In  this  study,  a 
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computer  program  was  written  to  automate  the  shifting  process. 
The  program  is  based  on  a  procedure  described  by  Vleeshouwers 
et  al  [72]  where  the  Levenberg-Marquardt  [28]  nonlinear  finite 
difference  scheme  was  used  to  fit  each  momentary  creep  com¬ 
pliance  curve  to  Struik's  compliance  model  given  by  equation 
2.19  as: 


D(t)  =  Doe 


(‘/lo)* 


(2.19) 


The  program  determines  the  values  of  Do,  to,  and  m  that  best 

fit  the  test  data.  The  variables  controlling  the  amount  of 
horizontal  and  vertical  shift,  on  a  logarithmic  scale,  become 
apparent  after  taking  the  log  of  both  sides  of  equation  2.19, 
yielding: 


logD(t)  =  logDo+{t/to)”]og^.  (2.44) 

Clearly,  Dq  controls  the  vertical  shift,  to  controls  the  hor¬ 
izontal  shift,  and  m  defines  the  shape  of  the  curve  on  a  log 
compliance  vs.  log  time  plot.  Both  Dp  and  to  are  functions  of 
temperature  and  aging  time,  whereas  m  is  a  constant  that 
defines  the  shape  of  the  compliance  curve. 
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The  first  pass  of  the  computer  program  determines  the 
values  of  Dq,  to,  and  m  that  best  fit  the  compliance  data  for 
each  aging  time.  The  m  values  of  all  the  compliance  curves 
from  a  test  sequence  are  compared,  and  an  average  m  value  is 
calculated.  The  second  pass  of  the  computer  recalculates  Dq 
and  to  with  the  m  value  forced  to  the  average  value.  The  amount 
of  horizontal  shift  ,  — loga,  required  to  move  each  compliance 
curve  to  the  referenced  curve  is  computed  using: 


-  log  a  =  log 


to(t«,r<r)  j 


(2.45) 


An  illustration  of  the  shifting  of  momentary  creep  compliance 
curves  is  shown  in  figure  2.17.  The  individual  compliance 
curves  in  figure  2.17  are  obtained  at  different  aging  times 
during  a  sequence  of  creep  tests  conducted  at  a  constant 
temperature  and  stress  level.  Note  that  the  last  aging  time 
in  the  test  sequence  was  chosen  as  the  referenced  creep  curve. 
The  amount  of  horizontal  shift,  -loga,  required  to  shift  each 
aging  test  to  the  referenced  creep  curve  are  also  shown  on 
figure  2.17.  Finally,  the  shift  rate,  \i,  was  determined  using 
a  linear  regression  routine  to  determine  the  slope  of  a  plot 
of  —loga  vs.  logt,,  shown  inserted  in  the  upper  left  corner  of 
figure  2.17.  The  aging  rate  calculation  was  found  to  be 
insensitive  to  the  average  value  of  m  used  to  fit  the  aging 
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curves  in  a  test  sequence.  The  sensitivity  was  checked  by 
varying  m  average  value  from  .34  to  .48  by  .02  increments  and 
calculating  the  aging  rate  for  each  value.  The  resulting  aging 
rates  were  all  within  ±0.6%  of  the  mean  aging  rate  calculated 
for  the  eight  m  values. 

The  small  amount  of  vertical  shifting  on  the  logarithmic 
scale  was  required  to  bring  the  various  aging  creep  compliance 
curves  into  alignment  with  the  referenced  aging  curve.  The 
amount  of  vertical  shift,  logB,  was  computed  using: 

log  B  =  log  -  ■  .  (2.46) 

The  small  vertical  shift  accounts  for  in  increase  in  initial 
compliance  value  as  the  material  physically  ages.  The  small 
vertical  shift  factors  used  in  the  data  reduction  procedure 
are  tabulated  in  chapter  4.0.  Additionally,  vertical  shift 
factors  are  used  to  correct  for  the  change  in  initial  compliance 
with  temperature.  Sullivan  [66]  reported  that  the  vertical 
shift  factors  account  for  the  increase  in  creep  compliance 
found  with  increasing  temperature. 

The  momentary  compliance  curve  resulting  from  the  shifting 
procedure  is  referenced  to  the  test  temperature  and  last  aging 
time  in  the  test  sequence.  The  effect  of  temperature  on  the 
shift  rate  is  determined  by  repeating  the  entire  creep  test 
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sequence  and  shifting  the  resulting  momentary  creep  curves  to 
determine  the  shift  rate  at  the  new  temperature.  The  number 
of  time  decades  the  compliance  curve  covers  is  limited  to  the 
duration  of  the  longest  aging  test.  In  order  expand  the  time 
range,  use  of  the  time-temperature  superposition  principle  is 
required. 

2. 2. 5. 3  TTSP  Master  Momentary  Compliance  Curve 

TTSP  Master  momentary  compliance  curves  were  formed  by 
shifting  momentary  compliance  curves,  obtained  at  different 
temperatures  with  the  same  aging  time,  using  the  time- 
temperature  superposition  principle.  An  illustration  of  the 
shifting  of  momentary  compliance  curves  to  form  a  TTSP  master 
momentary  compliance  curve  is  shown  in  figure  2.8.  The  master 
curve  in  this  figure  is  referenced  to  the  lowest  temperature 
at  the  referenced  aging  time  The  amount  of  horizontal 
shift,  -logat,  for  each  compliance  curve  is  a  function  of  both 
temperature  and  the  referenced  aging  time. 

Two  methods  were  used  to  determine  the  shifting  factors 
associated  with  temperature.  The  first  method  used  the  same 
computer  program  and  procedures  described  in  the  previous 
section  to  fit  each  compliance  curve  to  Struik's  equation. 
This  method  requires  the  calculation  of  an  average  m  value 
from  the  short-term  momentary  compliance  curves  at  various 
temperatures.  The  horizontal  shift  factors  were  calculated 
using: 
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f 


(2.47) 


-log  3,.  =  log 


toCr,^,  ^ 


The  calculated  shift  factors  are  referenced  to  a  particular 
reference  temperature  aging  time  and  the  average  m 
value.  Also,  using  this  method,  the  average  m  value  chosen 
predetermines  the  shape  of  the  master  compliance  curve.  The 
average  m  value  was  calculated  from  compliance  data  with  a 
limited  time  range  (~  3  decades)  and  may  or  may  not  be  the 
best  representation  of  the  actual  shape  of  the  master  compliance 
curve . 

The  second  method  utilized  a  graphical  procedure  to 
determine  the  shift  factors.  The  compliance  curves,  from 
various  temperatures  at  a  constant  aging  time,  were  all  plotted 
using  the  Freelance  graphical  package  by  Lotus  [24] .  The 
master  curve  was  generated  by  matching  the  slopes  of  the 
compliance  curves.  The  ruler  option  in  Freelance  was  utilized 
to  scale  the  x  and  y  axis  in  logarithmic  units.  The  shift 
factors  for  each  curve  were  easily  determined  by  subtracting 
the  final  curve  position  from  the  initial  position.  The  shifted 
compliance  data  was  then  fit  to  Struik's  equation  using  the 
nonlinear  finite  difference  program.  Comparisons  between  the 
two  master  curve  construction  methods  and  their  effect  on 
momentary  and  long-term  compliance  predictions  are  discussed 
in  section  4.0  and  5.0  respectively. 
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2. 2. 5. 4  Effttctiva  Tima  Theory 

In  this  section,  the  effective  time  theory  developed  by 
Struik  [59]  is  examined.  The  effective  time  theory  was  developed 
to  account  for  the  effect  of  aging  during  long-term  creep 
compliance  tests.  Although  time-temperature  superposition 
works  well  for  momentary  creep  compliance  curves  at  a  constant 
aging  state,  it  overestimates  the  creep  response  during 
long-term  tests.  Thus,  it  is  necessary  to  use  the  effective 
time  theory  to  account  for  the  gradual  stiffening  of  the 
material  during  physical  aging. 

In  order  to  avoid  confusion,  a  few  time  variables  associated 
with  long-term  creep  tests  will  be  defined.  First,  the  aging 
time  at  the  start  of  the  creep  test,  f,,  is  measured  starting 
after  the  quench  from  above  T,.  The  variable,  t,  is  used  to 
define  the  time  after  the  creep  load  is  applied.  Finally,  the 
aging  time  during  the  creep  test,  t',,  is  calculated  by  adding 
t  to  .  Each  of  these  time  variables  is  illustrated  in  figure 
2.18. 

Development  of  the  effective  time  theory  begins  by  noting 
that  the  shift  rate  ^  remains  constant  for  very  long  periods 
of  time.  Struik  [41]  suggests  that  the  shift  rate,  determined 
from  short-term  tests,  should  be  close  to  unity  (e.g.,  0.7  < 
^  <  1.0)  in  order  to  obtain  reliable  long-term  predictions. 
The  shift  rate  is  determined  using: 
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(2.48) 


^loga  loga,,-Ioga,^ 
^  dlogr.  log  r'.- log  f. 


where  a  is  the  acceleration  factor.  The  master  momentary  creep 
curve  is  referenced  to  the  aging  time  at  the  start  of  the 
long-term  test.  Thus,  the  value  of  loga,^  term  will  be  zero  at 
r, .  Substituting  this  relationship  into  equation  2.4  8  yields: 

=  loga,..  (2.49) 

Raising  both  sides  of  equation  2.49  to  the  tenth  power  yields: 


a(r'.,/.)  = 


(2.50) 


Abbreviating  the  acceleration  factor  as  a(t)  and  writing  it  in 
terms  of  t  and  yields: 


a(t)  = 


r,+t 


Y 


(2.51) 


At  aging  times  greater  than  r,,  the  effect  of  physical  aging 
shifts  the  referenced  momentary  compliance  curve  in  the  negative 
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(to  the  right  using  Struik' s  sign  convention)  direction.  This 
has  the  effect  of  extending  the  momentary  creep  compliance 
curve  along  the  time  axis  to  form  the  long-term  compliance 
curve  as  illustrated  in  figure  2.18.  The  momentary  compliance, 
D,^(X),  is  related  to  the  long-term  compliance,  D,^(t),  using  an 
effective  time  X.  The  relationship  is  given  by: 

D,/t)  =  DjX).  (2.52) 

Both  compliance  representations  are  referenced  to  the  aging 
time,  at  the  start  of  the  creep  test.  As  the  creep  test 
progresses  from  time  t  to  t+dt,  the  relaxation  processes  are 
all  slowed  down  by  the  same  amount.  The  time  interval  dt  is 
related  to  the  effective  time  interval,  dX,  by: 

ca  =  a(t)dt.  (2.53) 

Integration  of  equation  2.53  yields  the  effective  time  X  given 
by: 


X(t)=  f‘a(^)d^. 

Jo 


(2.54) 
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where  ^  is  a  dummy  integration  variable  on  the  t  time  scale. 
An  analytical  form  for  X  is  obtained  by  combining  equations 
2.51  and  2.54  and  performing  the  integration.  This  yields 
the  effective  time  expression  used  by  Struik  [59]  given  as: 


where 


for  ^1=1 


for  |1  <  1 


0<a=  1  -^i<  1 


(2.55) 


The  two  variables  that  influence  the  effective  time  and 
therefore  the  long-term  compliance  predictions  are  the  aging 
rate  [L,  determined  from  short-term  tests,  and  the  aging  time 
r,,  at  the  start  of  the  long-term  test.  The  long-term  compliance 
predictions,  associated  with  the  real  time  axis,  are  made  by 
substituting  the  effective  time,  computed  using  equation  2.55, 
into : 


D(t)  =  Doe 


(2.56) 
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In  addition  to  the  dependence  on  effective  time,  the  long-term 
compliance  prediction  also  depends  on  how  well  the  variables 
Do,  to,  and  m  define  the  momentary  master  compliance  curve  at 
.  Therefore,  obtaining  a  correct  representation  of  the 
momentary  creep  response  of  the  material  using  TTSP  is  essential 
for  accurate  long-term  predictions. 

2.3  Composite  Material  Compliance  Relationships 

The  composite  material  system  of  interest  in  this  study 
consisted  of  Radel  X  matrix  material  surrounding  IM7  fibers. 
The  fibers  were  embedded  in  the  matrix  to  form  a  thin  lamina 
sheet  approximately  0.005  inch  (0.013  cm)  thick.  The  thin 
unidirectional  lamina  sheets  were  stacked,  with  the  fibers  all 
orientated  in  the  same  direction,  to  form  a  12-ply  unidirec¬ 
tional  laminate.  After  curing,  the  completed  unidirectional 
laminate  panel  is  classified  as  an  orthotropic  material. 
Therefore,  nine  independent  material  properties  are  required 
for  full  characterization.  However,  since  the  thickness  of 
the  panel  is  quite  small  in  comparison  to  the  other  dimensions 
it  is  commonly  classified  as  a  two  dimensional  orthotropic 
material.  Thus,  under  in-plane  loading  (CF33  =  X23  =  Xs,  =  0)  only  four 
independent  properties  are  required. 

The  coordinate  system  used  to  define  the  directionality 
of  the  material  properties  followed  a  convention  used  by  Jones 
[29]  and  is  illustrated  in  figure  2.19.  In  this  system,  the 
loading  or  global  axis  is  the  (X-Y)  axis  and  the  material  or 
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local  axis  is  referred  to  as  the  (1-2)  axis.  The  angle  between 
the  coordinate  systems  is  referred  to  as  0  and  is  measured  in 
a  positive  direction  from  the  X  axis  to  the  1  axis.  The  four 
material  properties  associated  with  creep  testing  are  defined 
in  terms  of  the  compliance  matrix  in  the  material  principal 
directions.  The  elastic  constitutive  equations  relating 
stress,  strain,  and  compliance  for  the  two  dimensional  ort¬ 
hotropic  unidirectional  laminate  are: 


El' 

S,2 

O' 

£2 

'  = 

*^21 

S'22 

0 

* 

yn 

0 

0 

V 

(2.57) 


This  matrix  equation,  written  in  condensed  notation,  is  given 
by: 


{e}u  =  [51,2{o},2  (2.58) 

where 

{e},2  =  in-plane  local  strains 
[5],2  =  principal  compliance  matrix 
{o},2  =  in-plane  local  stresses. 

The  compliance  terms  are  related  to  common  engineering  prop¬ 
erties  by  the  following  equations: 
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Determination  of  the  principal  compliance  terms  requires 
testing  unidirectional  laminates  at  different  angle  orienta¬ 
tions  between  the  local  material  axis  and  the  global  load  axis. 
It  is  therefore  necessary  to  relate  local  and  global  coordinate 
systems  through  transformation  equations.  The  transformations 
between  the  local  (1-2)  axis  and  global  (X-Y)  axis  for  stress 
and  strain  are: 


{o}.2  =  [r,]{a}^  (2.60) 

{e},2  =  [7'J{e}^  (2.61) 

where  the  transformation  matrices  are  defined  as: 
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Also  for  completeness,  the  inverse  transformations  between  the 
global  (X-Y)  axis  and  local  (1-2)  axis  for  stress  and  strain 
are 


=  (2.62) 
{e}.2  (2.63) 


where  the  inverse  transformation  matrices  are  defined  as: 


cos*  6 

sin*0 

-2cos0sin0 

cos*0 

sin*0 

-cos0sin0 

[T,]-'  = 

sin*0 

cos*0 

2cos0sin0 
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sin*0 

cos*0 

cos  0  sin  0 

cos  6  sin  9 

-cos0sin0 

cos*  0- sin*  0 

2  cos  0  sin  0 

-2cos0sin0 

cos*  0- sin*  0 

The  compliance  matrix  in  the  material  principal  direction 
can  also  be  transformed  to  the  global  coordinate  system.  This 
is  accomplished  by  substituting  equations  2.58  into  2.63 
yielding : 


{eU  =  [7'J-‘[5],2{a},2 


(2.64) 


which  is  followed  by  a  substitution  of  {o},2  from  equation  2.60. 
The  result  is: 


(2.65) 
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or  written  in  condensed  notation  as: 


=  (2.66) 

where 

mi2  =  [rj-'[5],,[7,].  (2.67) 


[5]i2  is  referred  to  as  the  transformed  or  rotated  compliance 

matrix.  It  represents  the  material  compliance  in  the  global 
coordinate  directions.  The  expanded  rotated  compliance  terms 
are  computed  by  completing  the  matrix  multiplication  in  equation 
2.67  yielding: 


5,1=5,,  cos^  0  +  (25,2  +  sin*  0  cos*  0  +  522  sin^  0  (2.68u ) 

5,2  =  5,2(5111“  0  +  cos“  0)  +  (5, ,  +  522  -  ^66)  sin^  0  cos*  0  (2.686 ) 

5  22  =  5, ,  sin“  0  +  (25,2  +  sin*  0  cos*  0+522  cos“  0  (2.68c ) 

5,„  =  (25,,  -25,2 -5^4)  sin  0cos^0- (2522 -25,2“  5jg)sin^0  cos  0  i2.6Sd) 

526  =  (25,,  -25,2"  5jg)  sin*  0  cos  0- (2522 -25,2“  5g6)  sin  0cos^0  (2.68e) 

5^6  =  2(25„  +  522  -  45,2  -  ^66)  sin*  0  cos*  0  +  566(sin“  0  +  cos“  0).  (2.68/) 


The  four  principal  compliance  terms  5,„5,2,522,  and  were 

determined  from  the  results  of  creep  tests  conducted  at  known 
aging  times.  The  compliance  was  calculated  by  dividing  the 
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known  stress  level  applied  to  the  specimen  by  the  measured 
strain  response.  The  strain  response  was  determined  using 
foil  resistance  strain  gages  mounted  on  the  specimens.  Further 
details  on  the  use  the  strain  gages  are  contained  in  section 
3.4.  Other  methods  of  strain  measurement  such  as  extensometers, 
Moire  interferometry,  and  laser-interferometry  are  possible, 
but  were  not  used  in  this  study. 

A  specific  orientation  of  the  unidirectional  fibers,  with 
respect  to  the  global  axis,  was  required  in  order  to  determine 
each  principal  compliance  term.  The  5„  compliance  term  was 
obtained  using  a  specimen  with  the  fibers  aligned  along  the 
global  X  axis.  This  orientation  is  referred  to  as  a  0  degree 
specimen.  S„  was  calculated  by  substituting  the  axial  strain 
measured  in  the  specimen  into  the  equation: 


The  S22  compliance  term  was  determined  using  a  specimen  with 

the  fibers  oriented  perpendicular  to  the  global  X  axis.  This 
orientation  is  referred  to  as  a  90  degree  specimen.  S22  was 
calculated  by  substituting  the  axial  strain  along  the  load 
axis  into  the  equation: 
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(2.70) 
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The  S12  compliance  term  was  obtained  by  using  the  measured 

strain  in  the  transverse  direction  from  the  0  degree  specimen 
test.  S,2  was  calculated  using  the  equation: 


(2.71) 


The  compliance  tensor  was  assumed  symmetric  such  that  Si2  =  S2t. 

This  assumption  has  been  confirmed  by  other  investigators 
(47, 66]  for  similar  two  dimensional  orthotropic  materials  under 
in-plane  loading.  The  symmetric  assumption  was  not  verified 
in  this  study  due  to  the  difficulty  of  accurately  measuring 
the  transverse  strain  in  a  90  degree  specimen.  The  small  loads 
applied  to  the  90  degree  specimen  produced  very  small  transverse 
strains  in  the  0  degree  direction.  Thus,  it  was  impractical 
to  accurately  determine  52i. 

The  shear  compliance  term,  is  not  as  easily  determined 

as  the  other  compliance  terms.  Several  test  methods  have  been 
proposed  to  determine  Possible  methods  include  the 
rail-shear  test,  picture-frame  test,  and  off-axis  tensile 
specimen  test.  Among  these,  the  off-axis  tensile  specimen 
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test  has  gained  wide  acceptance.  The  10  degree  off-axis  test, 
proposed  by  Chamis  and  Sinclair  [9],  has  been  used  by  several 
investigators  at  VPI  &  SU  [25,68]  to  characterize  the  shear 
compliance  term.  A  disadvantage  of  this  test  is  the  requirement 
for  a  specimen  gage-length  to  gage-width  ratio  of  at  least  14. 
The  relatively  long  and  narrow  specimen  is  required  to  minimize 
the  effect  of  the  end  constraint  produced  by  the  specimen 
grips.  Unfortunately,  the  required  length  for  0.5  inch  (1.27 
cm)  wide  specimens,  used  in  this  study,  was  too  long  to  fit 
in  the  available  ATS  oven.  The  test  configuration  selected 
for  shear  characterization  utilized  a  45  degree  off-axis 
specimen.  Pindera  and  Herakovich  [52]  proposed  this  config¬ 
uration  to  minimize  the  effect  of  end  constraints.  They  found 
that,  for  even  highly  anisotropic  composites,  the  effect  of 
end  constraints  becomes  negligible  for  off-axis  angles  greater 
than  45  degree.  The  shear  compliance.  See,  is  determined  by 
measuring  the  shear  strain  in  the  45  degree  specimen  during 
uniaxial  tension  loading.  The  shear  strain  in  the  principal 
material  axis,  7,2,  is  determined  using  the  strain  readings 
from  a  three  strain  gage  rosette.  The  orientation  of  the 
strain  gage  rosette  on  the  45  degree  specimen  is  shown  in 
figure  2.20.  The  angle  P,  is  defined  as  the  angle  between  the 
X,  global  load  axis,  and  a  line  parallel  to  the  gage  orientation. 
The  following  standard  strain  transformation  equation: 
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(2.72) 


E^=^+^«>s2P+^sin2P 

can  be  used  to  relate  the  measured  gage  strain  with  the  global 
axis  strains.  These  relationships,  for  the  gages  shown  in 
figure  2.20  are: 
for  gage  1  with  =  +45° 


_e,-+E,  Yx, 

2  2 


(2.73) 


for  gage  2  with  =  0° 


(2.74) 


and  finally  for  gage  3  with  =  -45° 


2  2’ 


(2.75) 


Solving  for  the  global  strains  E,,E,,  and  y^y  terms  of  the 
gage  strains,  yields  the  following  relationships: 
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(2.76a) 

^  ^gl  ^g2  ^g3 

(2.76Z>) 

y^=^„-^g3- 

(2.76c) 

Using  equation  2.61,  the  material  axis  shear  strain  is  related 
to  the  global  axis  strains  by: 

7,2  =  2  sin  e  cos  0(6,  -  e,)  +  (cos^  0  -  sin^  0)7^.  (2.77) 

This  can  be  simplified  using  trigonometric  double  angle 
relationships  to  the  equation: 


Yi2  =  (£,  ~  e*)  sin  20 + 7^  cos  20. 


(2.78) 


Upon  the  substitution  of  equations  2.76a-c  into  equation  2.78, 
the  shear  strain  in  the  material  axis  can  be  written  in  terms 
of  the  measured  gage  strains  as: 


Yi2  =  (e,/  -  2£,2  +  sin  20  +  (e,,  -  )  cos  20.  (2.79) 


Finally,  the  shear  compliance  can  be  computed  by  substituting 

the  stress  strain  relationship  712  =  566X12  from  equation  2.57  and 
the  stress  transformation  equation  T,2  =  -sin0cos0a,  from  equation 
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2.60  into  equation  2.79,  yielding: 


(£,;  -  2e,2  +  e-ts)  sin  20  +  (e,/  -  e,j)  cos  20 
-sin0cos0o. 


(2.80) 


When  a  45  degree  specimen,  0  =  +45°,  is  used  for  shear  char¬ 
acterization  the  shear  compliance  equation  reduces  to 


2(-e,/+2eg2-ej5) 


(2.81) 
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3  EXPERIMENTAL  METHODS 


This  chapter  describes  the  experimental  equipment  and 
methods  used  during  the  characterization  of  the  Radel  X/IM7 
composite  material.  The  accuracy  of  the  long-term  creep 
predictions,  with  simultaneous  physical  aging,  ultimately 
depends  on  the  quality  of  test  data  used  in  developing  the 
prediction  models.  The  term  "quality  data"  is  used  to  imply 
the  data  was  obtained  with  the  best  possible  accuracy  using 
the  available  experimental  equipment.  The  experimental  methods 
used  were  closely  scrutinized  to  insure  quality  data  was 
collected.  Included  in  this  chapter  are  descriptions  of  testing 
equipment,  calibration  techniques,  specimen  preparation, 
strain  measurement  techniques,  data  collection,  aging  creep 
test  procedures  and  data  reduction  procedures  .  Each  description 
is  presented  to  document  actual  experimental  methods  used 
during  this  study. 

3 . 1  Testing  Equipment 

Creep  testing  utilized  an  Applied  Test  System  (ATS)  series 
2240  lever  arm  tester  with  a  20:1  load  ratio  and  a  8000  lb 
(35.6  N)  maximum  load  capacity.  Mounted  on  the  frame  was  an 
ATS  3610  oven  used  to  provide  the  temperature  environment  for 
the  tests.  The  oven  temperature  was  controlled  using  an  ATS 
2931  oven  controller.  The  frame,  oven,  and  controller  are 
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shown  in  figure  3.1.  Friction  type  grips  were  used  to  apply 
load  to  the  specimen.  The  frictional  force  was  increased  by 
inserting  a  fine  steel  mesh  between  the  specimen  and  the  grip 
face.  Additionally,  possible  slippage  of  the  specimen  in  the 
grips  was  eliminated  by  placing  a  1/8  inch  (0.3175  cm)  steel 
pin  thru  the  grips  and  specimen.  The  grips  were  attached  to 
the  ATS  creep  frame  using  3/8  inch  (0.9525  cm)  steep  pins  in 
a  conventional  pin-thru-the-grip  arrangement.  An  illustration 
of  how  the  specimen  was  attached  to  the  load  frame  is  shown 
in  figure  3.2. 

3.2  Calibration  Techniques 
3.2.1  Oven  Temperature 

The  temperature  distribution  in  the  ATS  3610  oven  was 
checked  prior  to  creep  testing.  The  temperature  was  measured 
at  3  locations  equally  spaced  along  a  8  inch  (20.32  cm)  length 
located  in  line  with  the  load  path.  Iron-constantan  (ANSI 
symbol  "J")  thermocouples  were  used  to  measure  temperatures. 
Initial  measurements,  at  oven  temperatures  up  to  450°F  (232°C)  , 
indicated  =  9°F  (5°C)  maximum  temperature  variation,  from  the 
top  to  bottom  of  the  8  inch  (20.32  cm)  test  length.  Owens 
[51]  reported  a  similar  temperature  variation,  »  1°F  (0.5°C) 
per  inch  of  gage  length  at  400°F  (204°C)  ,  in  the  same  ATS  oven. 
This  initial  temperature  variation  along  the  specimen  seemed 
high.  Sullivan  [66]  described  a  similar  temperature  variation 


75 


problem  encountered  using  Satec  test  chambers.  He  found  that 
placing  an  aluminum  cylinder  around  the  specimen  and  running 
preheated  nitrogen  through  the  oven  reduced  the  temperature 
variation  along  a  4  inch  (10.16  cm)  length  from  about  18°F 
(10°C)  to  5°F  (3°C)  .  The  improved  variation  in  the  Satec  oven 
is  close  to  the  initial  variation  in  the  ATS  oven.  Further 
improvement  is  the  ATS  temperature  variation  was  accomplished 
by  using  a  simple  aluminum  air  duct  placed  in  the  bottom  the 
ATS  oven.  This  air  duct  directed  the  heated  air  to  the  front 
of  the  oven  and  improved  the  mixing  of  the  air  within  the  oven. 
A  repeated  heating  test  revealed  that  the  temperature  variation 
had  decreased  to  less  than  4°F  (2°C) ,  over  the  8  inch  (20.32 
cm)  length.  This  small  temperature  variation  was  deemed 
acceptable  for  the  proposed  creep  testing.  Additionally,  it 
is  important  to  note  that  the  temperature  variation  along  the 
1/2  inch  (1.27  cm)  strain  gage  length  was  negligible. 

Initially,  the  desired  oven  temperature  was  set  using  an 
ATS  series  3821  oven  temperature  controller.  This  controller 
utilized  a  rheostat  knob  and  an  analog  needle  gage  for  setting 
and  monitoring  the  oven  temperature.  Unfortunately,  the 
pointing  needle  in  the  gage  was  broken.  This  made  it  nearly 
impossible  to  repeat  the  same  temperature  setting.  This  problem 
was  alleviated  by  installing  an  ATS  series  2931  oven  controller. 
This  controller  utilized  digital  thumbwheels  to  set  the  desired 
oven  temperature  and  a  type  E  thermocouple  to  sense  the  oven 


76 


temperature.  The  controller  was  calibrated  by  comparing  the 
digital  temperature  set  on  the  controller  to  oven  temperature 
readings  from  a  Fisher  Scientific  type  14-983-lOE  precision 
mercury  thermometer  and  an  Omega  type  J  thermocouple  ther¬ 
mometer.  The  thermocouples  and  the  thermometer  were  placed 
within  a  1/2  inch  (1.27  cm)  of  the  location  of  the  strain  gages 
in  the  test  oven.  The  mercury  and  E  thermocouple  thermometer 
readings  were  within  1 . 8°F  (1°C)  over  the  calibrated  temperature 
range  between  338-446°F  (170-230°C)  .  The  controller  thumbwheel 
settings  required  to  obtain  a  desired  oven  temperature  were 
recorded  in  9°F  (5°C)  increments  over  the  oven  temperature 
operating  range.  During  testing,  the  mercury  thermometer  had 
to  be  removed  from  the  oven.  Long-term  temperature  monitoring 
was  accomplished  by  connecting  the  J  thermocouple  to  a  GRAPHTEC 
strip  chart  recorder,  model  DX  SR6335,  with  a  J-1  thermocouple 
module.  The  strip  chart  record  was  used  to  check  for  power 
interruptions  and  verify  that  the  oven  temperature  remained 
constant  during  long-term  tests.  The  series  2931  controller 
was  found  to  hold  the  oven  temperature  within  »  1  °F  (.5°C) 
after  stabilizing  at  the  set  temperature.  The  oven  temperature 
was  also  found  to  be  slightly  affected  by  the  ambient  room 
temperature.  This  affect  required  an  adjustment  of  the  oven 
controller  setting  by  up  to  ±2°F  (+i.'’'7)  to  account  for  seasonal 
temperature  variations  and  periods  wht  n  the  air  conditioner 
in  the  lab  was  not  working  properly. 
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3.2.2  Quanch  Rate  Detexaination 

Booij  and  Palmen  [6]  reported  that  the  quench  rate  used 
to  cool  a  specimen  from  above  T,  to  the  aging  temperature 
influences  the  linear  viscoelastic  properties  of  the  resulting 
glassy  material  for  a  very  long  time.  In  order  to  compare 
viscoelastic  properties  at  temperatures  below  Tg,  it  is 
essential  to  know  the  thermal  history  of  the  specimen  as  it 
passed  through  Tg  during  cooling.  Very  slow  cooling  of  the 
specimen  produces  a  smaller  specific  volume  than  a  fast  cooling. 
Thus,  the  initial  free  volume  state  of  the  specimen  also  varies 
with  the  quench  rate.  A  fast  quench  rate  produces  a  larger 
differential  between  the  trapped  free  volume  state  and  the 
equilibrium  free  volume  state  than  is  obtained  with  a  very 
slow  cooling  rate.  Thus,  physical  aging  effects  (i.e.  free 
volume  relaxation)  are  more  prevalent  in  a  rapidly  quenched 
material  containing  a  large  free  volume  differential  from  the 
equilibrium  state. 

The  method  used  to  quench  specimens  from  above  T,  to  the 
aging  temperature  varies  between  investigators.  Booij  and 
Palmen  [6]  quenched  their  samples  using  liquid  nitrogen.  Kong 
[30]  quenched  his  graphite/epoxy  laminate  specimens  in  an 
ice-water  bath.  A  procedure  used  by  Sullivan  [66]  utilized 
high  pressure  air  to  quench  the  specimen  and  was  adopted  for 
use  in  this  study.  This  method  provided  the  only  practical 
means  to  quench  the  material  when  strain  gages  and  wires  are 
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attached  to  the  specimen.  The  injection  of  high  pressure  air 
into  the  oven  had  two  effects.  First,  the  air  cooled  the 
specimen  providing  the  desired  high  quench  rate.  Secondly, 
the  air  cooled  the  test  oven  to  the  desired  aging  temperature. 
Both  of  these  effects  were  characterized  using  a  90  degree 
specimen  instrumented  with  a  J  thermocouple  embedded  in  the 
material.  The  instrumented  specimen,  mounted  in  the  test  oven, 
was  initially  stabilized  at  a  temperature  Tg  +  18°F  (10°C), 
441°F  (227°C)  .  The  quenching  procedure  was  started  by  first 
turning  off  the  oven  controller.  Next,  the  oven  door  was 
opened  and  the  high  pressure  air  was  injected  into  the  oven 
around  the  specimen.  After  a  set  time,  the  oven  door  was 
closed  and  the  specimen  temperature  monitored  to  determine  the 
stabilization  temperature.  This  procedure  was  repeated  with 
the  duration  of  the  air  quench  varied.  A  plot  of  the  stabi¬ 
lization  temperature  vs.  quench  time  duration  is  shown  in 
figure  3.3.  This  plot  was  used  to  determine  the  quench  time 
required  to  obtain  a  desired  aging  temperature.  The  specimen 
temperature  was  also  monitored  as  a  function  of  time  during 
the  air  quench  using  the  GRAPHTEC  strip  recorder.  The  specimen 
temperature  vs.  time  data  for  five  different  quench  cycles  was 
plotted  and  is  shown  in  figure  3.4.  A  linear  regression  fit 
to  the  test  data  yielded  a  quench  rate  of  9.5°F/sec  (5.3°C/sec)  . 
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Another  concern  resulting  from  the  quenching  procedure  was 
the  formation  of  microcracks  in  the  matrix  material.  This 
concern  was  investigated  by  polishing  the  edge  of  a  90  degree 
specimen  and  photoreplicating  it  before  and  after  all  heat 
treatments  of  the  specimen.  A  few  small  cracks  were  found  in 
matrix  rich  regions  of  the  specimen  during  the  strain  gage 
curing  procedure.  The  maximum  length  of  the  cracks  was  about 
4  times  the  fiber  diameter  of  0.197  mil  (5  microns).  The 
microcracks  usually  connected  two  fibers  through  the  matrix 
rich  region.  Subsequent  repeated  quenches  and  creep  tests  of 
the  specimen  revealed  no  further  increase  in  crack  formation 
or  noticeable  growth  in  the  initially  formed  cracks.  Quenching 
of  laminates  formed  with  various  ply  angle  orientations  will 
likely  increase  the  formation  of  microcracks  due  to  the  high 
thermal  expansion  mismatch  between  plies. 

3.2.3  Specimen  Stress 

The  uniaxial  stress  applied  to  the  specimen  is  calculated 
by  dividing  the  measured  applied  load  by  the  specimen  cross- 
sectional  area.  Accurate  measurement  of  the  specimen  geometry 
and  applied  load  are  required  to  determine  an  accurate  applied 
stress.  The  specimen  geometry  was  measured  using  a  Starrett 
No.  722  digital  micrometer.  The  applied  load  on  the  specimen 
was  more  difficult  to  accurately  obtain  and  required  calibration 
of  the  loading  mechanism  of  the  creep  frame. 
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The  calibration  of  the  loading  mechanism  was  initially 
attempted  using  a  simple  balancing  procedure.  The  first  step 
of  this  procedure  was  to  balance  the  loading  arm  in  a  horizontal 
position  without  any  input  load  applied  to  the  arm.  The  term 
"null  force"  is  used  to  describe  the  amount  of  weight  required 
to  balance  the  arm  with  zero  input  load  applied.  The  null 
force  is  the  load  on  the  specimen  due  to  the  initial  unbalance 
of  the  loading  mechanism.  Next,  a  small  amount  of  lead  weight 
was  measured  and  added  as  an  input  load.  Following  this  step, 
a  balancing  load,  approximately  20  times  the  input  load,  was 
added  to  the  load  train  side  of  the  lever  arm.  This  weight 
was  tweaked  until  the  lever  arm  was  balanced  in  the  horizontal 
position.  The  actual  load  ratio  was  then  calculated  by  dividing 
the  balancing  load  by  the  input  load.  Unfortunately,  the  load 
ratio  didn't  remain  constant  and  depended  upon  the  magnitude 
of  the  balancing  load.  This  indicated  that  the  frictional 
forces  induced  by  the  three  knife  edge  contact  points  on  the 
lever  arm  must  be  included  in  the  frame  calibration. 

The  actual  calibration  procedure  used  during  this  study 
consisted  of  three  steps.  The  first  step  was  to  calibrate  a 
500  lb  (2224  N)  capacity  Baldwin  SR-4  type  U-1  load  cell.  This 
calibration  step  was  accomplished  by  connecting  the  load  cell 
to  a  Vishay  2120  strain  gage  conditioner  and  recording  the 
output  voltage  as  known  weights  were  applied  to  the  cell.  The 
conditioner  used  an  excitation  voltage  of  1.994  volts  and  a 
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calibration  voltage  of  .583  volts  per  1000  units.  The  load 
cell  was  zeroed  while  loaded  by  its  own  weight.  This  insured 
that  only  the  applied  loads  to  the  cell  were  measured.  The 
output  voltage  and  applied  load  data  was  plotted  and  fit  using 
a  linear  regression  routine  and  is  shown  in  figure  3.5.  The 
resulting  equation  relating  the  force  in  the  load  cell  to  the 
output  voltage  is: 

/orcg(/6)  =  73.97618498*(Vo/rage)-0.05480718  .  (3.1) 

The  significant  digits  displayed  in  the  equation  are  carried 
only  to  insure  accurate  calculation  of  the  linear  regression 
fit  of  the  data. 

The  second  calibration  step  consisted  of  using  the  cali¬ 
brated  load  cell  to  determine  the  actual  force  applied  to  the 
specimen.  The  load  cell  was  carefully  installed  in  the  load 
train  path  where  the  specimen  would  be  located.  The  load  cell, 
loaded  by  its  own  weight,  was  zeroed  and  the  proper  excitation 
and  calibration  voltages  were  checked.  Next,  the  support  under 
the  lever  arm  was  removed  and  the  null  force  voltage  was 
measured.  Conversion  of  the  null  voltage  to  a  force,  using 
equation  3.1,  yielded  the  force  on  the  specimen  with  zero  input 
load  applied  to  the  lever  arm.  The  calibration  process  continued 
with  gradual  increases  in  the  amount  of  input  load  as  the 
corresponding  output  voltages  were  recorded.  A  Mettler  PJ6000 
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balance  was  used  to  measure  the  weight  of  the  input  load  to 
an  accuracy  of  ±0.1  grams.  Using  equation  3.1,  the  output 
voltage  data  was  converted  to  pounds  force  and  plotted  verses 
the  input  load  measured  in  grams.  This  plot  is  shown  in  figure 
3.6.  The  results  of  a  linear  regression  fit  to  the  data  yielded 
the  following  relationship  between  input  load  and  force  measured 
by  the  load  cell: 

LoadCellForceilb.)  =  Qm56\62\2*InputLoad{gm.)  +  6.59Ae2>966e  .  (3.2) 

The  6.594689666  constant  represents  the  null  force  with  the 
load  cell  attached. 

The  final  calibration  step  consisted  of  balancing  the  lever 
arm  in  a  way  that  the  specimen  would  be  unloaded  when  zero 
input  load  was  applied.  This  was  accomplished  by  replacing 
the  load  cell  with  the  top  mounting  grip,  specimen,  and  strain 
gage  wires.  Next,  the  lever  arm  was  balanced  to  a  level 
position  by  moving  a  weight  attached  to  the  lever  arm,  as 
illustrated  in  the  initial  balance  diagram  in  figure  3.2. 
Finally,  the  lower  grip  and  attachment  hardware  was  attached 
to  the  specimen.  The  weight  of  this  lower  assembly  rested  on 
the  lower  frame  plate  until  the  load  was  applied  as  shown  in 
figure  3.2.  Upon  load  application,  the  lower  assembly  travels 
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slightly  upward  until  the  lower  nut  contacts  the  lower  frame 
plate  to  counteract  to  applied  load.  The  resulting  calibration 
equation  after  this  final  balancing  step  is: 

SpecimenForce(lb.)  =  Q.QA56\61\2*InputLoad{gm.).  (3.3) 

An  estimate  of  the  error  in  the  specimen  load  was  made  by 
summing  the  standard  error  estimate  of  the  two  dependent  linear 
regression  fits  as  suggested  in  [53] .  This  yielded  an  applied 
specimen  force  with  +0.17  lb  (0.75  N)  error  in  the  0  to  150 
lb  (667  N)  calibration  range. 

3 . 3  Specimen  Preparation 

The  Radel  X/IM7  composite  material  used  in  this  study  was 
supplied  by  NASA  Langley.  A  grey  scale  imaging  program  was 
used  to  estimate  the  composite  fiber  volume  fraction.  The 
program  separately  counts  the  number  of  pixels  contained  in 
the  fibers  and  matrix  material  from  a  scanned  video  image  of 
a  micrograph  of  the  edge  of  a  90  degree  specimen.  Using  a 
weighted  average  from  three  separate  scans  of  the  micrograph 
the  fiber  volume  fraction  was  estimated  at  0.65.  The  material 
was  received  in  11  inch  (28  cm)  by  11  inch  (28  cm)  unidirectional 
panels  with  a  12  ply  thickness.  A  50  ton  (445  KN)  hot  press, 
programmed  to  follow  the  consolidation  cycle  shown  in  figure 
3.7,  was  used  to  make  the  panels.  Upon  receipt  at  VPI  &  SU 
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the  panels  were  C-scanned  to  locate  any  irregularities  in  the 
material.  Panel  regions  with  suspected  irregularities  were 
avoided  during  the  fabrication  of  specimens.  A  0  degree 
reference  axis  was  produced  on  the  panel  by  pealing  off  a  thin 
strip  of  material  along  an  edge  parallel  to  the  fibers.  A 
water-cooled  diamond  wheel  abrasive  saw  was  used  to  cut  the 
panels  into  8  inch  (20.32  cm)  long  by  1/2  inch  (1.27  cm)  wide 
off-axis  specimens.  Standard  machining  methods  were  used  to 
cut  the  specimens  at  0,  90,  45,  and  30  degree  angles  from  the 
reference  fiber  direction. 

Since  the  same  specimens  would  be  used  to  characterize 
physical  aging  at  various  temperatures,  a  means  of  rejuvenating 
the  specimen,  i.e.,  erasing  the  previous  physical  aging  history, 
was  required.  Booij  and  Palmen  [6]  recommended  raising  the 
material  to  a  temperature  of  Tg  +  18°F  (10°C)  to  eliminate  the 
effects  of  the  previous  thermal  history.  They  stated  that  the 
volume  equilibrium  state  is  reached  in  a  few  seconds  at  this 
temperature.  Other  investigators  have  used  rejuvenating 
temperatures  ranging  from  Tg  +  12.6°F  (7°C)  [70]  to  Tg  +  72°F 
(40°C)  [66]  .  The  important  point  is  that  the  temperature  amount 
above  Tg  is  sufficient  to  erase  the  previous  thermal  history 
during  the  time  period  the  material  is  held  at  the  temperature. 

The  Tg  of  the  Radel  X/IM7  composite  was  determined  using 
a  DuPont  2100  Differential  Scanning  Calorimetry  (DSC)  at  a 
10°C/min.  heating  rate.  Various  methods  of  defining  the  Tg 


85 


from  the  DSC  trace  exist  [141]  .  The  midpoint  of  the  glass 
transition  region,  utilizing  a  half  height  calculation  method, 
was  used  in  this  study.  Using  this  method,  T,  is  determined 
at  the  heat  flow  value  halfway  between  the  onset  and  the  end 
of  the  glass  transition  region.  The  DCS  trace  for  the  Radel 
X/IM7  material  is  shown  in  figure  3.8.  The  Tg  value  calculated 
was  424°F  (217. 8°C).  Osiroff  [50]  reported  a  similar  Tg  value 
of  419°F  (215°C)  for  a  Radel  X/T650-42  composite  material.  A 
temperature  of  437°F  (225°C)  was  chosen  as  a  rejuvenating 
temperature  for  the  Radel  X/IM7  material.  This  temperature 
corresponded  to  the  end  of  the  glass  transition  region  defined 
by  the  DSC  trace  shown  in  figure  3.8.  Preliminary  tests 
verified  the  material  was  rejuvenated  by  holding  at  this 
temperature  for  15  minutes. 

The  Radel  X/IM7  material  was  also  subjected  to  a  Thermal 
Gravimetric  Analyzer  (TGA)  test  to  determine  if  the  material 
would  degrade  at  the  maximum  testing  temperature.  As  shown 
in  figure  3.9,  no  significant  weight  loss  occurred  up  to 
temperatures  of  788°F  (420°C)  .  This  temperature  was  well  beyond 
the  437°F  (225°C)  temperature  used  for  rejuvenation  of  the 
specimens . 

3.4  Strain  Gage  Installation 

Creep  strain  measurement  was  accomplished  using  resistance 
foil  strain  gages.  Strain  gages  were  chosen  over  extensometers 
for  several  reasons.  Most  notable  were  the  cheaper  cost  and 
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availability  of  strain  gages  over  extensometers .  Another 
reason  was  the  versatility  strain  gages  offer  over  extensometers 
in  measuring  strain  in  directions  other  than  the  normal  axial 
and  transverse  directions.  This  versatility  was  particularly 
useful  for  measuring  strains  in  the  three  directions  required 
to  compute  the  shear  compliance  term  ■^66* 

The  strain  gages  used  were  from  the  WK  gage  series  made 
by  Micro-Measurements.  The  normal  operating  temperature  range 
of  the  WK  gages  is  between  -452°  to  550°F  (-269°  to  290°C)  . 
WK-06-125TM-350  general  purpose  2-element  90  degree  tee  rosette 
gages  were  used  on  the  0  and  90  degree  specimens.  WK-06- 
125RA-350  general  purpose  3-element  45  degree  rectangular 
rosette  gages  were  used  on  the  45  and  30  degree  specimens. 
The  350-Ohm  resistance  gages  were  chosen  over  120-Ohm  gages 
in  order  to  reduce  the  effect  of  localized  heating  of  the 
material  under  the  gages  et  al  [69].  The  K-alloy  gages  were 
fully  encapsulated  with  high-endurance  leadwires  pre-attached. 
The  leadwires  were  attached  to  bondable  terminals  placed  away 
from  the  gages.  This  isolated  any  local  damage  caused  by  the 
hot  soldering  tip  away  from  the  gage  location.  Micro- 
Measurements  430-FST  four  conductor  cable  was  used  to  extend 
the  leadwires  to  the  strain  gage  signal  conditioning  amplifier 
unit.  The  individual  conductor  wires  in  this  cable  were  twisted 
together  and  covered  with  a  braided  shield  and  Teflon  jacket 
to  minimize  line  noise.  The  length  of  each  cable  was  a  constant 
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six  feet  to  insure  each  gage  experienced  the  same  leadwire 
resistance.  The  cable  was  rated  for  use  up  to  a  temperature 
of  500°F  (260°C)  . 

Strain  gages  were  placed  on  both  sides  of  each  specimen. 
During  testing,  both  an  active  and  a  thermal  compensating 
"dummy"  specimen  wf  placed  in  the  oven.  The  active  and  dummy 
specimens  were  of  the  same  fiber  orientation  and  experienced 
the  same  thermal  history.  The  four  gages,  two  each  from  the 
active  and  dummy  specimens,  were  wired  in  a  full  bridge 
Wheatstone  circuit  [12] .  A  diagram  of  the  circuit  is  shown 
in  figure  3.10.  Specific  placement  of  the  active  and  dummy 
gages  in  the  circuit  allowed  electronic  averaging  of  bending 
effects  and  cancelation  of  thermal  effects.  Thus,  the  output 
signal  is  entirely  due  to  the  applied  tension  load. 

Actual  mounting  of  the  strain  gages  on  the  specimens  was 
done  following  the  procedures  recommended  by  Micro-Measurements 
[39]  .  Both  the  strain  gage  and  the  bondable  terminal  were 
glued  to  the  specimen  using  M-bond  610  high  temperature 
adhesive.  An  initial  problem  encountered  was  that  the  pre¬ 
attached  leadwires  would  bond  to  the  specimen.  This  provided 
an  electrical  conductance  path  that  shorted  out  the  strain 
gage.  This  problem  was  solved  by  adding  a  layer  of  mylar  JG 
high  temperature  tape  between  the  pre-attached  lead  wires  and 
the  specimen  surface  before  curing.  This  kept  the  lead  wires 
from  adhering  to  the  specimen.  The  adhesive  was  cured  at  347°F 
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(175°C)  for  2  hours  in  a  model  838F  Fisher  Scientific  Isotemp 
programmable  oven.  The  oven  temperature  was  increased  and 
decreased  at  a  rate  9°F  (5°C)  per  minute  as  illustrated  in  the 
strain  gage  cure  cycle  shown  in  figure  3.11.  Following  this 
cure  cycle,  the  specimens  were  removed  from  the  clamping 
fixtures  and  the  high  temperature  mylar  tape  was  removed. 
Next,  the  specimens  were  subjected  to  the  two  level  post  cure 
cycle  also  shown  in  figure  3.11.  The  first  level,  a  12  hour 
hold  at  212°F  (100°C),  was  required  as  a  drying  cycle  to  remove 
absorbed  moisture  from  the  specimens.  Initial  specimens  were 
found  to  form  surface  blisters  when  heated  directly  to  the 
upper  post  cure  temperature.  Investigators  at  NASA  Langley 
had  observed  the  same  phenomena  and  recommended  using  the 
drying  cycle  to  avoid  the  blister  formation  associated  with 
trapped  moisture.  The  second  level,  a  8  hour  hold  at  437°F 
(225°C)  ,  was  used  as  a  high  temperature  post  cure  cycle  to 
insure  the  M-bond  610  strain  gage  adhesive  was  fully  cured. 

After  the  post  cure  cycle,  the  specimens  were  placed  inside 
a  desiccator  until  they  were  actually  wired  and  used  in  the 
aging  tests.  Monitoring  the  weight  of  a  sample  specimen  during 
this  storage  period  revealed  no  significant  increase  or  decrease 
in  weight.  Desiccant  material  was  also  placed  inside  the  ATS 
oven  to  minimize  any  moisture  absorption  during  the  testing 
period. 
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The  final  step,  in  preparation  for  testing,  was  to  solder 
the  pre-attached  lead  wires  to  the  bondable  terminals  and  the 
extension  cables.  This  simple  sounding  task  turned  out  to  be 
physically  quite  difficult.  The  pre-attached  leadwires  were 
quite  fragile  and  easily  pulled  off  the  strain  gage.  In 
addition,  the  Micro  Measurements  570-28R-20  high  temperature 
solder  was  found  to  rapidly  oxidize  the  soldering  tip  making 
wetting  and  flowing  of  the  solder  difficult.  After  completion 
of  the  soldering  proce  >s,  the  resistance  of  each  gage  was 
checked  at  the  end  of  the  extension  cables.  This  provided  a 
good  spot  check  of  each  gage  to  insure  that  the  solder  joints 
were  satisfactory  completed.  Next,  the  gages  were  wired  in 
the  full  bridge  circuit  and  connected  to  the  Vishay  signal 
conditioning  amplifier  unit.  The  circuit  was  then  checked  for 
proper  electrical  balance.  After  passing  these  screening 
tests,  the  solder  joints  were  cleaned  with  M-Line  rosin  solvent 
to  remove  any  residual  flux  material.  Following  cleaning,  the 
specimens  were  ready  for  mounting  in  the  test  oven. 

The  8  inch  (20.31  cm)  long  active  specimen  was  mounted 
with  1.5  inches  (3.81  cm)  inserted  into  each  grip.  This 
resulted  in  a  test  length  of  5  inches  (12.7  cm)  and  a  10  to 
1  length  to  width  aspect  ratio.  The  grips  were  placed  in  a 
special  fixture  that  eliminated  any  load  transfer  to  the 
specimen  while  the  Allen  screws  were  tightened  on  each  grip. 
After  removal  from  the  fixture,  a  6  inch  (15.2  cm)  metal  ruler 
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was  placed  between  the  grips  to  avoid  any  bending  of  the 
specimen  during  transfer  to  the  test  frame.  The  dummy  specimen, 
used  for  thermal  compensation,  was  initially  supported  from 
the  upper  3/8  inch  (0.9525  cm)  pin  in  the  load  train.  Later, 
a  method  of  mounting  the  specimen  from  the  lower  grip  was 
devised  in  order  to  eliminate  any  possible  loading  of  the  dummy 
specimen  by  the  weight  of  the  lead  wire  extension  cables. 

Once  mounted  in  the  test  frame,  each  strain  gage  circuit 
was  checked  for  excessive  zero  drift  (i.e.  drift  in  the  bridge 
output  voltage  with  no  load  applied) .  Possible  explanations 
for  zero  drift  include  a  gradual  change  in  the  gage  zero 
resistance,  oxidation  of  the  solder  joints,  temperature 
changes,  moisture  changes,  or  a  slight  load  applied  to  the 
specimen.  The  initial  zero  drift  for  the  three  gages  on  the 
45  degree  specimen  at  392°F  (200°C)  is  shown  in  figure  3.12. 
The  zero  drift  during  this  14  hour  test  seemed  excessive  for 
all  three  gages.  The  oven  was  therefore  cooled  in  order  to 
check  the  solder  joints  on  the  strain  gages.  The  pre-attached 
leadwires  seemed  loose  in  several  joints,  possibility  due  to 
oxidation  of  the  solder,  and  were  resoldered.  Another  zero 
drift  check  at  392°F  (200°C)  revealed  a  dramatic  decrease  in 
voltage  drift  as  shown  in  figure  3.13.  At  this  point,  it  was 
decided  to  check  the  effect  of  mounting  the  dummy  specimen 
from  the  bottom  grip  instead  to  the  top  grip.  Mounting  from 
the  bottom  grip  would  eliminate  any  loading  from  the  lead  wire 
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extension  cables  that  were  attached  below  the  strain  gages. 
The  drift  measurements,  after  this  remounting,  were  observed 
to  change  direction  for  gages  1  and  2,  as  shown  in  figure  3.14. 
It  was  impossible  to  determine  if  the  changes  were  due  to 
unloading  of  the  dummy  specimen  or  due  to  a  slight  change  in 
the  balance  of  the  loading  system  that  would  place  a  small 
load  on  the  active  specimen.  The  zero  drift  rate  was  observed 
to  always  decrease  with  time.  Micro  Measurements  [40]  reports 
the  same  decrease  in  drift  rate  with  time  for  most  gage  alloys. 

3 . 5  Data  Collection 

Initial  tests  using  a  Vishay  2100  series  signal  conditioner 
amplifier  were  plagued  with  excessive  line  noise  in  the  strain 
gage  output  voltage.  The  source  of  the  line  noise  was  traced 
to  the  circulating  fan  within  the  ATS  oven.  An  oscilloscope 
was  used  to  determine  that  the  noise  was  concentrated  at  60 
Hz  with  a  peak  magnitude  of  about  20  mv.  This  voltage  equated 
to  approximately  ±15  micro  strain  change  in  the  strain  readings 
on  a  90°  test  specimen.  The  line  noise  was  eliminated  by 
switching  to  a  Vishay  2210  series  signal  conditioner  amplifier. 
This  unit  incorporated  a  noise  filtering  system  that  was  set 
to  eliminate  line  noise  above  lOHz.  A  comparison  of  the  strain 
measurements  obtained  from  the  two  amplifiers  is  shown  in 
figure  3.15.  The  Vishay  2210  amplifier  clearly  reduced  the 
line  noise  and  was  used  exclusively  thereafter  in  this  study. 
Additionally,  all  test  equipment,  except  the  test  oven,  were 
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plugged  into  a  Trip  Lite  BC-750a  battery  backup  system  to  guard 
against  short  power  interruptions.  This  unit  also  protected 
the  equipment  against  voltage  spikes  and  suppressed  line  noise 
from  the  power  distribution  system. 

All  strain  gage  measurements  were  taken  with  the  amplifier 
gage  factor  set  at  2.0.  Adjustments  for  the  actual  gage  factor 
and  temperature  induced  variations  were  included  during  the 
data  reduction  process  using  the  procedures  recommended  by 
Micro  Measurements  [42]  .  The  strain  gages  were  calibrated 
using  a  shunt  calibration  procedure  [12,44].  The  shunt  was 
applied  at  the  connector  box  where  the  extension  cable  was 
attached  to  the  conditioner  unit.  The  equivalent  strain,  due 
to  the  applied  shunt,  was  calculated  using 


where 


P  = _ -1 _ 

nFiR^+R,) 

^UNT  -  equivalent  shunt  Strain 
Rg  =  resistance  of  gage 
R,  =  resistance  of  shunt 
n  =  number  of  active  gages 
F  =  assumed  gage  factor. 


(3.4) 


Precision  87,150  ohm  and  174,650  ohm  calibration  resistors 
were  used  to  simulate  1000  and  500  micro  strain  respectively 
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in  the  full  bridge  circuit.  The  bridge  excitation  voltage 
was  determined  using  procedures  recommended  by  Micro  Mea¬ 
surements  [41] .  The  composite  material  was  characterized  as 
a  poor  heat-sink  material.  The  recommended  power  density  for 
this  type  material  is  about  0.086  watts/in^  (0.138  KW/m^)  for 
each  gage.  The  larger  grid  area  WK-06-125TM-350  gages,  used 
on  the  0  and  90  degree  specimens,  were  excited  with  1.50  volts. 
The  smaller  grid  area  WK-06-125RA-350  gages,  used  on  the  45 
and  30  degree  specimens,  were  excited  with  0.96  volts.  The 
estimated  gain  required  to  obtain  the  desired  bridge  output 
voltage  with  the  shunt  applied  was  calculated  using: 

4£o 

where 

G  =  estimated  gain 
Eg  =  bridge  outputvoltage 
n  =  numberof  active  gages 
F  =  assumed  gage  factor 
E^uffr  =  equivalent  shunt  Strain 
V  =  bridge  excitationvoltage . 

The  estimated  gain  was  set  on  the  signal  conditioning  amplifier 
and  delicately  adjusted  in  order  to  obtain  the  proper  cali¬ 
bration  output  voltage  for  each  particular  gage. 
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The  actual  output  voltages  from  the  strain  gages  were 
recorded  using  an  IBM  PC  computer  with  a  Data  Translation 
2805/5716  digitial/analog  (D/A)  board.  A  BASIC  program  written 
for  data  acquisition  named  API  [46]  was  modified  for  use  during 
this  study.  One  modification  fixed  a  timing  problem  that 
occurred  at  midnight  when  the  day  changed.  Another  modification 
incorporated  a  procedure  to  take  three  readings  from  the  D/A 
board  and  average  them  to  obtain  the  reading  at  the  particular 
time  step.  The  averaging  process  was  installed  as  a  possible 
method  to  minimize  the  line  noise  while  using  Vishay  2100 
signal  conditioning  amplifier.  Unfortunately,  the  averaging 
process  had  little  effect  on  the  line  noise  concentrated  at 
60  Hz  but  was  left  in  the  program  to  smooth  the  data.  The 
time  required  to  take  three  data  samples  for  the  averaging 
process  was  negligible  since  the  DT  board  has  an  A/D  throughput 
of  2500  samples  per  second.  The  resolution  of  the  2805/5716 
D/A  board  was  16  bits.  This  allowed  the  +10  to  -10  voltage 
range  to  be  divided  with  a  resolution  of  0.0003  volts.  This 
voltage  resolution  equated  to  strain  resolution  of  less  than 
0.2  HE  for  the  strain  per  volt  calibration  ranges  used  in  this 
study. 

Ideally,  strain  data  should  be  taken  on  a  logarithmic  time 
spacing  for  creep  testing.  This  feature  could  not  be  easily 
added  to  the  API  program.  Instead,  the  data  was  taken  at 
increasing  time  intervals  and  later  reduced  to  approximate 
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logarithmic  spacing  using  a  post-processing  program.  The  data 
was  collected  at  time  intervals  of  1,  10,  100,  and  1000  seconds. 
The  number  of  points  collected  at  each  time  interval  were  350, 
100,  100  and  up  to  1450  respectively.  This  data  collection 
scheme  allowed  choosing  the  final  data  set  at  time  intervals 
approximately  equal  to  0.1  log  time  units  for  a  period  up  to 
16.9  days. 

3 . 6  Aging  Creep  Test  Procedures 

Both  short-term  and  long-term  creep  tests  were  conducted 
during  this  study.  The  results  of  the  short-term  (momentary) 
creep  tests  were  used  to  characterize  the  aging  rate  of  the 
material  as  described  in  section  2.2.5. 1.  The  results  of  the 
long-term  creep  tests  were  used  to  verify  the  physical  aging 
models.  The  procedures  used  to  conduct  both  test  types  were 
identical,  except  for  the  duration  of  the  test. 

The  testing  sequence  was  started  by  heating  the  specimens 
to  the  rejuvenation  temperature  of  437°F  (225°C)  and  holding 
for  a  minimum  of  15  minutes.  Initial  tests  revealed  that  the 
90  degree  orientations  were  particularity  susceptible  to 
warping  when  heated  above  Tg.  The  warpage  was  minimized  by 
placing  a  six  inch  metal  ruler  between  the  grips  before  heating 
above  Tg.  The  ruler  provided  additional  rigidity  to  the  specimen 
and  solved  the  warpage  problem.  As  a  precaution,  the  ruler 
was  used  on  all  other  angle  ply  orientations  to  minimize  the 
chance  of  warping  occurring.  Before  the  quench,  the  oven 
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controller  was  set  to  90°F  to  keep  the  oven  from  heating  during 
the  cooling  process .  After  the  desired  air  quench  was  conducted, 
the  metal  ruler  was  removed  and  the  oven  door  was  closed.  The 
oven  temperature  typically  overshot  the  desired  aging  tem¬ 
perature  by  about  5°F  (3°C)  while  the  oven  temperature  sta¬ 
bilized.  The  aging  start  time  and  date  was  consistently  started 
when  the  desired  aging  temperature  was  reached  after  this 
slight  temperature  overshoot.  Upon  reaching  the  desired  aging 
temperature,  the  oven  controller  was  reset  to  the  aging  tem¬ 
perature.  Next,  a  sequence  of  short-term  tests  were  conducted 
at  aging  times,  t,,  of  1,  3,  9,  27,  and  81  hours  after  the 
quench.  These  particular  aging  times  were  selected  to  allow 
sufficient  recovery  time  between  tests  and  to  provide  a  man¬ 
ageable  testing  schedule.  Before  starting  the  first  creep 
test  after  the  quench,  the  amp  balance,  excitation  voltage, 
circuit  balance,  and  shunt  calibration  voltage  for  each 
Wheatstone  circuit  was  checked  and  properly  adjusted. 

The  API  data  acquisition  program  was  initiated  for  each 
creep  and  recovery  test  in  order  to  create  a  separate  data 
file  for  each  short-term  test  segment.  The  first  entry  in  the 
data  file  was  the  time  and  date  the  local  test  segment  started. 
This  information  was  used  to  relate  the  local  time  data, 
referenced  to  the  start  of  the  local  test  segment,  to  the 
global  time  scale,  referenced  to  the  aging  start  time  and  date 
for  the  entire  isothermal  test  sequence. 
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3.7  Data  Reduction  Procedures 

The  data  reduction  process  was  initiated  by  calculating 
the  aging  time,  t,,  for  each  creep  test  segment.  This  was 
accomplished  by  subtracting  the  global  aging  start  time  and 
date  from  the  local  segment  starting  time  and  date  stored  in 
each  data  file.  Next,  several  adjustments  were  made  to  the 
strain  gage  data.  First,  the  strain  data  for  each  creep  test 
was  adjusted  by  subtracting  the  extrapolated  residual  recovery 
voltages  from  the  previous  recovery  period  as  explained  in 
section  2.2.5. 1.  The  extrapolated  recovery  voltages  for  each 
gage  were  calculated  by  fitting  the  last  50%  of  the  recovery 
data  with  a  linear  equation.  An  example  of  the  linear  fit, 
i.e.,  the  recovery  voltages  for  the  45° specimen  at  365°F  (185°C) , 
after  the  t',=9  hr  creep  test,  is  shown  in  figure  3.16.  The 
global  time  scale  was  used  in  the  recovery  data  so  that  the 
extrapolated  recovery  voltages  could  be  linked  to  the  time 
scale  in  next  creep  test.  Also,  the  total  extrapolation  of 
the  recovery  data  into  the  next  creep  test  period  was  no  more 
than  16%  of  the  total  recovery  period.  The  subtraction  of  the 
recovery  voltages  also  eliminated  any  accumulated  drift  in  the 
strain  gages  by  effectively  zeroing  the  voltages  before  each 
creep  test  segment. 

The  next  adjustment,  applied  to  each  strain  gage  reading, 
corrected  the  2.0  gage  factor  used  during  data  collection  to 
the  quoted  manufacturers  gage  factor  adjusted  to  the  actual 
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test  temperature  using  the  procedures  recommended  by  Micro 
Measurements  [42]  .  Finally,  the  strain  gage  data  was  corrected 
for  transverse  sensitivity  effects  using  the  procedures  rec¬ 
ommended  by  Micro  Measurements  [43]  .  The  transverse  sensitivity 
correction  was  not  required  for  the  90  degree  specimen  since 
the  ratio  of  transverse  to  axial  strain  was  very  small.  The 
predicted  error  in  the  strain  measurements  without  the 
transverse  sensitivity  was  less  than  1%  for  the  S22  direction 
and  approximately  150%  for  the  ^21  direction.  Since  the 
fiber/transverse  coupling  compliance  /  ■S'2if  is  highly  susceptible 
to  transverse  sensitivity  error  in  the  90  degree  specimen,  the 
preferred  method  to  determine  the  coupling  compliance  is  to 
use  a  0  degree  specimen.  The  0  degree  specimen  produces  larger 
strains  in  the  transverse  direction  and  permits  a  more  accurate 
determination  of  the  fiber/transverse  coupling  compliance,  5,2 
term. 

Ideally,  the  load  is  applied  instantaneously  as  a  step 
function  at  the  start  of  a  creep  test.  In  practice,  the  load 
is  actually  applied  over  a  certain  period  of  time.  During 
this  study,  the  load  was  applied  by  lowering  the  table  supporting 
the  weight  by  manually  turning  the  loading  handle  counter¬ 
clockwise.  Typically,  the  loading  process  was  accomplished 
in  about  4  seconds.  The  creep  response  resulting  from  the 
ramp  loading  will  obviously  be  different  from  that  of  a  step 
loading.  The  difference  in  creep  responses  will  become  neg- 
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ligible  after  a  certain  time,  usually  referred  to  as  the  initial 
transient  period.  Lou  and  Schapery  [36]  recommended  waiting 
a  period  equal  to  five  times  the  initial  loading  period  to 
allow  the  compliance  response  from  a  ramp  load  to  converge  to 
the  step  loading  response.  During  this  study,  the  calculation 
of  the  compliance  response  was  started  20  seconds  after  the 
load  was  applied.  The  sensitivity  of  the  aging  rate  to  the 
20  second  transient  period  was  checked  by  recalculating  the 
aging  rates  using  a  60  second  transient  period.  The  variation 
in  aging  rate  between  the  20  and  60  second  transient  periods 
ranged  from  -0.3%  to  +3.3%  for  the  five  test  temperatures  used 
with  the  90  degree  specimen.  Since  the  variation  was  small, 
the  start  of  the  creep  compliance  data  was  computed  after  the 
20  second  transient  period. 
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4  CHARACTERIZATIOK  RESULTS  AND  DISCUSSION 


This  chapter  contains  the  results  and  discussion  of  the 
momentary  creep  tests  conducted  to  characterize  the  effect  of 
physical  aging  on  the  Radel  X/IM7  composite.  The  physical 
aging  theory  developed  by  Struik  [59]  was  successfully  applied 
to  the  orthotropic  composite  material.  The  compliance  response 
of  the  composite  was  characterized  in  each  of  the  principal 
material  directions. 

4 . 1  Room  Temperature  Material  Response 

Room  temperature  tensile  tests  were  conducted  to  determine 
the  basic  material  properties  of  the  "as  received"  composite 
material.  The  results  were  used  to  estimate  the  linear  elastic 
range  for  the  90  and  0  degree  orientations.  The  tests  were 
conducted  using  an  Instron  4204  test  machine  with  a  MTS  model 
632.13B-20  extensometer .  The  specimens  were  held  using  wedge 
action  friction  grips.  The  active  specimen  length  was  5  in. 
(12.7  cm)  with  the  strain  measured  over  a  0.5  in.  (1.27  cm) 
gage  length.  The  ultimate  strength  and  stress-strain  behavior 
were  recorded  while  the  specimen  was  loaded  at  a  high  rate. 

The  resulting  stress-strain  curve  for  the  90  degree  specimen 
is  shown  in  figure  4.1.  A  linear  regression  routine  was  used 
to  calculate  a  modulus  value  of  1.07  Msi  (7.38  GPa) .  The 
stress  in  the  gage  section  at  failure  was  4.04  Ksi  (27.8  MPa) . 
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The  specimen  failed  near  the  grip  where  the  thickness  was 
slightly  reduced,  yielding  an  ultimate  failure  stress  of  4.46 
Ksi  (30.8  MPa)  .  A  comparison  of  the  test  results  to  commercially 
available  data  is  given  in  table  4.1.  The  commercial  data 
in  table  4.1  was  determined  for  a  composite  with  T650-42  fiber 
instead  of  IM7.  The  two  graphite  fiber  materials  have  similar 
material  properties.  Table  4.2  provides  a  list  of  the  properties 
of  T650-42  and  IM7  fiber  materials  for  comparison.  The  tensile 
modulus  obtained  during  the  room  temperature  testing  was  within 
8.5%  of  the  published  value,  but  the  tensile  strength  was  about 
35%  low.  The  low  tensile  strength  may  be  attributed  to  a 
slight  misalignment  of  the  specimen  or  the  absence  of  specimen 
end  tabs  which  help  lower  the  load  concentration  effect  near 
the  grips.  After  analyzing  the  results  of  Gramoll  et  al  [25], 
a  stress  level  equal  to  approximately  10%  of  the  ultimate  room 
temperature  strength  was  chosen  as  a  starting  stress  level  for 
linear  viscoelastic  testing.  Verification  tests  were  used  to 
show  that  the  material  remained  in  the  linear  viscoelastic 
range  at  this  stress  level. 

The  0  degree  stress-strain  curve  is  shown  in  figure  4.2. 
A  linear  response  was  obtained  up  to  a  stress  level  of  about 
185  Ksi  (1275  MPa) .  At  this  stress  level,  the  extensometer 
began  slipping  and  noticeable  cracking  of  the  fibers  was  heard. 
The  curve  was  linearly  extrapolated,  dashed  line  in  figure 
4.2,  to  the  ultimate  failure  stress  of  261  Ksi  (1800  MPa). 
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Again,  a  linear  regression  routine  was  used  to  determine  a 
fiber  direction  tensile  modulus  of  23.1  Msi  (159.3  GPa) .  The 
tensile  modulus  was  within  2.5%  of  the  published  value,  but 
similar  to  the  90  degree  specimen,  the  ultimate  tensile  strength 
in  the  0  degree  specimen  was  low  by  about  21%.  The  stress 
level  used  for  creep  testing  of  the  0  degree  specimen  was 
selected  at  a  value  equal  to  20%  of  the  failure  stress  in  order 
to  produce  measurable  strains  in  the  transverse  direction. 

Room  temperature  testing  also  included  ultimate  tensile 
tests  of  both  45  and  30  degree  specimen  orientations.  The 
specimens  were  held  using  the  friction  grips  that  were  utilized 
during  creep  testing.  This  insured  that  the  same  end  rotation 
constraint  was  used  on  the  unbalanced  off-angle  laminates  for 
both  ultimate  tension  and  creep  testing.  The  ultimate  specimen 
stress  of  the  45  and  30  degree  orientations  were  6.32  Ksi  (43.6 
MPa)  and  13.2  Ksi  (90.9  MPa)  respectively.  Both  specimens 
failed  in  shear  parallel  to  the  fibers  with  the  fracture 
originating  at  the  grip. 

4.2  90  Degree  Specimen  Tests,  ^22  Determination 

The  90  degree  specimen  orientation  was  used  to  determine 
the  transverse  direction  compliance,  S22,  in  the  linear  vis¬ 
coelastic  range  using  the  methods  explained  in  section  2.3. 
The  90  degree  specimen  orientation  was  tested  first  to 
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demonstrate  that  the  composite  physically  ages.  This  orien¬ 
tation  is  resin  dominated  and  was  expected  to  age  similarly 
to  a  neat  polymer. 

4.2.1  Linearity  and  Sequencing  Effects 

Initial  tests  of  the  90  degree  specimen  were  used  to  verify 
that  the  proposed  400  psi  (2.76  MPa)  applied  stress  level  would 
produce  a  linear  viscoelastic  response.  A  series  of  creep 
tests  were  conducted  at  different  stress  levels,  at  an  aging 
time  of  3  hours,  with  the  temperature  held  constant  at  392°F 
(200°C)  .  Between  each  test,  the  specimen  was  heated  above  T, 
and  air  quenched  back  to  the  aging  temperature.  The  compliance 
response  of  the  material  at  the  four  stress  levels  tested  is 
shown  in  figure  4.3.  The  compliance  curves  have  the  same  shape 
with  a  slight  offset  in  the  initial  compliance  values.  The 
slightly  increased  curvature  of  the  490  psi  (3.38  MPa)  curve 
was  attributed  to  a  slight  increase,  =®1-2°F  (1°C)  ,  in  temperature 
during  this  test.  The  consistency  of  the  compliance  curves 
over  the  stress  range  tested  indicated  that  the  material 
remained  in  the  linear  viscoelastic  range  at  the  test  tem¬ 
perature.  The  400  psi  (2.76  MPa)  stress  level  was,  therefore, 
used  for  all  the  aging  tests  of  the  90  degree  specimens. 

The  validity  of  the  creep  test  sequencing  procedure  was 
checked  using  the  procedure  described  in  section  2. 2. 5.1.  Both 
compliance  curves  shown  in  figure  4.4  were  obtained  at  338°F 
(170°C)  with  an  aging  time  of  9  hours.  The  compliance  curve 


104 


with  asterisk  symbols  was  the  third  creep  test  in  an  aging 
sequence  test.  The  compliance  curve  with  circle  symbols  was 
from  an  unsequenced  test  where  the  initial  two  creep  tests  of 
the  sequence  were  omitted-  The  compliance  curves  have  the 
same  shape  with  a  1.2%  offset  in  the  initial  compliance  mag¬ 
nitude.  This  slight  offset  can  be  attributed  to  the  error 
introduced  during  the  fitting  and  subtraction  of  the  recovery 
response  from  the  previous  creep  test  in  the  aging  sequence. 
The  overall  agreement  between  the  two  compliance  curves  dem¬ 
onstrates  the  applicability  of  using  the  aging  test  sequence 
to  characterize  the  aging  effect  in  the  linear  stress  range. 

4.2.2  Aging  Shift  Rate  Determination 

The  effect  of  physical  aging  on  the  transverse  direction 
compliance,  S22,  was  determined  using  a  sequence  of  short-term 
momentary  creep  tests.  The  tests  were  conducted  at  five 
temperature  levels  using  the  procedures  described  in  section 
3.6.  The  compliance  curves  resulting  from  the  short-term  tests 
are  shown  in  figures  4. 5-4. 9.  The  symbols  in  these  figures 
depict  the  actual  test  data  with  the  solid  line  representing 
the  numerical  fit  of  Struik's  equation  (eqn.  2.19)  to  the  data. 
An  average  m  value  of  0.4167  was  used  in  the  fitting  program. 
This  value  was  calculated  by  averaging  the  individual  m  values 
from  the  338,  365,  and  195°F  {170,  185  and  195°C)  aging  sequences . 
The  m  value  from  the  individual  creep  curves  showed  considerable 
scatter  varying  from  a  low  of  0.2921  to  a  high  of  0.4894. 
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The  scatter  is  most  likely  attributed  to  the  variability  of 
the  shape  of  the  compliance  curve  during  the  short  time  span 
of  the  momentary  creep  tests. 

The  momentary  creep  compliance  curves  were  shifted  to  form 
master  momentary  compliance  curves.  The  master  curve  data, 
referenced  to  the  81  hour  aging  time,  is  shown  in  figure  4.10. 
The  shifting  process  was  accomplished  using  the  procedures 
described  in  section  2. 2. 5. 2.  The  amount  of  horizontal  shift 
— loga  and  the  vertical  shift  iogB  for  each  creep  curve  are  given 
in  table  4.3. 

The  aging  shift  rate,  |i,  was  determined  for  each  aging 
test  sequence  using  a  linear  regression  routine.  A  plot  of 
the  aging  data,  -loga  versus  logt,,  and  the  calculated  aging 
shift  rate,  slope  of  the  curve,  for  each  test  sequence  are 
shown  in  figure  4.11.  The  aging  shift  rate  increased  as  the 
temperature  decreased,  until  a  temperature  of  365°F  (185°C) 
was  reached.  At  this  point,  the  aging  shift  rate  started  to 
decrease  as  the  temperature  decreased.  The  effect  of  tem¬ 
perature  of  the  aging  shift  rate  is  discussed  further  in  section 
4.5. 

4.2.3  TTSP  Master  Momentary  Compliance  Curve 

TTSP  master  momentary  compliance  curves  were  formed  by 
shifting  momentary  compliance  curves,  obtained  at  various 
temperatures  with  the  same  aging  time,  using  the  time- 
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temperature  superposition  principle.  The  90  degree  long-term 
test  was  started  at  an  aging  time  of  9  hours  with  the  temperature 
at  338°F  (170°C)  .  Several  methods  were  examined  for  obtaining 
the  TTSP  master  momentary  compliance  curve  at  this  reference 
state. 

The  first  method  numerically  fit  the  test  data  from  various 
temperatures  at  an  aging  time  of  9  hours  to  Struik's  equation 
(eqn.  2.19)  using  the  average  m  value  of  0 . 4167 .  The  horizontal 
shift  factors  were  calculated  using  equation  2.47.  The 
resulting  transverse  direction  master  momentary  compliance 
curve,  referenced  to  338°F  (170°C)  and  r,  =  9  hours,  is  shown 
in  figure  4.12. 

The  second  method  utilized  a  graphical  procedure,  described 
in  section  2. 2. 5. 3,  to  determine  the  shift  factors  required 
to  form  the  master  curve.  The  resulting  curve,  again  referenced 
to  338°F  (170°C)  and  /,  =  9  hours,  is  shown  in  figure  4.13.  An 
examination  of  the  shift  factor  summaries,  given  in  figures 
4.12  and  4.13,  shows  that  first  and  second  methods  give 
essentiality  the  same  shift  factors. 

The  final  shift  method  utilized  the  graphical  procedure 
to  obtain  the  master  curve  referenced  to  338°F  (170°C)  at  an 
aging  time  of  81  hours  instead  of  9  hours.  The  resulting 
curve,  referenced  to  338°F  (170°C)  and  t,  =  81  hours,  is  shown 
in  figure  4.14.  The  aging  shift  rate,  |i,  at  a  temperature  of 
338°F  (170°C)  ,  was  0.77.  Using  this  aging  rate  and  equation 
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2.48,  the  shift  factor  required  to  shift  the  =  81  hours  curve 
to  r,  =  9  hours  was  determined.  An  additional  small  vertical 
shift  was  also  included  during  shifting  procedure.  The 
resulting  TTSP  master  momentary  curve  is  shown  along  with  the 
curves  resulting  from  the  first  two  methods  in  figure  4.15. 
The  curves  show  essentially  the  same  shape  with  a  small  variation 
between  10^  to  10*  seconds.  The  effect  of  the  small  variation 
in  the  compliance  curves  on  long-term  predictions  is  discussed 
in  section  5.0. 

4.3  45  Degree  Specimen  Tests,  5^  Determination 

The  45  degree  specimen  orientation  was  used  to  determine 
the  shear  direction  compliance,  in  the  linear  viscoelastic 
range  using  the  methods  explained  in  section  2.3.  Exper¬ 
imentally  obtaining  the  shear  direction  compliance  data  pro¬ 
vided  the  unique  challenge  of  getting  twelve  strain  gages 
working  at  the  same  time  in  a  high  temperature  environment. 

4.3.1  Linearity  and  Sequencing  Effects 

The  45  degree  specimens  were  tested  at  an  applied  stress 
level  of  800  psi  (5.52  MPa).  Using  equation  2.60,  the  applied 
stress  was  found  to  produce  a  normal  stress  component  of  400 
psi  (2.76  MPa)  in  the  matrix  material.  This  was  the  same 
matrix  normal  stress  used  during  testing  of  the  90  degree 
specimens.  A  series  of  creep  tests  were  conducted  to  determine 
if  the  material  remained  in  the  linear  viscoelastic  range  at 
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the  applied  stress  level.  The  compliance  curves  from  tests 
at  different  stress  levels,  at  an  aging  time  of  3  hours  with 
the  temperature  held  constant  at  365°F  (185°C) ,  are  shown  in 
figure  4.16.  The  excellent  overlap  of  the  compliance  curves 
shows  that  the  material  remained  in  the  linear  viscoelastic 
range  up  to  the  800  psi  (5.52  MPa)  stress  level.  The  800  psi 
(5.52  MPa)  test  was  repeated  and  with  the  results  also  shown 
in  figure  4.16  as  triangular  symbols.  The  excellent  overlap 
with  the  previous  800  psi  (5.52  MPa)  test  demonstrates  the 
repeatability  of  the  creep  test  results. 

The  effect  of  sequencing  creep  tests  and  using  Boltzmann's 
superposition  principle  to  subtract  the  recovery  strain  from 
previous  creep  tests  was  also  checked  for  the  45  degree  specimen 
orientation.  Compliance  curves  from  tests  conducted  with  and 
without  sequencing  at  two  different  temperatures  and  aging 
times  are  shown  in  figure  4.17.  The  sequenced  and  unsequenced 
compliance  curves,  from  both  test  conditions,  have  the  same 
shape  with  a  1.8%  maximum  offset  between  them.  The  agreement 
between  sequenced  and  unsecpaenced  test  results  verifies  the 
applicability  of  using  the  aging  test  sequencing  procedure. 

4.3.2  Aging  Shift  Rate  Determination 

The  effect  of  physical  aging  on  the  shear  direction  com¬ 
pliance,  Sjft,  was  determined  at  three  temperature  levels.  The 
upper  two  temperature  levels,  399  and  392°F  (204  and  200°C)  , 
that  were  used  with  the  90  degree  specimen  orientation,  were 
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not  tested.  They  were  eliminated  in  view  of  time  constraints 
and  the  desire  to  reduce  the  risk  of  failing  the  specimen. 
The  compliance  curves  resulting  from  the  aging  sequence  tests 
are  shown  in  figures  4.18-4.20.  Again,  the  symbols  depict  the 
test  data  with  the  solid  line  representing  the  numerical  fit 
of  Struik's  equation  to  the  data.  The  individual  m  values 
from  the  three  aging  sequences  were  averaged  resulting  in  an 
average  m  value  of  0.4564.  This  value  was  used  in  the  computer 
fitting  program.  The  m  values  from  the  individual  creep  curves 
varied  from  a  low  of  0.3518  to  a  high  of  0.5082. 

The  momentary  creep  compliance  curves  were  shifted  to  form 
master  momentary  compliance  curves.  The  master  curve  data, 
referenced  to  the  81  hour  aging  time,  is  shown  in  figure  4.21. 
The  amount  of  horizontal  shift  -loga  and  the  vertical  shift 
logB  for  each  creep  curve  are  given  in  table  4.4. 

A  plot  of  the  aging  data,  —loga  versus  logt,,  used  to  calculate 
the  aging  shift  rate,  ^i,  for  each  test  sequence  are  shown  in 
figure  4.22.  The  aging  shift  rate,  i.e.,  slope  of  the  curve, 
increased  as  the  test  temperature  decreased.  A  comparison  of 
the  aging  shift  rates  determined  from  the  90  and  45  degree 
specimen  tests  is  presented  in  section  4.5. 

4.3.3  TTSP  Master  Momentary  Compliance  Curve 

TTSP  master  momentary  compliance  curves  were  formed  using 
the  same  methods  described  in  section  4.2.3  for  the  90° 
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specimens.  Again,  the  goal  was  to  obtain  the  compliance  curve 
at  the  long-term  test  reference  aging  time  of  9  hours  with  the 
temperature  at  338°F  (170°C)  . 

The  first  method  numerically  fit  the  test  data  from  various 
temperatures  at  an  aging  time  of  9  hours  to  Struik's  equation 
(eqn.  2.19)  using  the  average /n  value  of  0.4564.  The  horizontal 
shift  factors  were  calculated  using  equation  2.47.  The 
resulting  transverse  direction  TTSP  master  momentary  compliance 
curve,  referenced  to  338°F  (170°C)  and  /,  =  9  hours,  is  shown 
in  figure  4.23.  Note  that  fitted  curve  is  extrapolated  beyond 
the  test  data  at  times  greater  than  10^'^  seconds. 

The  TTSP  master  momentary  compliance  curve  resulting  from 
graphical  shifting  of  the  r,  =  9  hours  compliance  curves  is 
shown  in  figure  4.24,  along  with  the  required  shift  factors. 
Similar  to  the  90  degree  test  results,  the  numerical  and 
graphical  methods  produce  equivalent  shift  factors  for  the 
same  referenced  aging  time. 

Finally,  the  graphical  shift  procedure  was  used  to  obtain 
the  master  curve  referenced  to  338°F  (170°C)  at  an  aging  time 
of  81  hours  instead  of  9  hours.  The  resulting  curve,  referenced 
to  338°F  (170°C)  and  =  81  hours,  is  shown  in  figure  4.25. 
The  aging  shift  rate,  \i,  of  0.92  was  used  in  the  calculations 
to  determine  the  horizontal  shift  required  to  move  the  r,  =  81 
hours  curve  to  a  ?,  =  9  hours  reference  time.  The  shifted  curve 
is  shown  in  figure  4.26  along  with  the  curves  that  were 
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determined  by  the  two  other  methods.  Initially/  the  three 
curves  are  nearly  identical  initially.  However,  as  the  time 
increases  beyond  10^'^  seconds,  the  curves  start  to  separate. 
This  time  was  observed  to  correspond  to  the  end  of  the  shifted 
compliance  data  shown  in  figures  4.23-4.24.  The  separation 
of  the  curves  beyond  10^'^  seconds  results  from  the  extrapolation 
of  the  curves  beyond  the  test  data.  Obtaining  a  better 
compliance  representation  beyond  10^'^  seconds  would  require 
testing  at  higher  temperatures.  Unfortunately,  nonlinear 
effects,  aging  to  equilibrium,  and  possible  specimen  rupture 
may  occur  as  the  test  temperature  approaches  Tg. 

4.4  0  Degree  Specimen  Tests 

The  0  degree  specimen  orientation  was  used  to  determine 
the  fiber  direction  compliance,  S^,  and  the  fiber/transverse 
coupling  compliance,  5,2.  The  0  degree  orientation  is  dominated 
by  the  properties  of  the  elastic  graphite  fiber.  Previous 
testing  of  graphite  fiber  composites  [26,50,68,74]  revealed 
little  time-dependent  behavior  in  the  0  degree  fiber  orien¬ 
tation.  As  expected,  neither  5„  or  5,2  showed  appreciable 
time-dependent  behavior  at  temperatures  from  338°F  (170°C)  to 
399. 2°F  (204°C)  .  Thus,  both  5„  and  5,2  were  considered  temperature 
and  aging  independent . 
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4.4.1  5„  Determination 

The  0  degree  specimen  orientation  was  tested  at  338°F  (170°C) 
for  linear  behavior  at  stress  levels  of  5000  psi  (34.5  MPa) 
and  7500  psi  (51.7  MPa).  The  resulting  5„  compliance  curves 
were  nearly  identical  and  are  shown  in  figure  4.27.  The 
compliance  data  from  an  aging  sequence  test  conducted  at  338°F 
(170°C)  is  shown  in  figure  4.28.  No  physical  aging  effects 
were  evident  from  the  test  results.  The  data  points  from  the 
338°F  (170°C)  aging  tests  were  averaged  yielding  a  fiber 
compliance  value  of  5i,  =  40.230  (1/Gsi)  (0.005835  (1/GPa) )  . 
The  average  compliance  value  is  shown  as  a  solid  line  in  figure 
4.28. 

The  compliance  data  from  a  second  aging  sequence  test, 
conducted  at  399. 2°F  (204°C) ,  is  shown  in  figure  4.29.  A 
slight  negative  time  dependence  (contraction  under  tension) 
of  8  HE  over  a  10^  second  creep  test  was  noticed  at  the  higher 
temperature.  Griffith  et  al  [26]  also  noticed  a  slight  negative 
time  dependence  during  his  testing  of  0  degree  specimens  of 
T300/934  graphite/epoxy  composite  material.  He  attributed  the 
phenomena  to  the  viscoelastic  nature  of  the  strain  gage 
adhesive.  The  small  magnitude  of  the  strain  change,  only  8 
[i£,  could  also  be  attributed  to  normal  experimental  error 
associated  with  the  strain  measurement. 
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4.4.2  5)2  Determination 


The  fiber /transverse  coupling  compliance  term,  5,2,  was 

obtained  by  measuring  the  strain  in  the  transverse  direction 
on  the  0  degree  specimen,  as  explained  in  section  2.3.  The 
5i2  compliance  curves  from  the  linearity  check  on  the  0  degree 
specimen  are  shown  in  figure  4.30.  A  linear  response  was 
obtained  between  the  5000  psi  (34.5  MPa)  and  7500  psi  (51.7 
MPa)  applied  stress  levels.  Since  the  5i2  compliance  value  is 
negative,  the  log  of  the  absolute  value  of  the  S12  compliance 
term  was  plotted. 

The  compliance  data  from  an  aging  sequence  test,  conducted 
at  338°F  (170°C),  is  shown  in  figure  4.31.  Similar  to  5„,  the 
S12  compliance  term  showed  no  aging  effects  or  time-dependent 
behavior  at  338°F  (170°C)  .  An  average  compliance  value,  5,2  = 
-13.991  (1/Gsi)  (-0.002029  (1/GPa) ) ,  was  calculated  and  is 
shown  as  a  solid  line  in  figure  4.31.  The  compliance  data 
shown  in  figure  4.32  was  from  an  aging  sequence  test  conducted 
at  399. 2°F  (204°C)  .  The  apparent  time-dependence  of  the  5,2 
term  results  from  a  small  change  in  strain,  9  n£  over  the  10^ 
second  creep  test,  at  the  higher  test  temperature.  It  was 
impossible  to  determine  if  this  small  strain  change  was  a 
viscoelastic  material  response  or  a  strain  measurement  anomaly 
at  the  higher  test  temperature. 
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4.5  Temperature  Dependence  of  Aging  Shift  Rate 

The  aging  shift  rate,  Ji,  of  the  composite  material  was 
found  to  have  a  temperature  dependence  similar  to  that  of 
polymer  materials.  The  aging  shift  rate  values  calculated  at 
various  test  temperatures  with  specimen  fiber  orientations  of 
90,  45,  and  30  degrees  are  shown  in  figure  4.33.  Although 
some  scatter  exists  between  specimen  orientations,  a  general 
trend  in  the  data  appears  and  is  represented  by  the  solid  curve 
in  figure  4.33.  The  aging  shift  rate  appears  to  increase  from 
aging  equilibrium  (zero  aging  shift  rate)  at  Tg  to  a  plateau 
region  where  the  aging  shift  rate  is  about  unity.  Tests  at 
temperatures  below  338°F  (170°C)  were  not  conducted  because  of 
difficulties  in  measuring  the  small  creep  response.  It  is 
expected  that  the  aging  shift  rate  would  decrease  with  lower 
test  temperatures  as  the  material  approaches  the  secondary  |3 
transition  temperature.  Below  this  temperature  physical  aging 
ceases  because  the  free  volume  can  no  longer  change  as  explained 
in  section  2.2.1.  Sullivan  [66]  noticed  a  slight  decrease  in 
the  aging  shift  rate  of  a  composite  material,  made  of  Dow 
Derakane  470-36  resin  with  glass  fibers,  at  temperatures  of 
Tg  -  126°F  (Tg  -  70°C)  .  The  P  transition  temperature  is  usually 
determined  using  Dynamic  Mechanical  Analysis  (DMA) .  The  test 
is  conducted  over  a  wide  temperature  range  obtained  by  heating 
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the  specimen  in  an  oven  and  cooling  it  with  a  liquid  nitrogen 
injection  system.  The  ^  transition  temperature  was  not 
determined  during  this  study. 

The  decrease  in  |i  observed  at  338°F  (170°C)  in  the  90  degree 
specimen  was  not  found  in  the  45  or  30  degree  specimen 
orientations.  A  more  extensive  characterization  program  with 
multiple  specimen  tests  for  each  fiber  orientation  would  be 
required  to  determine  if  the  decrease  in  the  90  degree  p.  value 
at  338°F  (170°C)  was  a  point  anomaly  or  a  material  trend  for 
the  90  degree  specimen  orientation. 

4 . 6  Temperature  Dependence  of  Shift  Factors 

The  horizontal  shift  factors  required  to  form  a  TTSP  master 
momentary  compliance  curve  are  functions  of  both  the  test 
temperature  and  the  referenced  aging  time.  Both  the  Arrhenius 
and  WLF  models,  which  are  usually  used  to  represent  temperature 
shift  factor  data,  are  based  on  polymer  viscoelasticity  behavior 
without  physical  aging.  Any  future  models  developed  must  also 
include  the  aging  state  of  the  material  as  a  variable.  The 
shift  factors  computed  for  compliance  curves  obtained  at 
different  temperatures  with  the  same  aging  time  have  con¬ 
tributions  due  to  temperature  and  contributions  due  to  aging. 
The  aging  contribution  is  not  the  same  for  each  compliance 
curve  because  the  aging  shift  rate,  ^i,  has  a  temperature 
dependence . 
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The  horizontal  shift  factors,  — logaT,  that  are  used  to  form 

the  TTSP  master  momentary  compliance  curves  shown  in  figures 
4.13-4.14  and  4.24-4.25  are  plotted  versus  test  temperature 
in  figures  4.34  and  4.35  respectively.  An  observation  from 
these  figures  is  that  the  shift  factors  for  the  81  hour 
referenced  aging  time  appear  to  be  larger  in  magnitude  than 
for  the  9  hour  referenced  aging  time.  This  may  be  due  to  the 
decrease  in  the  aging  rate  at  temperatures  greater  than  the 
referenced  temperature  of  338°F  (170°C)  .  Continued  investi¬ 
gation  of  the  aging  dependence  of  the  horizontal  shift  factors 
offers  an  excellent  opportunity  for  follow  on  work  to  this 
study. 


4.7  Quench  Rate  Effect 

The  effect  of  physical  aging  has  been  shown  to  significantly 
change  the  compliance  properties  of  a  specimen  that  is  quenched 
from  above  Tg  to  an  aging  temperature  below  Tg.  The  high 
pressure  air  quenching  procedure  used  in  this  study  yielded 
a  quench  rate  of  9.5°F/sec  (5.3°C/sec)  .  This  rate  was  determined 
using  the  procedures  described  in  section  3.2.2.  However, 
normal  composite  processing  procedures  would  use  a  slower  rate 
to  cool  the  material.  An  obvious  question  is:  How  is  physical 
aging  process  affected  by  a  slower  cooling  rate? 

An  effort  to  study  the  cooling  rate  effect  on  physical 
aging  was  made  by  testing  the  90  and  45  degree  specimens  at 
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a  slower  cooling  rate.  A  slower  cooling  rate  was  obtained 
by  shutting  the  oven  controller  off  and  allowing  the  closed 
oven  to  cool  to  the  aging  temperature.  This  procedure  yielded 
an  oven  cooling  rate  of  0.056°F/sec  (0 . 031°C/sec) ,  or  about 
170  times  slower  than  the  air  quench.  Aging  sequence  tests 
were  conducted  on  both  the  90  and  45  degree  specimens  after 
they  were  oven  cooled  to  383°F  (195°C)  .  The  shifted  master 
compliance  curves,  referenced  to  t,  =  81  hours,  are  shown  in 
figures  4.36  and  4.37  for  the  90  and  45  degree  specimens 
respectively.  The  compliance  curves  are  essentially  the  same 
shape  with  the  oven  cooled  data  shifted  slightly  to  longer 
times.  The  aging  shift  rate,  for  the  compliance  direction 
decreased  from  0.66  for  the  air  quench  to  0.59  for  the  oven 
cool.  Similarly,  (i,  for  the  compliance  direction  decreased 
from  0.78  for  the  air  quench  to  0.73  for  the  oven  cool.  The 
decrease  in  aging  shift  rate  with  decreasing  cooling  rate  is 
consistent  with  the  observations  of  Booij  and  Palmen  [6]  during 
their  study  of  the  effects  of  thermal  history  on  the  viscoe¬ 
lasticity  of  ABS  samples. 

Clearly  the  effect  of  physical  aging  continues  at  the 
slower  oven  cooling  rate.  Theoretically  an  infinitely  slow 
cooling  rate  would  produce  a  state  of  free  volume  equilibrium 
where  physical  aging  would  not  be  possible.  However,  time 
constraints  associated  with  the  manufacture  of  composite  parts 
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require  removal  of  the  part  from  the  oven  as  soon  as  the  curing 
or  consolidation  process  is  complete.  Physical  aging  effects 
should  therefore  be  expected  in  actual  composite  structures. 
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5  LONG-TERM  CREEP  RESULTS  AND  DISCUSSION 


In  the  previous  chapter,  the  short-term  testing  used  to 
characterize  the  effect  of  physical  aging  on  the  compliance 
response  of  the  composite  material  in  each  of  four  material 
principal  directions  was  described.  In  this  chapter,  the 
results  of  long-term  tests,  conducted  to  verify  the  influence 
of  physical  aging  on  the  Radel  X/IM7  composite  material,  are 
presented.  Long-term  compliance  results  are  compared  with 
theoretical  predictions  determined  using  Struik's  effective 
time  theory  and  the  short-term  characterization  data. 

5.1  90  Degree  Specimen  Test,  S22  Comparison 

A  comparison  of  the  transverse  direction  long-term  com¬ 
pliance  predictions  with  the  observed  test  data  is  shown  in 
figure  5.1.  The  long-term  test  was  conducted  at  a  temperature 
of  338°F  (170°C)  with  a  starting  age  time  of  9  hours.  The  test 
data  is  depicted  with  asterisk  symbols,  while  the  predictions 
obtained  using  each  of  the  TTSP  master  momentary  compliance 
curves  from  figure  4.15,  along  with  the  effective  time  theory, 
are  shown  as  various  line  types.  All  of  the  long-term  pre¬ 
dictions  with  aging  show  the  same  shape  as  the  test  data.  The 
predictions  consistently  under-predicted  the  compliance  test 
data.  The  same  under-prediction  trend  was  noticed  in  long-term 
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compliance  data  presented  by  Sullivan  [65].  The  cause  of  the 
under-prediction  is  believed  to  be  related  to  drift  of  the 
strain  gage  voltages  during  long-term  tests.  The  effect  of 
strain  gage  drift  was  effectively  canceled  out  during  short-term 
tests  by  the  subtraction  of  the  recovery  voltages  from  the 
previous  test.  The  accumulation  of  the  drift  of  the  strain 
gage  voltages  could  be  responsible  for  the  increased  strain 
readings  during  long-term  tests.  Small  changes  in  the  bridge 
balance  and  excitation  voltages  were  noticed  at  the  end  of 
long-term  tests.  Unfortunately,  numerous  factors  can  effect 
the  drift  of  strain  gage  voltages  over  long  time  durations. 
The  factors  include  the  initial  balance  of  the  specimen,  quality 
of  the  solder  connections,  supply  voltage  variations,  and 
specimen  thermal  history.  Measurement  of  the  specimen  strain 
using  optical  techniques  along  with  strain  gages  offers  a 
possible  future  means  to  determine  the  strain  gage  drift  effect. 

The  sensitivity  of  the  long-term  compliance  predictions 
to  the  short-term  data  fitting  procedure  is  apparent  in  figure 
5.1.  The  variation  of  the  long-term  predictions  from  the 
experimental  data  range  from  7.3%  for  the  =  0.41671  curve 
to  3.5%  for  the  TTSP  curve  which  was  obtained  at  an  aging  time 
of  81  hours  and  shifted  to  the  referenced  aging  time  of  9 
hours.  The  closeness  of  the  predictions  makes  it  difficult 
to  determine  if  one  data  reduction  technique  is  more  reliable 
than  another. 
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Also  included  in  figure  5.1  is  the  compliance  prediction 
calculated  using  the  real  time  instead  of  effective  time  for 
the  TTSP  =  9  hr.  data  reduction  technique.  The  prediction 
without  aging  effects  included  over-predicts  the  compliance 
data  by  50%.  The  corresponding  prediction  with  aging  effects 
under-predicts  the  compliance  data  by  5.7%.  Clearly,  the 
effect  of  physical  aging  must  be  included  in  order  to  use 
short-term  momentary  compliance  data  to  predict  the  long-term 
compliance  response. 

5.2  45  Degree  Specimen  Tests  /  *^66  Comparison 

Long-term  shear  direction  compliance  tests  were  conducted 
at  a  temperature  of  338°F  (170°C)  at  starting  age  times  of  9 
and  3  hours.  As  in  figure  5.1,  the  test  data  is  again  depicted 
with  asterisk  symbols,  while  the  compliance  predictions  are 
shown  with  various  line  types  in  figures  5.2  and  5.3  for  the 
9  and  3  hour  starting  age  times  respectively.  Noted  with 
respect  to  both  figures  is  the  region  where  the  predictions 
were  made  using  the  extrapolation  of  the  characterization  data. 
The  sensitivity  of  the  long-term  prediction  to  the  short-term 
fitting  procedure  becomes  more  apparent  in  the  extrapolated 
region.  The  best  fit  to  the  test  data  was  the  prediction 
obtained  by  graphically  shifting  the  TTSP  data  to  the  aging 
time  of  the  test.  The  predictions  from  this  method  were  only 
6.3%  and  7.3%  from  the  test  data  at  starting  age  times  of  9 
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and  3  hours  respectively. 

The  effect  of  physical  aging  on  the  long-term  compliance 
prediction  becomes  evident  at  smaller  starting  aging  times. 
The  compliance  prediction  without  aging  for  the  9  hour  starting 
age  time,  shown  in  figure  5.2,  is  65%  from  the  actual  test 
data  at  the  end  of  the  test.  The  prediction  without  aging  for 
a  test  at  a  3  hour  starting  age  time,  shown  in  figure  5.3,  is 
greater  than  165%  from  the  actual  test  data  at  the  time 
corresponding  to  the  end  of  the  9  hour  aging  time  test.  As 
the  aging  time  before  the  start  of  the  creep  test  was  increased, 
the  apparent  effect  of  physical  aging  was  shifted  to  longer 
test  times.  Significant  variation  between  "with  aging"  and 
"without  aging"  predictions  starts  at  about  10^‘^  and  10^  seconds 
respectively  for  3  and  9  hours  aging  start  times.  Clearly, 
the  aging  time  at  the  start  of  the  creep  test  and  the  test 
duration  determine  if  physical  aging  significantly  changes  the 
creep  response  of  the  material  during  the  test. 

The  ambiguous  thermal  history  of  the  specimens  used  in 
previous  test  programs  [13,25,26,68]  makes  it  impossible  to 
determine  if  physical  aging  effects  were  significant.  Typi¬ 
cally,  specimens  from  other  test  programs  were  tested  in 
succession  at  increasing  stress  levels  at  a  constant  temperature 
or  at  a  constant  stress  level  with  the  temperature  increased 
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between  tests.  Since  the  tests  were  usually  conducted  in 
succession,  the  physical  aging  effects  were  most  likely  blended 
in  with  temperature  and  stress  effects. 

5.3  30  Degrao  Specimen  Test,  S„  Comparison 

The  load  direction  compliance,  referred  to  as  S„,  was 

determined  using  equation  2.68a  and  the  four  principal  com¬ 
pliance  terms  5n, 5i2, 522*  and  The  TTSP  momentary  compliance 
curves  used  in  the  computations  were  determined  using  the 
graphical  procedure  at  an  aging  time  of  9  hours  and  referenced 
to  a  temperature  of  338°F  (170°C)  .  Aging  shift  rate,  \x.,  values 
of  0.77  and  0.93  were  substituted,  respectively  for  the  ^22  and 
S«  principal  compliance  directions,  into  equation  2.55  during 
the  calculation  of  the  effective  time  for  long-term  predictions. 
The  test  data  and  predictions  for  the  load  direction  long-term 
compliance  at  a  temperature  of  338°F  (170°C)  and  a  starting 
age  time  of  9  hours  is  shown  in  figure  5.4.  Predictions  are 
shown  for  fiber  angles  of  29.5,  30  and  30.5  degrees  to  dem¬ 
onstrate  the  sensitivity  of  the  compliance  predictions  to  the 
possible  misalignment  of  the  fiber  orientation  during  specimen 
fabrication  or  the  misalignment  of  the  strain  gage.  The  29.5 
degree  prediction  matches  the  initial  compliance  value  the 
best,  but  under-predicts  the  compliance  at  10®  seconds  by  9.1%. 
This  under-prediction  of  the  compliance  was  also  seen  in  the 
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522  and  5^^  long-term  predictions.  The  compliance  prediction 
obtained  without  physical  aging  effects  over-predicted  the 
test  compliance  value  at  10*  seconds  by  137%. 

5.4  Overall  S„  Comparison 

Predictions  of  the  load  direction  compliance,  s„, 

for  fiber  orientation  angles  ranging  from  0  to  90  degrees  at 
a  temperature  of  338°F  (170°C)  and  an  aging  start  time  of  9 
hours  are  shown  in  figure  5.5.  The  90,  45,  30  and  0  degree 
fiber  angle  long-term  load  direction  compliance  test  data  is 
also  shown  in  figure  5.5  with  asterisk  symbols.  The  gradual 
stiffening  of  the  material  from  physical  aging  is  apparent  in 
both  the  predictions  and  the  test  data.  The  predictions  and 
test  data  compare  quite  well  considering  the  possible  sources 
of  experimental  error  and  specimen  to  specimen  material  property 
variations.  Ideally,  multiple  specimens  of  each  fiber 
orientation  would  be  tested  so  a  statistical  distribution  of 
the  compliance  response  could  be  determined. 

The  shape  of  the  compliance  curve  beyond  the  end  of  the 
test  data  at  10®  seconds  can  be  estimated  by  determining  the 
double-logarithmic  slope  of  the  long-term  creep  curves.  The 
shape  of  the  double-logarithmic  slope  primarily  depends  on  the 
aging  shift  rate  |1.  A  p.  value  of  0.85  yields  a  nearly  constant 
double-logarithmic  slope  value  producing  a  straight  line  on 
the  double-logarithmic  plot  of  compliance  versus  time.  Values 
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of  p,  greater  than  0.85  produce  a  slope  that  slightly  decreases 
with  time.  When  the  value  of  |i  is  less  than  0.85  the  slope 
slightly  increases  with  time.  The  proximity  of  the  ^  values 
for  both  the  S22  and  ^66  compliance  directions  to  0.85  explains 
the  nearly  straight  line  prediction  of  the  compliance  at  longer 
test  times. 
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6  SUMMARY,  RECQMMEMDATIONS,  AND  CONCLUSIONS 


€ .  1  Summary 

The  goal  of  this  study  was  to  determine  the  effect  of 
physical  aging  on  the  creep  response  of  Radel  X/IM7  composite 
material  in  the  linear  viscoelastic  region.  Three  program 
objectives  were  defined  to  meet  this  goal.  The  first  objective 
was  to  develop  the  experimental  procedures  to  conduct  physical 
aging  experiments.  The  second  objective  was  to  apply  the 
experimental  procedures  to  characterize  the  influence  of  sub 
Tg  physical  aging  on  the  creep  response  of  the  thermoplastic 
composite  material  in  each  of  the  material  principal  directions. 
The  third  objective  was  to  use  short-term  momentary  compliance 
test  results  and  the  effective  time  theory  to  predict  the 
long-term  compliance  response  of  unidirectional  laminates. 

All  three  objectives  were  accomplished.  Experimental 
problems  were  encountered  and  solved  during  the  study. 
Momentary  compliance  results  were  used  to  characterize  the 
aging  effect  by  determining  the  aging  shift  rate  of  the  material. 
Momentary  creep  tests  revealed  significant  changes  in  the 
compliance  response  of  the  composite  as  it  physically  aged. 
Finally,  long-term  predictions  using  the  effective  time  theory 
varied  from  the  test  data  by  less  than  10%  at  the  end  of  the 
long-term  test  period. 
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The  Radel  X/IM7  thermoplastic  composite  material  charac¬ 
terization  process  was  completed  by  conducting  momentary  creep 
tests  of  90,  45,  and  0  degree  unidirectional  laminates.  The 
effect  of  physical  aging  on  the  compliance  response  was 
determined  in  each  of  the  four  principal  material  directions. 
The  resin  dominated  transverse  direction  compliance,  and 
shear  direction  compliance,  S«6,  were  found  to  be  time  and  aging 
dependent.  The  fiber  direction  compliance,  S,i,  and  the 
fiber/transverse  coupling  compliance,  5,2,  were  found  to  be 
virtually  independent  of  time  and  aging  effects. 

The  aging  shift  rate,  |i,  of  the  composite  was  found  to 
have  a  temperature  dependence.  The  aging  shift  rate  appeared 
to  increase  from  a  predicted  zero  value  at  the  glass  transition 
temperature  to  a  plateau  region  starting  at  a  temperature  about 
63°F  (35°C)  below  T,.  It  was  not  possible  during  this  study 
to  determine  if  the  variation  between  the  aging  shift  rate  of 
the  various  specimen  orientations  was  experimental  scatter  or 
a  function  of  the  fiber  angle  of  the  specimen. 

The  Time-Temperature  Superposition  Principle  was  suc¬ 
cessfully  used  to  determine  TTSP  master  momentary  compliance 
curves  at  a  specific  referenced  aging  time.  The  horizontal 
shift  factors,  -loga-r,  were  found  to  depend  on  temperature  and 
the  referenced  aging  time  of  the  master  curve.  Long-term 
predictions  using  the  TTSP  master  momentary  compliance  curves. 
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without  including  effective  time  theory  to  account  for  the 
physical  aging  process,  significantly  overestimated  the  com¬ 
pliance  response. 

The  rate  the  specimen  was  cooled  from  above  T,  to  the  aging 
temperature  was  found  to  slightly  affect  the  aging  shift  rate. 
The  aging  shift  rate  was  found  to  decrease  when  the  specimen 
cooling  procedure  was  changed  from  an  air  quench  to  an  oven 
cool  method.  Continually  decreasing  the  cooling  rate  to  an 
infinitely  slow  rate  should  theoretically  produce  a  state  of 
equilibrium  free  volume  that  would  eliminate  physical  aging 
effects  in  the  material.  Although,  achieving  the  equilibrium 
free  volume  may  require  a  cooling  rate  that  is  impractical  for 
industrial  composite  manufacturing  applications 

Momentary  compliance  data  and  the  effective  time  theory 
were  used  to  accurately  predict  long-term  S22  and  compliance 
responses  of  the  90  and  45  degree  specimens.  Additionally, 
the  four  principal  compliance  terms  were  substituted  into  a 
laminate  transformation  equation  to  predict  the  load  direction 
compliance  for  a  30  degree  specimen.  The  long-term  load 
direction  compliance  prediction  matched  the  test  data  within 
normal  error  ranges  associated  with  slight  misalignments  in 
the  fiber  orientation  or  strain  gage  direction. 

6 . 2  Recommendations 

This  research  project  successfully  addressed  some  of  the 
questions  concerning  the  effect  of  physical  on  composite 


129 


materials  and  subsequently  created  a  few  new  ones.  Many  of 
the  uncertainties  in  this  study  came  from  using  test  results 
from  one  specimen.  Ideally,  each  creep  test  would  be  repeated 
with  multiple  specimens  so  that  a  reliable  statistical  dis¬ 
tribution  of  the  compliance  properties  could  be  determined. 
Obviously,  testing  more  specimens  would  require  a  much  longer 
time  period  and/or  the  use  of  multiple  test  fixtures.  Future 
studies  should  concentrate  on  one  fiber  angle  orientation  and 
examine  numerous  test  variables  in  more  detail. 

During  this  study,  the  effect  of  moisture  was  minimized 
by  drying  the  material  and  keeping  it  in  an  environment  with 
desiccant  to  maintain  constant  humidity.  The  absorption  and 
desorption  of  moisture  will  affect  the  compliance  response  of 
the  composite  in  a  manner  similar  to  temperature.  Both  the 
moisture  and  temperature  environment  will  experience  cyclic 
variation  during  the  design  life  of  the  composite  structure. 
The  effect  of  physical  aging  on  the  properties  of  the  material 
after  complex  moisture,  thermal,  and  stress  histories  needs 
further  investigation. 

The  physical  aging  theory  should  be  extended  so  that 
nonlinear  stress  effects  in  composite  materials  are  included. 
Sequencing  of  creep  tests  is  no  longer  possible  in  the  nonlinear 
stress  range.  Thus,  the  specimen  must  be  reheated  above  Tg 
and  quenched  back  to  the  aging  temperature  between  each  creep 
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test.  The  study  of  nonlinear  effects  will  significantly 
increase  the  amount  of  test  time  required  to  characterize  the 
material . 

Development  of  improved  strain  measuring  methods  which  do 
not  require  elaborate  time-consuming  specimen  preparation  are 
needed  in  order  to  conduct  a  larger  quantity  of  tests  in  a 
short-time  frame.  Strain  measurement  procedures  utilizing 
laser  interferometry  may  provide  required  long-term  stability 
and  ease  of  use  required  to  measure  the  strains  in  a  high 
temperature  environment.  Another  alternative  is  to  use  dynamic 
mechanical  analyzer  (DMA)  equipment  to  measure  the  creep 
response.  The  strain  or  displacement  is  measured  using  a 
built-in  sensor  in  the  equipment.  This  alleviates  problems 
associated  with  strain  gages.  DMA  equipment  usually  includes 
a  liquid-nitrogen  cooling  attachment  to  the  oven.  This  feature 
can  be  used  to  produce  various  cooling  rates  in  the  oven.  This 
would  facilitate  the  study  of  cooling  rate  effects  on  physical 
aging.  Unfortunately,  DMA  equipment  is  quite  expensive  and 
is  rarely  devoted  to  long-term  aging  studies. 

Finally,  the  prediction  capability  needs  to  be  extended 
from  unidirectional  laminates  to  general  laminates  and  verified 
with  experimental  testing.  This  could  be  accomplished  by 
modifying  incremental  numerical  procedures  developed  by  Dillard 
et  al  [13]  or  Gramoll  et  al  [25]  to  include  physical  aging 
effects . 
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6 . 3  Conclusions 


The  creep  response  of  Radel  X/IM7  composite  material  has 
been  found  to  be  significantly  affected  by  physical  aging  in 
the  temperature  range  between  338°F  (170°C)  and  399. 2°F  (204°C)  . 
Physical  aging  effects  become  more  predominant  during  long-term 
tests  where  the  test  duration  is  much  longer  than  the  aging 
time  at  the  start  of  the  test.  Prediction  methods  based  solely 
on  the  time-temperature  superposition  principle  will  signif¬ 
icantly  over-predict  the  creep  response  of  the  material  if  the 
effect  of  physical  aging  is  ignored. 
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Table  4.1  Comparison  of  Commercial  Data  of  Radel  X/T650-42 
and  Radel  X/IM7  Composite  Materials. 


Property 

Radel  X/T650-42* 

(62%  fiber  volume) 

Radel  X/IM7 

(test  data) 

Tensile  Strength  90° 

6.91  Ksi 

(47.6  MPa) 

4.46  Ksi 

(30.8  MPa) 

Tensile  Modulus  90° 

1.17  Msi 

(8.07  GPa) 

1.07  Msi 

(7.38  GPa) 

Tensile  Strength  0° 

330  Ksi 

(2275  MPa) 

261  Ksi 

(1800  MPa) 

Tensile  Modulus  0° 

23.7  Msi 

(163.4  GPa) 

23.1  Msi 

(159.3  GPa) 

®  Amoco  Performance  Products 


142 


Table  4.2  Comparison  of  Commercial  T650-42  and  IM7  Fiber 

Properties. 


1  Property 

T650-42'* 

Fiber 

IM7‘’  1 

Fiber 

Tensile  Strength 

730  Ksi 

(5037  MPa) 

683  Ksi 

(4713  MPa) 

Tensile  Modulus 

42  Msi 

(290  GPa) 

41  Msi 

(283  GPa) 

Density 

0.064  Ib/in^ 

(1.78  g/cm^) 

0.064  Ib/in^ 

(1.78  g/cm^) 

" " ' 

1  Ultimate 

1  Elongation 

1.7% 

1.6% 

®  Thornel,  Advance  Composite  Systems 
Hercules,  Graphite  Materials  &  Composite  Structures  Marketing 
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Table  4.3  Shift  Factors  Required  for  Master  Momentary 
Compliance  Curves  Referenced  to  =  81  hours. 


Temperature 

(hour) 

-log  a 

logB 

399°F  (204°C) 

1 

-0.7914 

0.0168 

3 

-0.5713 

0.0095 

9 

-0.3584 

-0.0016 

27 

-0.1306 

-0.0106 

81 

0.0000 

0.0000 

392°F  (200°C) 

-0.9535 

-0.0015 

3 

-0.6805 

0.0000 

9 

-0.4183 

0.0009 

27 

-0.1968 

0.0024 

81 

0.0000 

0.0000 

383°F  (195°C) 

1 

-1.2458 

-0.0124 

3 

-0.9400 

-0.0090 

9 

-0.6105 

-0.0067 

27 

-0.3043 

-0.0066 

81 

0.0000 

0.0000 

365°F  (185°C) 

1 

-1.5795 

-0.0071 

3 

-1.1668 

-0.0065 

9 

-0.8030 

-0.0051 

27 

-0.3647 

-0.0036 

81 

0.0000 

0.0000 

338°F  (170°C) 

1 

-1.4950 

-0.0142 

3 

-0.9979 

-0.0121 

9 

-0.7098 

-0.0074 

27 

-0.3158 

-0.0037 

81 

0.0000 

0.0000 
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Table  4.4  Shift  Factors  Required  for  Master  Momentary 
Compliance  Curves  Referenced  to  =  81  hours. 


Temperature 

t,  (hour) 

-log  a 

logB 

|  383°F  (195°C) 

1 

-1.4755 

-0.0041 

3 

-1.0519 

-0.0119 

9 

-0.6838 

-0.0083 

27 

-0.2987 

-0.0038 

81 

0.0000 

0.0000 

365°F  (185°C) 

1 

-1.6068 

-0.0243 

3 

-1.1803 

-0.0200 

9 

-0.7435 

-0.0181 

27 

-0.3220 

-0.0129 

81 

0.0000 

0.0000 

338°F  (170°C) 

1 

-1.7643 

0.0021 

3 

-1.2829 

-0.0052 

9 

-0.8979 

-0.0021 

27 

-0.4168 

0.0013 

81 

0.0000 

0.0000 
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Figure  2.3  Stress  and  Strain  Response  for  a 
Constant  Strain  Rate  Test. 
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6  Representation  of  Maxwell  and  Kelvin 
Mechanical  Analog  Models. 
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Log  Time 

Figure  2.8  Time-Temperature  Superposition  Principle  Illustration. 


Figure  2.9  Visualization  of  Two-Dimensional 

Hard- Sphere  Model  of  Free  Volume  [75] . 
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Figure  2.11  Typical  Compliance  Response  of  Thermoplastic  and  Thermoset  Materials. 


Figure  2.12  The  Origin  of  Aging  Explained 

from  Free  Volume  Concepts  [59] . 
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►  tog  t^,days 


NOTE:  Log  a  is  the  shift  relative  to  the  creep  curve 
measured  at  a  t^of  1  day. 


Figure  2.14  Calculation  of  |i  from  the  Shifting  of  the 
Creep  Curves  of  Rigid  PVC  Quenched  from 
90  °C  to  0,  20,  and  40  °C  [59] . 
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Figure  2.15  The  Aging  Shift  Rate,  p.  ,  vs .  Temperature 
for  Various  Polymers  [59] . 
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llustration  of  the  Sequence  of  Creep  and  Recovery  Tests 
or  Determining  Aging  Effects  [59] . 
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Figure  2.17  Illustration  of  Aging  Shift  Rate  Determination  by 
Shifting  Momentary  Creep  Compliance  Curves. 
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Figure  2.18  Illustration  of  Effective  Time  Variables. 


gure  2.19  Coordinate  Systems  Used  to  Define 
Material  Properties. 
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Figure  2.20  Strain  Gage  Rosette  Placement  on  45  degree  Specimen. 
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Figure  3.2  Illustration  of  Specimen  Attachment  to  Load  Frame. 
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Figure  3.-  Stabilization  Temperature  vs.  Quench  Time. 


69 
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Figure  3.4  Specimen  Temperature  vs.  Time  Data  Used  for 
Quench  Rate  Calculation. 
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VOLTAGE 

Figure  3.5  Calibration  of  Output  Voltage  vs.  Specimen  Load. 
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Figure  3.7  Hot  Press  Cure  Cycle. 
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Sample;  HADEL  X/IM7  T  P  A  HASTIER. 30 

Size:  24.4766  mg  I  /A  Operator:  paul  s  vail 

Method:  lOC/min  to  800  Run  Date:  2B-Sep-90  10: 41 

Comment:  Purged  in  air 
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Figure  3.10  Illustration  of  Full  Bridge  Wheastone  Circuit. 
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Figure  3.11  Strain  Gage  Curing  Cycles. 
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Figure  3.12  Initial  Zero  Drift  Measurements. 
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Figure  3.13  Zero  Drift  Measurements  after  Resoldering. 
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Figure  3.14  Zero  Drift  Measurements  after  Remounting  Dummy  Specimen. 
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Figure  3.15  Comparison  of  Line  Noise  Between  the  Vishay  2100  and  2210 
Signal  Conditioner  Amplifiers. 
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Figure  4.1  Stress  Strain  Curve  for  Transverse  Direction. 
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Figure  4.2  Stress  Strain  Curve  for  Fiber  Direction. 
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Figure  4.3  Linearity  Check  of  Transverse  Direction  Compliance  at  a  Constant 

Temperature  of  392  F  (200  C)  and  te  =  3  hours. 


TRANSVERSE  DIRECTION  COMPLIANCE.  S22  -  -0.75 

TEMPERATURE  =  338  F  (1 70  C).  te  =  9  hr.,  400  psi  (2.76  MPa) 


(BdO/l.)  30NVndl/\l00  001 


m 

IT) 

00 

O) 

O) 

C) 

o 

d 

>1— 

lO 

CO 

lO 

o> 

in 

GO 

o 

O) 

• 

cJ 

CO 

csi 

CO 

CM 

cvi 

(IS0/1-)  aoNvndi^oo  ooi 


186 


LOG  TIME  (sec) 

Figure  4.4  Check  of  Test  Sequencing  Effect  on  Transverse  Direction  Compliance 

at  a  Temperature  of  338  F  (170  C)  and  te  =  9  hours. 
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Figure  4.5  Transverse  Direction  Compliance  for  a  Constant  Temperature 

of  338  F  (170  C)  at  Various  Aging  Times. 
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Figure  4.6  Transverse  Direction  Compliance  for  a  Constant  Temperature 

of  365  F  (185  C)  at  Various  Aging  Times. 
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Figure  4.7  Transverse  Direction  Compliance  for  a  Constant  Temperature 

of  383  F  (195  C)  at  Various  Aging  Times. 
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Figure  4.8  Transverse  Direction  Compliance  for  a  Constant  Temperature 

of  392  F  (200  C)  at  Various  Aging  Times. 
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Figure  4.11  Aging  Shift  Rate  for  the  90  degree  Test 

Sequence  at  Various  Temperatures. 
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Figure  4.12  Transverse  Direction  TTSP  Master  Momentary  Compliance  Curve 
Computer  Fit  with  m  avg  =  0.4167  and  Referenced  to  T=338F  (170C) ,  te=9  hours. 
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Figure  4.13  Transverse  Direction  TTSP  Master  Momentary  Compliance  Curve 
Graphically  Fit  and  Referenced  to  T=338F  (170C),  te=9  hours. 
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Figure  4.15  Comparison  of  Transverse  Direction  TTSP  Master  Momentary 
Compliance  Curves  Referenced  to  T=338F  (170C),  te=9  hours. 
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Figure  4.16  Linearity  Check  of  Shear  Direction  Compliance  at  a  Constant 

Temperature  of  365  F  (185  C)  and  te  =  3  hours. 
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Figure  4.17  Check  of  Test  Sequencing  Effect  on  Shear  Direction  Compliance. 
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Figure  4.18  Shear  Direction  Compliance  for  a  Constant  Temperature 

of  338  F  (170  C)  at  Various  Aging  Times. 


(BdO/l-)  30NVndl/MO0  001 


lO 

in 

d 

I 


CD 

d 

I 


in 

CO 

d 

I 


o 

I 


in 

d 

I 


GO 

d 

I 


(ISO/I.)  30NVndl/MO0  001 


in 


CO 


o 

CD 

>52- 

LU 


h- 

(D 

O 


CM 


201 


Figure  4.19  Shear  Direction  Compliance  for  a  Constant  Temperature 

of  365  F  (185  C)  at  Various  Aging  Times, 
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Figure  4.20  Shear  Direction  Compliance  for  a  Constant  Temperature 

of  383  F  (195  C)  at  Various  Aging  Times. 
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Figure  4.21  Shear  Direction  Master  Momentary  Compliance  Curves 
at  Various  Temperatures  Shifted  to  the  81  hour  Aging  Time. 


Figure  4.22  Aging  Shift  Rate  for  the  45  degree  Test 

Sequence  at  Various  Temperatures. 


2 


SHEAR  DIRECTION  COMPLIANCE.  S66 


205 


Figure  4.23  Shear  Direction  TTSP  Master  Momentary  Compliance  Curve 
Computer  Fit  with  m  avg  =  0.4564  and  Referenced  to  T=338F  (170C),  te=9  hours 
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Figure  4.24  Shear  Direction  TTSP  Master  Momentary  Compliance  Curve 

Graphically  Fit  and  Referenced  to  T=338F  (170C),  te=9  hours 
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Figure  4.25  Shear  Direction  TTSP  Master  Momentary  Compliance  Curve 

Graphically  Fit  and  Referenced  to  T=338F  (170C),  te=81  hours. 
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Figure  4.26  Comparison  of  Shear  Direction  TTSP  Master  Momentary 

Compliance  Curves  Referenced  to  T=338F  (170C),  te=9  hours. 
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Figure  4.27  Linearity  Check  of  Fiber  Direction  Compliance  at  a  Constant 

Temperature  of  338  F  (170  C)  and  te  =  3  hours. 
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Figure  4.28  Fiber  Direction  Compliance  for  a  Constant  Temperature 

of  338  F  (170  C)  at  Various  Aging  Times. 


211 


LOG  TIME  (sec) 

Figure  4.29  Fiber  Direction  Compliance  for  a  Constant  Temperature 

of  399.2  F  (204  C)  at  Various  Aging  Times. 
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Figure  4.30  Linearity  Check  of  Fiber/Transverse  Coupling  Compliance  at 
Constant  Temperature  of  338  F  (170  C)  and  te  =  3  hours. 
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Figure  4,31  Fiber/Transverse  Coupling  Compliance  for  a  Constant  Temperature 

of  338  F  (170  C)  at  Various  Aging  Times. 
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Figure  4.32  Fiber/Transverse  Direction  Compliance  for  a  Constant  Temperature 

of  399.2  F  (204  C)  at  Various  Aging  Times. 
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Figure  4.33  Temperature  Dependence  of  the  Aging  Shift  Rate. 
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Figure  4.34  Temperature  Dependence  of  the  Transverse  Direction 

TTSP  Horizontal  Shift  Factors. 
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Figure  4.35  Temperature  Dependence  of  the  Shear  Direction 

TTSP  Horizontal  Shift  Factors. 
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Figure  4.36  Comparison  of  Air  Quenched  and  Oven  Cooled  Transverse  Direction 
Master  Momentary  Compliance  Curves  at  383  F  (195  C)  with  te=81  hours. 
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TRANSVERSE  DIRECTION  COMPLIANCE.  S22 
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Figure  5.1  Transverse  Direction  Long-Term  Compliance  Predictions  and 

Test  Data  at  338  F  (170  C)  and  te  =  9  hours. 
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Figure  5.2  Shear  Direction  Long-Term  Compliance  Predictions  and 
Test  Data  at  338  F  (170  C)  and  te  =  9  hours. 
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Figure  5.3  Shear  Direction  Long-Term  Compliance  Predictions  and 
Test  Data  at  338  F  (170  C)  and  te  =  3  hours. 


Figure  5.4  Load  Direction  Long-Term  Compliance  Predictions  for  the 
30  degree  Specimen  at  338  F  (170  C)  and  te  =  9  hours. 
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Figure  5.5  Load  Direction  Long-Term  Compliance  Predictions  and 
Test  Data  at  338  F  (170  C)  and  te  =  9  hours. 
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